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Abstract 

This manuscript is the outcome of the subgroup "PDFs, shadowing and pA 
collisions" from the CERN workshop "Hard Probes in Heavy Ion Collisions at 
the LHC". In addition to the experimental parameters for pA collisions at the 
LHC, the issues discussed are factorization in nuclear collisions, nuclear parton 
distributions (nPDFs), hard probes as the benchmark tests of factorization in 
pA collisions at the LHC, and semi-hard probes as observables with potentially 
large nuclear effects. Also, novel QCD phenomena in pA collisions at the LHC 
are considered. The importance of the pA program at the LHC is emphasized. 
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1. WHY pA COLLISIONS : INTRODUCTION AND SUMMARY 

Kari J. Eskola, Walter Geist, and Jianwei Qiu 

Hard Probes [1,2] are high-energy probes of the Quark Gluon Plasma (QGP) which are produced 
in the primary partonic collisions. The production of hard probes involves a large transfer of energy- 
momentum at a scale Q » Aqcd- Such hard probes include the production of Drell-Yan dileptons, 
massive gauge bosons, heavy quarks, prompt photons, and high-py partons observed as jets and high-p^ 
hadrons. In pp collisions, it is well established that the inclusive cross sections of these hard processes 
can be computed through collinear factorization, i.e. using short-distance cross sections of parton-parton 
scatterings and well-defined universal parton distribution functions (PDFs). While the partonic sub- 
cross sections and the scale evolution of the PDFs are calculable in perturbative QCD (pQCD), the PDFs 
contain nonperturbative information which must be extracted from the measured cross sections of various 
hard processes. 

Provided that leading-power collinear factorization is applicable also in AA collisions, the cross 
sections of hard probes can be used as benchmark cross sections against which the signals and properties 
of the QGP can be extracted. Therefore, it is of extreme importance that the applicability of factorization 
be tested in pA interactions, especially at the LHC where truly hard probes will finally become available 
for nuclear collisions. 

In this document, we divide the hard probes into two classes according to their scales: hard probes 
and semi-hard probes. The hard probes refer to observables with hard partonic subprocesses of energy 
exchange Q greater than several tens of GeV, up to 100 GeV and higher. For these hard probes, the 
process-dependent nuclear effects (i.e., power corrections) should remain negligible. The semi-hard 
probes correspond to probes at moderate scales, Q of few GeV up to 10. . . 20 GeV, where the process- 
dependent nuclear effects may already be sizable in pA collisions and even larger in AA collisions. This 
document mainly addresses the hard and semi-hard probes in pA collisions from the factorization (large 
scale) point of view but we also discuss (Sec. 7) other interesting new physics which can be explored in 
pA collisions at the LHC. 

Section 2 summarizes the main machine parameters for the pA runs at the LHC as well as the 
acceptance regions of the ALICE, CMS and ATLAS detectors. In addition to protons, the LHC machine 
can accelerate light beams as dilute as the deuteron and as compact as the helium nucleus. The machine 
can also accelerate a variety of heavier nuclei ranging from oxygen to lead. All three complementary 
detectors, ALICE, CMS and ATLAS, have wide rapidity coverage for multiplicity measurements, the ca- 
pability of precision centrality determination in pA collisions and excellent jet reconstruction efficiency 
for a range of transverse momentum, pr, and pseudorapidity, 77. 

Experimental data become more valuable and powerful when there are theoretical calculations to 
compare with. Because of the nature of the strong interactions, our ability to do calculations and to make 
predictions for hard probes in hadronic collisions completely relies on factorization in QCD perturbation 
theory. For the hard probes, the benchmark tests, factorization can be expected to hold in pp, pA and 
A A collisions. For the semi-hard probes which are sensitive to the properties of nuclear matter via 
multiple parton scatterings, there is, however, a critical difference between pA and AA collisions: while 
factorization may be valid for such semi-hard observables in pA collisions, it can fail in AA collisions. 
Factorization in nuclear collisions is summarized in Section 3. 

Nuclear parton distribution functions (nPDFs) are an essential ingredient in understanding the 
magnitude of the nuclear effects on the factorized hard probe cross sections. While AA collisions typi- 
cally are too complex for detailed verification of the nPDF effects on the semi-hard probe cross sections, 
pA collisions provide a much cleaner environment for studies of the nPDFs. The pA collisions can also 
probe the spatial dependence of the nPDFs - an important issue in the studies of hard probes in different 
AA centrality classes as well as in understanding the origin of the nuclear effects (shadowing) on the 
nPDFs. We devote Section 4 to discussion of several aspects of the nPDFs, ranging from global DGLAP 



fits to studies of the origin of shadowing and diffractive phenomena, and PDFs in the gluon saturation 
region at small values of x and Q. 

Section 5, "The benchmark tests in pA", deals with hard probes for which the nonfactorizable 
nuclear effects, i.e., process-dependent power corrections, should remain negligible. For these large-Q 
processes, the factorizable nuclear effects (nuclear effects on the partem distributions) also remain small. 
We emphasize that it is very important to measure these processes in pA collisions at the LHC: universal 
nuclear effects that are internal to the nucleus would, if observed, confirm the practical applicability of 
factorization in hard nuclear processes and would establish a reliable reference baseline for AA. It would 
be especially important to verify at what scale the observed nuclear dependence becomes universal and 
is completely determined by that of the nPDFs. On the other hand, observation of any significant nuclear 
effects in pA would be a truly interesting and important signal of sizable power corrections in the cross 
sections and thereby also of new multiple collision dynamics of QCD. 

For nuclear collisions, the very large-scale processes with Q > 100 GeV are in practice only 
available at the LHC, the ultimate hard probe machine. As shown in Table 1, due to the high cms-energy 
and luminosity - higher than in A A - the counting rates of the hard probes at the scale Q ~ 100 GeV 
in pA collisions will be large enough to give sufficient statistics for reliable measurements of the nuclear 
effects. For example, in a month (assumed to be 10 6 s) we expect about 3.3 • 10 8 jets with Et > 100 GeV 
and \r]\ < 2.5 to be collected in pPb collisions at yfs = 8.8 GeV with the maximum machine luminosity 
1.4 • 10 30 cm _2 s _1 and 2.6 • 10 7 jets with the maximum ALICE luminosity 1.1 • 10 29 cm _2 s _1 (for the 
luminosities, see Sec. 2.1). About half of these jets fall into the interval 100 GeV < E T < 120 GeV. The 
corresponding number of Z° bosons at \yz\ < 2.4 is 6.8 • 10 6 in a month with the maximum machine 
luminosity and 5.4 • 10 5 with the maximum ALICE luminosity. 

Table 1 : Expected counting rates (in 1/s, with two significant digits) of jets with transverse energy Et > 100 GeV and rapidity 
\r]\ < 2.5 (see Sec. 5.3), and heavy gauge bosons at rapidities \y\ < 2.4 in pPb collisions at ^fs = 8.8 TeV (see Sec. 5.1). 
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Section 6, "Processes with potentially large nuclear effects", discusses the physics of the semi- 
hard probes. The theoretical reference cross sections for these processes inevitably contain some model 
dependence, causing unknown theoretical uncertainties in the corresponding cross sections in AA col- 
lisions. For these semi-hard probes, it is even more important to experimentally confirm the reference 
cross sections in pA collisions because of the substantial non-universal nuclear dependence. That is, the 
reliable reference cross sections for semi-hard probes in AA can be obtained only if the corresponding 
pA measurements are performed. A concrete example of this is the clear suppression of hadrons with 
Pt ^ 10 GeV recently observed at RHIC in central AuAu collisions at yfs = 200 GeV relative to the 
pp data [3-5]. Only with help of the p(d)A data is it possible to verify that the origin of the suppression 
indeed is the strongly interacting medium and not the nuclear effects, the process-dependent power cor- 
rections and/or nPDFs, in the production of 10 ... 20 GeV partons. A significant suppression in p(d)A 
collisions would indicate that the nuclear effects, the process-dependent power corrections in particular, 
would play a major role in suppressing the hadron spectrum, especially so in ^4^4 collisions. In this 
case, the observed suppression in AuAu collisions would be only partly due to the dense medium. On 
the other hand, a significant enhancement (Cronin effect) in p(d)A collisions would also be a signal of 
sizable nuclear effects which would cause an even stronger enhancement of 10 ... 20 GeV parton pro- 
duction in AuAu collisions. In this case, the suppression observed in the AuAu data would, in fact, be 
stronger than expected based on the pp reference data. The new data on single hadron px distributions 
in dAu collisions at y/s = 200 GeV at RHIC [6-8] confirm the latter case by showing the absence of 



suppression or even the appearance of a Cronin-type enhancement up to p? ~ 8 GeV, in contrast to the 
strong suppression observed in AuAu collisions. The experience at RHIC clearly demonstrates that in 
order to measure the detailed properties of the strongly interacting matter (Equation of State, degrees 
of freedom, degree of thermalization, etc.) produced at the LHC, one needs precise knowledge of the 
absolute reference cross sections where the nuclear effects are constrained by those measured in pA 
collisions. 

Apart from the connection to the QGP signals, we emphasize the importance of measuring the 
magnitude of the nuclear effects on hard and semi-hard probes in pA relative to pp. It is only through 
these measurements, where we can separate the universal nuclear dependence of the nPDFs from the 
process-dependent power corrections due to rescattering in nuclear matter, that we can differentiate the 
contributions to the nuclear dependence. Then it is possible to determine when and where new QCD 
phenomena, such as gluon saturation at the scale of a few GeV, set in. 

Finally, in Section 7, "Novel QCD phenomena in pA collisions at the LHC", we emphasize the fact 
that pA collisions per se offer lots of new interesting QCD phenomena to explore, such as the breaking 
down of the leading twist QCD for a wide range of momentum transfer due to increase of the gluon fields 
at small x, related phenomena of pr broadening of the forward spectra of leading jets and hadrons, strong 
quenching of the low pr hadron production in the proton fragmentation region, and coherent diffraction 
into three jets. 

In summary, the pA program at the LHC serves a dual role. It is needed to calibrate the AA mea- 
surements for a sounder interpretation. It also has intrinsic merits in the framework of a more profound 
understanding of QCD (e.g. shadowing vs. diffraction, nonlinear QCD and saturation, higher twists...). 
Thus, one should foresee three key steps in this program: (1) (pp and) pPb runs with reliable deter- 
mination of centrality at the same collision energy as PbPb interactions; (2) a systematic study of pA 
collisions, requiring a variety of collision energies and nuclei as well as a centrality scan; and (3) an in- 
terchange of p- and ^4-beams for asymmetric detectors. The expected physics output from measurements 
of hard and semi-hard probes in pA collisions at the LHC may be summarized as follows. We are able to 

• test the predictive power of QCD perturbation theory in nuclear collisions by verifying the appli- 
cability of factorization theorems and the universality of the nPDFs through the hardest probes, 
only available at the LHC — ► Sees. 3 and 5. 

• measure the nuclear effects (internal to the nucleus) in the nPDFs over an unprecedented range of 
x and Q; investigate the interplay between the "EMC" effect, nuclear shadowing and saturation as 
well as the transverse-coordinate dependence of the nPDFs; and possibly, discover a new state of 
matter, the color glass condensate, by probing very soft gluons in heavy nuclei through the rapidity 
dependence — ► Sec. 4. 

• determine the nuclear dependence of the cross sections of the semi-hard probes in pA collisions 
and study QCD multiple partem scattering in nuclear matter and its corresponding dependences 
beyond the universal nuclear effects included in the nPDFs — ► Sec. 6. 

• extract excellent information on QCD dynamics in hadronization because normal nuclear matter 
acts like a filter for color neutralization and partem hadronization and explore potential new QCD 
phenomena in pA collisions, such as diffraction into three jets, double PDFs, etc. — ► Sec. 7. 

• use the hard and semi-hard probe cross sections as references for the QGP signals in A A collisions: 
the hard probes set the benchmark of the applicability of factorization while the semi-hard probes 
help to understand the size of the nuclear modifications not caused by the dense medium produced 
in A A collisions. 



2. THE EXPERIMENTAL PARAMETERS FOR pA AT THE LHC 
2.1 Constraints from the LHC machine 

Andreas Morsch 



2.11 pA collisions as part of the LHC ion programme 

The LHC Ion programme has been reviewed in a recent LHC Project Report [9]. Whereas collisions 
between light ions are already part of the LHC Phase II programme, pA and pA-Yike collisions (i.e. pA, 
dA, a A) are part of the LHC upgrade programme. Collisions between protons and ions are in princi- 
ple possible since two independent RF systems are already foreseen in the LHC nominal programme. 
However, pA collisions also require the availability of two independent timing systems. The necessary 
cabling is already included in the baseline layout. Acquisition of the dedicated electronics remains to be 
discussed. 

The question whether dA collisions will be available in the LHC is still under investigation be- 
cause the same source has to be used for the production of both protons and deuterons. Since proton 
availability is required for other CERN programs in parallel with the LHC pA runs, the final answer 
will depend on the time required for the source to switch back and forth between proton and deuteron 
operation. Production of a particles in parallel with proton operation is simpler and can be considered 
as an interesting alternative to pA and dA collisions. 

As an alternative, the PS division is presently studying a new layout for Linac 3, with possibly two 
sources and a switch-yard allowing simultaneous production of two different types of ions, including 
deuterons and a. 



2.12 Beam properties 

At the LHC, protons and ions have to travel in the same magnetic lattice and, hence, the two beams have 
the same rigidity. This has several consequences for the physical properties of the collision system. The 
condition of having the same rigidity results in beam momenta p which depend on the atomic mass A 
and charge Z of the beam particles. For the nominal LHC bending field one obtains 



p = 7Tev4 • (1) 
A 



The centre of mass energy is given by 



^= 14TeV Vtt ' <2) 

Since, in general, the two beams have different momenta, the rest system of the colliding particles 
does not coincide with the laboratory system. The rapidity of the centre of mass system relative to the 
laboratory is 

1 Z X A 2 

In Table 2 the properties of several collision systems of interest have been summarized together 
with their geometric cross sections. 



2.13 Luminosity considerations for pA collisions 

Two processes limit the proton and ion beam lifetimes and, hence, the maximum luminosity [10]: 

1 Beam-beam loss is due to hadronic interactions between the beam particles and proportional to the 
geometric cross section (see Table 2). For asymmetric beam intensities (different intensity proton 
and ion beams), the lifetime is proportional to the intensity of the higher intensity beam. 



Table 2: Geometric cross section, maximum centre of mass energy and rapidity shift for several collision systems. 
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2 Intra-beam scattering increases the beam emittance during the run. Since the timescale of emit- 
tance growth is proportional to AZ, it mainly affects the ion beam. 

It can be shown that maximizing the beam lifetime leads to solutions with asymmetric beam 
intensities (high-intensity proton beam, low-intensity ion beam). However, for the pA luminosity re- 
quirements published by the ALICE collaboration [11], it is rather the source intensity than the lumi- 
nosity lifetime which limits the luminosity. Since high intensity proton beams are readily available 
(jymax _ io 11 /bunch) but ion sources can deliver only relatively small intensities, mainly limited by 
space charge effect in the SPS, these conditions also push the realistic scenarios towards asymmetric 
beams. 

Table 3 shows the maximum intensities per bunch for Pb, Ar, and O. It has to be emphasized that 
these maximum intensities represent rough estimates and much more detailed studies are needed 
to confirm them. They are by no means to be quoted as official LHC baseline values. At these 
maximum intensities the transverse emittance growth time amounts to about 35 h. In the same table, the 
proton intensities needed to obtain the luminosities requested by ALICE are shown. In all cases these 
intensities are below iV™ ax . Hence, higher luminosities from which CMS and ATLAS could profit can 
be envisaged, as shown in the last column. 

Table 3: Maximum ion intensity per bunch, ALICE luminosity request, the proton intensity needed to fulfill the ALICE request 
and the maximum luminosity for N p = 10 11 for three different pA systems. 
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2.2 The ALICE detector 

Michiel Botje 

2.21 Introduction 

To study the properties of hadronic matter under conditions of extreme energy density the ALICE col- 
laboration has designed a general purpose detector optimised to measure a large variety of observables 
in heavy ion collisions at the LHC [12]. The apparatus will detect and identify hadrons, leptons and 
photons over a wide range of momenta. The requirement to study the various probes of interest in a 
very high multiplicity environment, which may be as large 8000 charged particles per unit of rapidity in 
central PbPb collisions, imposes severe demands on the tracking of charged particles. This requirement 
has led to a design based on high granularity, but slow, drift detectors placed in a large volume solenoidal 
magnetic field. Additional detectors augment the identification capabilities of the central tracking sys- 
tem which covers the pseudorapidity range \q\ < 0.9. A forward muon spectrometer (2.5 < r] < 4.0) 
provides measurements of the quarkonia state. Forward detectors provide global event properties by 
measurements of photon and charged particle multiplicities and forward calorimetry. 

In addition to running with Pb beams at the highest available energy of ^fs = 5.5 TeV, ALICE will 
take proton-proton (pp) and pA (pPb, maybe dPb or ctPb) reference data. Note that there are restrictions 
imposed on the pA programme by both the LHC machine and the ALICE detector. First, the maximum 
luminosity is limited because of pile-up in the TPC to typically (10 30 , 10 29 , 10 27 ) cm _2 s~ 1 for pp, pPb 
and PbPb collisions, respectively. Second, the energy per nucleon of ions in the LHC bending field 
scales to the nominal proton beam energy (7 TeV) with the ratio Z/A so that the centre of mass energy 
per nucleon and the rapidity shift of the centre of mass with respect to the LHC system are ^fs rs 
(9, 6, 6) TeV and Ay « (0.5, 0, 0) for pPb, dPb and aPb collisions, respectively, as given by Eqs. 
2 and 3. A comparison of nominal pPb data with pp (14 TeV) or PbPb (5.5 TeV) therefore needs a 
good understanding of the energy dependence. Furthermore, the rapidity shift of Ay rj 0.5 in pPb 
collisions needs a dedicated study of the ALICE acceptance. Whereas detailed simulations of PbPb and 
pp collisions are already well under way, a comprehensive study of the ALICE pA acceptance is not yet 
available. 

In the following we will briefly describe the main components of the ALICE detector. For a recent, 
more detailed review, see Ref. [13] and references therein. 

2.22 The ALICE detector 

The layout of the ALICE detector is shown in Fig. 1 . The central barrel is placed inside the L3 magnet 
which provides a solenoidal field of up to 0.5 T The choice of field is a compromise between low mo- 
mentum acceptance (~ 100 MeV/c), momentum resolution and tracking and trigger efficiency. ALICE 
intends to run mostly with a higher field of 0.5 T to study high pr probes but also with a lower field 
configuration of 0.2 T which is optimal for soft hadronic physics. 

The central tracking detector covers a range of \rj\ < 0.9 and the full range in azimuth. It consists 
of an inner tracking system (ITS) with six layers of high resolution silicon detectors located at a radial 
distance from the beam axis of r = 4, 7, 15, 24, 39 and 44 cm, respectively. The ITS provides tracking 
and identification of low px particles and improved momentum resolution for particles at higher mo- 
menta which traverse the TPC. An impact parameter resolution of better than 100 pm allows for vertex 
reconstruction of hyperons and heavy quarks. 

A time projection chamber (TPC) of 0.9 (2.5) m inner (outer) radius provides track finding, 
momentum measurement and particle identification via dE/dx with a resolution of better than 10%. 
Hadrons are identified in the TPC and ITS in the range ~ 100-550 MeV/c and up to 900 MeV/c for 
protons. The overall efficiency of the ALICE tracking is estimated to be better than 90% independent of 
p T down to 100 MeV/c. 




Fig. 1 : Layout of the ALICE detector. 



One of the unique features of ALICE is its particle identification capability. First of all, a transition 
radiation detector (TRD) with an inner (outer) radius of 3 (3.5) m provides electron identification. The 
TRD allows for the measurement of light and heavy meson resonances, the dilepton continuum and open 
charm and beauty via their semileptonic decay channels. The fast tracking capability of the TRD can be 
used to trigger on high px leptons and hadrons, providing a leading-particle trigger for jets. 

Time-of-flight counters (TOF), covering the TPC geometrical acceptance, are installed at a radius 
of 3.5 m and have a time resolution of better than 150 ps. This detector will be able to identify ir, K and 
protons in the semi-hard regime up to about 3 GeV/c. 

A single arm, high resolution electromagnetic calorimeter (PHOS) is located 5 m below the in- 
teraction point and covers a range of |r/| < 0.12 in pseudorapidity and A<p = 100° in azimuth. The 
calorimeter is optimised to measure ir°, rj and prompt photons in the energy range of one to several tens 
ofGeV. 

A RICH detector (HMPID) is positioned on top of the ALICE detector at about 4.5 m from the 
beam axis covering a range of \rj\ < 0.5 and A<p = 57°. The RICH will extend the particle identification 
capabilities of the central barrel to « 3 GeV/c for tt/K and « 5 GeV/c for K/p separation. 

The forward muon spectrometer consists of an absorber, starting at 1 m from the vertex, followed 
by a dipole magnet with 3 Tm field integral, placed outside the central magnet, and 10 planes of high 
granularity tracking stations. The last section of the muon arm consists of a second absorber and two 
more tracking planes for muon identification and triggering. This detector is designed to measure the 
decay of heavy quark resonances (J/ijj, tp', T, T' and T") with sufficient mass resolution to separate the 
different states. The coverage is 2.5 < i] < 4 with pr( single muon) > 1 GeV/c. 

Several smaller detector systems located at small angles will measure global event characteristics. 
In particular, the FMD detectors measure the charged particle multiplicity in the ranges —5.1 < rj < 
— 1.7 and 1.7 < rj < 3.4 and full azimuth. The FMD and ITS combined will then cover most of the 
phase space (—5.1 < r\ < 3.4) for multiplicity measurements. Initial Monte Carlo studies indicate that 



a centrality determination in pA collisions with a resolution of about 2 fm might be within reach from a 
measurement of the total charged multiplicity. 

To measure and to trigger on the impact parameter in heavy ion collisions, the neutron and proton 
spectators will be measured by zero degree calorimetry. The neutron and proton spectators are separated 
in space by the first LHC dipole and measured by a pair of neutron and proton calorimeters installed at 
about 90 m both left and right of the interaction region. A centrality determination in pA collisions can 
be provided by measuring black (p < 250 MeV, where p is the momentum in the rest system of nucleus 
A) and grey (p = 250-1000 MeV) protons and neutrons which are correlated with the impact parameter 
or, equivalently, with the number of collisions [14]. Initial simulations show that the ZDC acceptance for 
slow nucleons in pA collisions is quite large, see Fig. 2, so that such a centrality measurement in ALICE 
might be feasible. More detailed studies are currently under way. 




Fig. 2: Hit distribution of slow nucleons from pPb collisions at the front face of the ZDC. The squares at (—20) cm indicate 
the neutron (proton) calorimeter. The Pb beam is visible inside one of the beam pipes (shown as circles) at about x — — 9 cm. 



2.3 The CMS detector 

Walter Geist 

The CMS detector at the LHC was conceived for optimal measurements of hard processes in pp 
collisions at a center-of-mass energy yfs = 14 TeV and at a very high luminosity. Among the rare 
and hard processes are e.g. the production of Higgs- and SUSY-particles corresponding to mass/energy 
scales beyond 100 GeV; their decays yield charged leptons, photons, and weakly interacting particles as 
well as jets of hadrons. The CMS design is therefore driven by the need for 

( i) a superb muon-system, 

(ii) the best possible electromagnetic calorimeter (ECAL) compatible with (i), and 
( Hi) a high-performance tracking detector. The set-up is completed by 

(iv) a large acceptance hadronic calorimeter (HCAL). 




Compact Muon Solenoid 

Fig. 3: Layout of the CMS detector. 



Magnet. A direct consequence of these scientific goals is a strong (4 T) compact solenoid centered at 
the beam crossing; it bends the trajectories of charged particles in the plane perpendicular to the beam 
direction and thus facilitates the use of the position of the small diameter beam for establishing primary 
and secondary vertices. Its length makes it possible to reconstruct tracks with good momentum resolution 
up to pseudorapidities \rj\ < 2.4. The field integral leads to transverse momentum resolutions 5pt/pt 
smaller than about 15% atpx = 1 TeV/c. At the same time, particles with transverse momenta less than 
about 1 GeV / c spiral close to the beam line thereby reducing detector occupancies farther away. 



Tracker. The all-Si tracker system of CMS consists of 13 cylindrical layers (3 layers of pixel detectors 
and 10 layers of strip-modules) as well as of forward-backward discs (2 pixel discs and 12 Si-strip discs) 
perpendicular to the beams. They complete the geometric acceptance for tracks up to \rf\ = 2.4. The 
Si-detectors are characterized by fast response and very good spatial resolution. Four Si-strip layers 
are made of back-to-back ('double sided') modules for quasi-3-dimensional space points. Good spatial 
resolution and low occupancies (< 1%) in pp collisions at high luminosity result in track reconstruction 
efficiencies usually better than 98%, and better than 95% inside jets. The relative resolution for transverse 
momenta is expected to be Spr/pr = (15pt + 5)%, with px in units of TeV, for \rj\ < 1.6. Furthermore, 
this system, in particular the pixel part, provides a reliable separation of primary vertices and secondary 
vertices due to S-decays. 

Calorimeters. The electromagnetic and hadronic calorimeters surround the tracker system hermeti- 
cally. They are positioned inside the coil of the magnet. Their large rapidity coverage ensures efficient 
detection of 'missing' transverse energy in the case of weakly-interacting particle production. 

• ECAL: The barrel section of ECAL covers the rapidity range up to \rj\ = 1.48. The forward 
part extends the acceptance to \rj\ = 5, where preshower detectors are also foreseen. The fine 
segmentation of the PbW04 crystals in rapidity and azimuth matches the transverse shower size 
from single and photons. 

• HCAL: the hadronic calorimeter is also divided into a barrel part using copper-plastic tiles for 
| r/ 1 < 3 and a forward calorimeter which increases the acceptance up to 1 1]\ = 5 and is based upon 
a fibre concept. 

Muon detection. The muon chambers outside the magnet coil are grouped into stations of gaseous 
detector layers interleaved with the iron of the return yoke. They are protected by an absorber thickness 
of 10 interaction lengths in front of the first station and an additional 10 interaction lengths before the last 
station. Also, the muon system is organized into cylindrical and forward disc parts. It is a stand-alone 
system, capable of determining the muon momenta to very good precision. 

Trigger. The CMS triggers select rare processes from the overwhelming rate of soft reactions at very 
high luminosities. The total rate of events to be transferred to a storage medium is limited to about 
100 Hz. In order to cope with this low rate, the rate of individual triggers must be limited, typically by 
imposing a transverse momentum/energy cutoff. 

2.31 CMS and pA collisions 

While the first detailed investigations of CMS performance with PbPb collisions have been undertaken, 
this is not the case for pA collisions. However, no major problems are anticipated [15]. There are no 
a priori technical obstacles to running the CMS experiment at the lower pA- and ^^-luminosities (see 
Sec. 2.1) with e.g. lower ^/energy thresholds and/or a lower magnetic field, as long as occupancy 
problems do not obstruct track reconstruction substantially. Global measurements corresponding to an 
increased acceptance at lower transverse momenta are very useful, given that the p^-collision energy 
will exceed that previously available by more than order of magnitude. 

It should be emphasized that the TOTEM experiment will take data in conjunction with CMS. 
TOTEM needs a geometrical acceptance for secondary particles up to | rj \ = 7. Furthermore, a 'zero- 
degree' calorimeter covering the range of beam rapidities is considered for CMS heavy ion runs. 



2.4 The ATLAS detector 

Brian Cole, Helio Takai, and Stefan Tapprogge 

2. 41 The ATLAS detector 

The ATLAS detector is designed to study proton-proton collisions at the LHC design energy of 14 TeV 
in the center of mass. The extensive physics program pursued by the collaboration includes the Higgs 
boson search, searches for SUSY and other scenarios beyond the Standard Model, as well as precision 
measurements of processes within (and possibly beyond) the Standard Model. To achieve these goals 
at a full machine luminosity of 10 34 cm _2 s _1 , ATLAS will have a precise tracking system (Inner De- 
tector) for charged particle measurements, an as hermetic as possible calorimeter system, which has an 
extremely fine grain segmentation, and a stand-alone muon system. An overview of the detector is shown 
in Fig. 4. 




Fig. 4: The overall layout of the ATLAS detector 



The Inner Detector is composed of (1) a finely segmented silicon pixel detector, (2) silicon strip 
detectors (Semiconductor Tracker (SCT)) and (3) the Transition Radiation Tracker (TRT). The segmenta- 
tion is optimized for proton-proton collisions at design luminosity. The inner detector system is designed 
to cover a pseudorapidity of \q\ < 2.5 and is located inside a 2 T solenoid magnet. 

The calorimeter system in the ATLAS detector surrounding the solenoid magnet is divided into 
electromagnetic and hadronic sections and covers pseudorapidity |t/| < 4.9. The EM calorimeter is an 
accordion liquid argon device and is finely segmented longitudinally and transversely for \rj\ < 3.1. 
The first longitudinal segmentation has a granularity of 0.003 x 0.1 (At/ x A0) in the barrel and slightly 
coarser in the endcaps. The second longitudinal segmentation is composed of At/ x A(f> = 0.025 x 0.025 



Table 4: ATLAS detector acceptance for pA runs 
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cells and the last segment At] x A<fi = 0.05 x 0.05 cells. In addition, a finely segmented (0.025 x 0.1) 
pre-sampler system is present in front of the electromagnetic (EM) calorimeter. The overall energy 
resolution of the EM calorimeter determined experimentally is 10%/y/E © 0.5%. The calorimeter also 
has good pointing resolution, 60 mrad / \f~E for photons, and a timing resolution of better than 200 ps for 
energy showers larger than 20 GeV. 

The hadronic calorimeter is also segmented longitudinally and transversely. Except for the end- 
caps and the forward calorimeters, the technology utilized for the calorimeter is a lead-scintillator tile 
structure with a granularity of Ar) x Atp = 0.1 x 0.1. In the endcaps the hadronic calorimeter is imple- 
mented in liquid argon technology for radiation hardness with the same granularity as the barrel hadronic 
calorimeter. The energy resolution for the hadronic calorimeters is 50%/s/E © 2% for pions. The very 
forward region, up to r] = 4.9, is covered by the Forward Calorimeter implemented as an axial drift 
liquid argon calorimeter. The overall performance of the calorimeter system is described in Ref. [16]. 

The muon spectrometer in ATLAS is located behind the calorimeters, thus shielded from hadronic 
showers. The spectrometer is implemented using several technologies for tracking devices and a toroidal 
magnet system which provides a 4 T field to have an independent momentum measurement outside the 
calorimeter volume. Most of the volume is covered by Monitored Drift tubes (MDTs). In the forward 
region, where the rate is high, Cathode Strip Chamber technology is chosen. The stand-alone muon 
spectrometer momentum resolution is of the order of 2% for muons with p? in the range 10 . . . 100 GeV. 
The muon spectrometer coverage is \r]\ < 2.7. 

The trigger and data acquisition system of ATLAS is a multi-level system which has to reduce the 
beam crossing rate of 40 MHz to an output rate to mass storage of 0(100) Hz. The first stage (LVL1) is 
a hardware based trigger, making use of coarse granularity calorimeter data and dedicated muon trigger 
chambers only, to reduce the output rate to about 75 kHz, within a maximum latency of 2.5 fis. 

The performance results mentioned here were obtained using a detailed full simulation of the 
ATLAS detector response with GEANT and have been validated by an extensive program of test beam 
measurements of all components. 

2. 42 Expected detector performance for pA 

The proton-nucleus events are asymmetric, with the central rapidity shifted by almost 0.5 units. In 
proton-nucleus runs it is expected that all of the subsystems will be available for data analysis. For the 
expected luminosity of 10 29 cm _2 s _1 in pPb collisions, the occupancy of the detectors will be less than 
or at most comparable to the high luminosity proton-proton runs. The tracking will also benefit from not 
having up to 25 collisions (as in pp interactions at the design luminosity of 10 34 cm~ 2 s _1 ) in the same 
bunch crossing, but only a single event vertex. The detector acceptance for a pA run is the same as for 
pp and is shown in Table 4 (in the detector frame. See Sec. 2.1 for the rapidity shifts in the cms frame of 
the pA collisions). 

The muon spectrometer will also benefit from the low luminosity proton-nucleus mode. In the high 



luminosity proton-proton runs, a large number of hits in the spectrometer comes from slow neutrons from 
previous interactions. These are not present in low luminosity runs. The only expected backgrounds 
are muons from ir and K decays in flight. As has been shown for pp, a substantial number of these 
background muons are expected to be rejected when matching the track from the stand-alone muon 
system to a track in the Inner Detector is required. 

2.43 The physics potential of ATLAS in pA 

Measurements of proton-nucleus collisions at the LHC will not only provide essential control of the hard 
and semi-hard processes in heavy ion collisions, but these measurements can also address physics that 
is, in its own right, of fundamental interest. The ATLAS detector provides an unprecedented opportunity 
to study proton-nucleus collisions in a detector with both large acceptance and nearly complete coverage 
of the various final states that can result from perturbative QCD processes. 

An important physics topic that can be addressed by proton-nucleus measurements in ATLAS is 
the applicability of factorization in hard processes involving nuclei. No detector that has studied either 
nuclear deep inelastic scattering or proton-nucleus collisions has been able to simultaneously study many 
different hard processes to explicitly demonstrate that factorization applies. In ATLAS, we will have the 
opportunity to study single jet and jet — jet events, photon — jet events, heavy quark production, Z° 
and W ± production, Drell-Yan dilepton production, etc. and probe the physics of jet fragmentation in 
detail. With the wide variety of available hard processes that all should be calculable with the same 
partem distributions, we should be able to clearly demonstrate the success of factorization for the large- 
Q 2 processes. 

The ATLAS inner detector has been designed to have good tracking efficiency in the momentum 
range of 1 GeV up to 1 TeV over the pseudorapidity range \rj\ < 2.5. The ability to reconstruct jets in 
ATLAS over a large pseudorapidity region (up to |r/| < 4.9) and over the full azimuth would probe the 
onset of saturation effects in high pt jets. One aspect of saturation that is not yet well understood is the 
effect on hard scattering processes for Q 2 ~ Q 2 , where Q s is the scale at which saturation sets in. It has 
been suggested that the gluon distribution could be significantly modified well above Q s . If so, then 
ATLAS should be able to study the Q 2 evolution of the saturation effects using photon — jet, bb and 
jet — jet measurements. 

Detailed studies of the behavior of jet fragmentation as a function of pseudorapidity could be 
used to determine whether the jet fragmentation is modified by the presence of the nucleus or its large 
number of low-x gluons. The observation or lack thereof of modifications to jet fragmentation could 
provide sensitive tests of our understanding of formation time and coherence in the re-dressing of the 
hard scattered parton and the longitudinal spatial spread of low-x gluons in a highly Lorentz-contracted 
nucleus. 



3. QCD FACTORIZATION AND RESCATTERING IN pA COLLISIONS 

Jianwei Qiu 

Perturbative QCD has been very successful in interpreting and predicting hadronic scattering pro- 
cesses in high energy collisions, even though the physics associated with an individual hadron wave 
function is nonperturbative. It is the QCD factorization theorem [17] that provides prescriptions to sep- 
arate long- and short-distance effects in hadronic cross sections. The leading power contributions to a 
general hadronic cross section involve only one hard collision between two partons from two incoming 
hadrons of momenta p A and pg. The momentum scale of the hard collision is set by producing either a 
heavy particle (like W/Z or virtual photon in Drell-Yan production) or an energetic third partem, which 
fragments into a hadron h of momentum p' . The cross section can be factorized as [17] 

Eh d*AB-*W = J2j dx'f j/B (x') j dxf i/A (x) j dzD h/k (z)E h ^^(x PA ,x'p B ,^), (4) 
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where Y^ijk runs over an partem flavors and all scale dependence is implicit. The f^ A are twist-2 dis- 
tributions of parton type i in hadron A, and the D h / k are fragmentation functions for a parton of type 
k to produce a hadron h. For jet production, fragmentation from a parton to a jet, suitably defined, is 
calculable in perturbation theory and may be absorbed into the partonic hard part a. For heavy particle 
production, the fragmentation function is replaced by (5(1 — z). 

The factorized formula in Eq. (4) illustrates the general leading power collinear factorization the- 
orem [17]. It consistently separates perturbatively calculable short-distance physics into a, and isolates 
long-distance effects into universal nonperturbative matrix elements (or distributions), such as f i i A or 
D h j k , associated with each observed hadron. Quantum interference between long- and short-distance 
physics is power-suppressed by the large energy exchange of the collisions. Predictions of pQCD follow 
when processes with different hard scatterings but the same set of parton distributions and/or fragmenta- 
tion functions are compared [18]. 

With the vast data available, the parton distributions of a free nucleon are well determined by QCD 
global analyses [19,20]. Thanks to recent efforts, sets of fragmentation functions for light hadrons are 
becoming available though the precision is not as good as the parton distributions due to the limited data 
[21,22]. 

Studies of hard processes at the LHC will cover a very large range of longitudinal momentum 
fractions x of parton distributions: x > XTe y /(2 — xxe~ y ) with xt = 2pr/\/s for inclusive jet 
production in Eq. (4), where y and pt are the rapidity and transverse momentum of the produced jets, 
respectively. For the most forward or backward jets, or low pt Drell-Yan dileptons, x can be as small 
as 10~ 6 at yfs = 14 TeV. The number of gluons having such small x and transverse size Ar± oc 1/pr 
may be so large that gluons appear more like a collective wave than individual particles, and a new 
nonperturbative regime of QCD, such as gluon saturation or the color glass condensate [23], might be 
reached. 

The use of heavy ion beams allows one to enhance the coherent effects because of the larger 
number of gluons. In the small x regime of hadron-nucleus collisions, a hard collision of a parton of 
the projectile nucleon with a parton of the nucleus occurs coherently with all the nucleons at a given 
impact parameter. The coherence length (~ 0.1/x fm) far exceeds the nuclear size. To distinguish 
parton-nucleus multiple scattering from partonic dynamics internal to the nucleus, we classify multiple 
scattering in the following three categories: (a) interactions internal to the nucleus, (b) initial-state parton- 
nucleus interactions (ISI), and (c) final-state parton-nucleus interactions (FSI), as shown in Fig. 5 [18]. 

Interactions internal to the nucleus change the parton distributions of the nucleus, as shown in 
Fig. 5a. As a result, the effective parton distributions of a large nucleus are different from a simple sum 
of individual nucleon parton distributions. Since only a single parton from the nucleus participates in the 
hard collision to leading power, the effect of the initial-state interactions internal to the nucleus does not 




Fig. 5: Classification of parton multiple scattering in nuclear Fig. 7: Soft interactions of long-range fields that might en- 
medium: (a) interactions internal to the nucleus, (b) initial- hance the gluon density in nuclear medium and affect the 
state interactions, and (c) final-state interactions. FSI. 



change the hard collisions between two incoming partons. This preserves the factorized single scattering 
formula in Eq. (4) [24] except now the twist-2 parton distributions j i j^ L are replaced by the corresponding 
effective nuclear parton distributions, which are defined in terms of the same operators but on a nuclear 
state. Such effective nuclear parton distributions include the "EMC" effect and other nuclear effects so 
that they differ from the parton distributions of a free nucleon. However, they are still twist-2 distribution 
functions by definition of the operators of the matrix elements and are still universal. 

Because of the twist-2 nature of the nuclear parton distributions, their nuclear dependence can 
only come from the A-dependence of the nonperturbative input parton distributions at a low momentum 
scale Qo and the ^4-dependent modifications to the DGLAP evolution equations from C}o to Q due to 
the interactions between partons from different nucleons in the nucleus [25,26]. The nonperturbative 
nuclear dependence in the input distributions can be either parameterized with parameters fixed by fits to 
experimental data [27-29] or calculated with some theoretical inputs [30, 31]. 

Let hard probes be those whose cross sections are dominated by the leading power contributions 
in Eq. (4). That is, nuclear dependence of the hard probes comes entirely from that of nuclear parton 
distributions and is universal. Because of the wide range of x and Q covered, the hard probes in hadron- 
nucleus collisions at the LHC directly detect the partonic dynamics internal to the nucleus and provide 
excellent information on nuclear parton distributions. The knowledge of these distributions is very useful 
for understanding gluon saturation, a new nonperturbative regime of QCD [31, 32]. 

On the other hand, the hadronic cross sections receive power-suppressed corrections to Eq. (4) 
[33-35]. These corrections can come from several different sources, including the effect of partonic 
non-collinear momentum components and effects of non-vanishing invariant mass of the fragmenting 
parton k, as well as the effects of interactions involving more than one parton from each hadron, as 
shown in Fig. 5b and 5c. Although such multiple coherent scattering is formally a higher-twist effect and 
suppressed by powers of the large momentum scale of the hard collision, the corrections to the leading 
power factorized formula in Eq. (4) are proportional to the density of additional scattering centers and 
can be substantial due to a large density of soft gluons in the nucleus-nucleus collisions. For the ISI, the 



number of gluons available at the same impact parameter of the hard collision is enhanced due to the 
large nucleus, while the densities of soft partons available to the FSI can be even larger since they can 
come either from the initial wave functions of the colliding nuclei or be produced in the long-range soft 
parton interactions between different nucleons at the time when the hard collision took place. 

Consider the scattering of an elementary particle or a parton (a quark or a physically polarized 
gluon) of momentum xp in nuclear matter, as shown in Fig. 6. A hard scattering with momentum transfer 
Q can resolve states whose lifetimes are as short as 1/Q [18]. The off-shellness of the scattered parton 
increases with the momentum transfer simply because the number of available states increases with 
increasing momentum. Typically, the scattered parton (say, of momentum p{) is off-shell by mj < Q, 
where mj is the invariant mass of the jet into which parton fragments. Further interactions of the off- 
shell parton are suppressed by an overall factor of 1/mj ~ 1/Q 2 , since the effective size of the scattered 
parton decreases with momentum transfer, and by the strong coupling evaluated at scale mj ~ Q. Thus 
rescattering in nuclear collisions is suppressed by a s (Q)/Q 2 compared to single scattering. 

Counting of available states ensures that mj > ]C/i(-^fr) m /i ^ Aqcd> where J2h runs over 
all hadron types in the jet and (N^) is the corresponding multiplicity. On the other hand, if we are to 
recognize the jet, we must have mj <C Ej = p\, where Ej is the jet energy. In the rest frame of 
the nucleus, the scattered parton has a lifetime At ~ (Ej/mj)/mj. Thus, at high enough jet energy, 
At > Ra, the lifetime of the scattered parton will exceed the size of nuclear matter even though the 
parton itself is far off mass shell. Then the interactions of the scattered off-shell parton with nuclear 
matter may be treated by the formalism of pQCD [18]. 

In order to consistently treat the power-suppressed multiple scattering, we need a factorization 
theorem for higher-twist (i.e., power suppressed) contributions to hadronic hard scattering. It was shown 
in Ref. [33] that the first power-suppressed contribution to the hadronic cross section can be factorized 
into the form 

E h A y } = £ J dx'f j/B (x')j dzD h/k (z) (5) 
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where the partonic hard part a^l^j, is infrared safe and depends on the identities and momentum 
fractions of the incoming partons but is otherwise independent of the structure - in particular the size 
- of the hadron and/or heavy ion beams. The superscript "(4)" on the partonic hard part indicates the 
dependence on twist-4 operators [33,36]. The correlation functions Tt U iy A are defined in terms of 
matrix elements of twist-4 operators made of two pairs of parton fields of flavor i and i', respectively. 
For example, a correlation function between quark and gluon in a hadron h of momentum p can be 
expressed as 

T {qG y A ( Xl ,x 2 ,x 3 ) = j ' l! ' l ^! !h ^ : ' — 

x (p\TP(0)^F +a (y^F a + (y^4,(y^)\p) , (6) 

with quark field ip and gluon field strength F +a . 

Demonstrating factorization at the next-to-leading power is a starting point for a unified discussion 
of the power-suppressed effects in a wide class of processes. A systematic treatment of double scattering 
in a nuclear medium is an immediate application of the generalized factorization theorem [18]. Because 
of the infrared safe nature of the partonic hard part a , ^},^_ tk in Eq. (5), the nuclear dependence of the 
double scattering comes entirely from the correlation functions of two pairs of parton fields which can 



be linearly proportional to A 1 / 3 (or nuclear size) [36]. Therefore, if the scattered off-shell parton has a 
lifetime longer than the nuclear size, rescattering is enhanced by A 1 / 3 from the medium size and gets 
an overall suppression factor, a s (Q) A 1 / 3 X 2 /Q 2 , where A has the dimension of mass and represents 
the nonperturbative scale of the twist-4 correlation functions in Eq. (6). A semiclassical estimate gives 
A 2 ~ (F +a F a + )im 2 /TT [36]. 

For the ISI encountered by the incoming parton, A 2 is proportional to the average squared trans- 
verse field strength, (F +a F a + ), inside the nucleus, which should be more or less universal. From the 
data on Drell-Yan transverse momentum broadening, it was found [37] that A 2 SI » A^y ~ 0.01 GeV 2 . 
However, the numerical value of A 2 for the FSI does not have to be equal to the A 2 SI due to extra soft 
gluons produced by the instantaneous soft interactions of long-range fields between the beams at the 
same time when the jets were produced by two hard partons, as shown in Fig. 7. 

It was explicitly shown [38-40] that the corrections to hadronic Drell-Yan cross sections cannot be 
factorized beyond the next-to-leading power. However, it can be shown using the technique developed in 
Ref. [33] that the type of ^4 1 / 3 -enhanced power corrections to hadronic cross sections in hadron-nucleus 
collisions can be factorized to all powers, 

Eh ^V MP ° = £ fj/BiX ' ] ® Dh/k{z) ® T £)M (Xl) ® Eh ^a^^ X ' PB > 7>' (7) 
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where <g) represents convolutions in fractional momenta carried by the partons and T^ 2n \xi) represent 
the correlation functions of n pairs of parton fields with flavors i n = 1, 2, n. In Eq. (7), no power 
corrections are initiated by the higher-twist matrix elements of the incoming hadron B because such 
contributions do not have the A 1 / 3 -type enhancement in hadron-nucleus collisions. For example, the 
recently proposed reaction operator approach to multiple scattering [41,42] should fit into the type of 
factorization in Eq. (7) with the independent-scattering-center approximation to the T^ 2n \xi) and keep- 
ing the lowest order in a s to the partonic a^ 2n \ which can be obtained using a recursive method. To 
the contrary, in nucleus-nucleus collisions, even the nuclear size enhanced power corrections cannot be 
formally factorized beyond the next-to-leading power. 

Let semi-hard probes in hadronic collisions be those with a large momentum exchange as well 
as large power corrections. We expect that pQCD has good predictive power for semi-hard observables 
in hadron-nucleus collisions but the pQCD factorization approach does not work well for semi-hard 
observables in nucleus-nucleus collisions. 

In conclusion, the factorized single scattering formula in Eq. (4) remains valid for hard probes in 
nuclear collisions, except that the parton distributions are replaced by corresponding effective nuclear 
parton distributions which are independent of the hard scattering and thus universal. Hard probes in 
nuclear collisions at the LHC can provide excellent information on nuclear parton distributions and 
detect the partonic dynamics internal to the nucleus. 

However, the power-suppressed corrections to the single scattering formula can be substantial and 
come from several different sources. The effect of the off-shellness of the fragmenting parton k leads to a 
correction of the order (mj /pr) 2 with mj = Ylhi^f^^h- Both the ISI and FSI double scatterings give 
power-suppressed corrections proportional to a s A l / 3 \ 2 /p 2 -,, where A 2 SI is relatively small and almost 
universal and Ap SI is sensitive to the number of soft partons produced by the instantaneous collisions of 
long-range fields between nucleons in the incoming beams. 

Beyond double scattering (or next-to-leading power corrections), pQCD calculations might not be 
reliable due to the lack of factorization theorems at this level. However, pQCD factorization for the type 
of vl 1 / 3 -enhanced power corrections in hadron-nucleus collisions is likely to be valid to all powers. 



4. NUCLEAR PARTON DISTRIBUTION FUNCTIONS (nPDFs) 
4.1 Global DGLAP fit analyses of the nPDFs: EKS98 and HKM 

Kari J. Eskola, Heli Honkanen, Vesa J. Kolhinen, and Carlos A. Salgado 



4.11 Introduction 

Inclusive cross sections for hard processes A+B — > c+X involving a sufficiently large scale Q 3> Aqcd 
are computable by using collinear factorization. In the leading-twist approximation power corrections 
~ l/Q 2 are neglected and 



da(Q 2 , \/~s)ab^c+x = 



Z A f? /A ( Xl ,Q 2 ) + (A-Z A )f? /A ( Xl ,Q 2 ) 



Z B ff B (x2,Q 2 ) + {B- Z B )ff B (x 2 ,Q 2 ) 



dd-(Q 2 ,xi,x 2 )ij^ c +x (8) 



where A, B are the colliding hadrons or nuclei containing Z A and Zb protons respectively, c is the 
produced parton, x and X are anything, da(Q 2 ,xi,X2)ij^ c +x is the perturbatively calculable differential 
cross section for the production of c at the scale Q, and xi^ ~ Q/\fs are the fractional momenta of 
the colliding partons i and j. The number distribution function of the parton flavour i of the protons 
(neutrons) in A is denoted as ff^ A (f^ A ), and similarly for partons j in B. 

In the leading-twist approximation, multiple scattering of the bound nucleons does occur but all 

collisions are independent, correlations between partons from the same object A are neglected, and only 

p/A 

one-parton densities are needed. The parton distribution functions (PDFs) f\ are universal quanti- 
ties applicable in all collinearly factorizable processes. The PDFs cannot be computed by perturbative 
methods and, so far, it has not been possible to compute them from first principles, either. Thus, nonper- 
turbative input from data on various hard processes is needed for the extraction of the PDFs. However, 
once the PDFs are known at some initial (lowest) scale Qq 3> Aqcd, QCD perturbation theory predicts 
the scale evolution of the PDFs to other (higher) values of Q 2 in form of the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equations [43]. 

The method of extracting the PDFs from experimental data is well established in the case of the 
free proton: the initial nonperturbative distributions are parametrized at some and evolved to higher 
scales according to the DGLAP equations. Comparison with the data is made in various regions of the 
(x, Q 2 ) -plane. The parameters of the initial distributions, ff{x, Qq), are fixed when the best global fit 
is found. The data from deeply inelastic Ip scattering (DIS) are most important for these global DGLAP 
fits, especially the HERA data at small values of x and Q 2 . The sum rules for momentum, charge and 
baryon number give further constraints. In this way, through the global DGLAP fits, groups like MRST 
[44], CTEQ [19] or GRV [45] obtain their well-known sets of the free proton PDFs. 

The nuclear parton distribution functions (nPDFs) differ in magnitude from the PDFs of the free 
proton. In the measurements of the structure function, F 2 A = Z A F^ A + (A — Z A )F^ A , of nuclear 
targets in DIS (see e.g. Ref. [46] and references therein), and especially of the ratio 

A 2 \da lA /dQ 2 dx \F^(x,Q 2 ) 
R F2 (x,Q)= i d(jlD/dQ2dx « i F D {XjQ 2y W 

the following nuclear effects have been discovered in different Bjorken-x regions: 

• shadowing; a depletion, Rp 2 < 1, at x < 0.1, 

• antishadowing; an excess, Rp 2 > 1, at 0.1 < x < 0.3, 

• EMC effect; a depletion at 0.3 < x < 0.7, and 

• Fermi motion; an excess towards x — > 1 and beyond. 



The Q 2 dependence of Rp 2 is weaker and has thus been more difficult to measure. Data with high 
enough precision, however, exist. The NMC collaboration discovered a clear Q 2 dependence in the 
ratio da^ Sn / 'da^ c [47], i.e. the scale dependence of the ratio F 2 Sn /F 2 c , at x > 0.01. Since Fl 



■p(n)/A 



2~2q e q x lfg + fq^ n ^ A ] + 0(a s ), the nuclear effects in the ratio Rp 2 directly translate into nuclear 
effects in the parton distributions: f^ A ^ ff. 

The nPDFs, ff^ A , also obey the DGLAP equations in the large-Q 2 limit. They can be determined 
by using a global DGLAP fit procedure similar to the case of the free proton PDFs. Pioneering studies 
of the DGLAP evolution of the nPDFs are found in e.g. Ref. [48-51]. References for various other 
studies of perturbative evolution of the nPDFs and also to simpler (^-independent parametrizations of 
the nuclear effects in the PDFs can be found e.g. in Refs. [52, 53]. The nuclear case is, however, more 
complicated because of additional variables, the mass number A and the charge Z, and, because the 
number of data points available in the perturbative region is more limited than for the PDFs of the free 
proton. The DIS data play the dominant role in the nuclear case as well. However, as illustrated by Fig. 8, 
no data are available from nuclear DIS experiments below x<5-10~ 3 atQ 2 >l GeV 2 . This makes 
the determination of the nuclear gluon distributions especially difficult. Further constraints on the global 
DGLAP fits of the nPDFs can be obtained from e.g. the Drell-Yan (DY) process measured in fixed-target 
pA collisions [54,55]. Currently, there are two sets of nPDFs available which are based on the global 
DGLAP fits to the data: (i) EKS98 [27, 28] (the code in Refs. [56, 57]), and (ii) HKM [29] (the code in 
Ref. [58]). We shall compare the main features of these two analyses and comment on their differences 
below. 
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Fig. 8: The average values of x and Q of the DIS data from NMC [59-61] (triangles) and E665 [62, 63] (diamonds) in I A, and 
of X2 and M 2 of the DY dilepton data [54] (squares) in pA. The heavy quark mass scales are shown by the horizontal dashed 
lines. The initial scale Qq is mj? in EKRS and 1 GeV 2 in HKM. For the rest of the figure, see the text in Sec. 4.13. 



4. 12 Comparison ofEKS98 and HKM 

EKS98 and overview of constraints available from data. The parametrization EKS98 [27] is based 
on the results of the DGLAP analysis [28] and its follow-up [27]. We refer to these together as EKRS 



here. In the EKRS approach, the nPDFs, the parton distributions of the bound protons, f j , are defined 
through the modifications of the corresponding distributions in the free proton, 



where the PDFs of the free proton are from a known set such as CTEQ, MRS or GRV. As in the case 
of the free nucleons, for isoscalar nuclei the parton distributions of bound neutrons are obtained through 
isospin symmetry, /"^ = f^it and f™!t = . Although exact only for isoscalar and mirror nuclei, 
this is expected to be a good first approximation for all nuclei. 

To simplify the determination of the input nuclear effects for valence and sea quarks, the following 
flavour-independent initial modifications are assumed: R£ v (x,Qq) ~ R% (x,Qq) ps Ry(x, Qq), and 
R^(x,Qq) ps i?j(x,Qo) ~ R£( x :Qo) ~ Rs(x, Qo)- Thus only three independent initial ratios, Ry, 
Rg and Rq are to be determined at Qq = m 2 = 2.25 GeV 2 . Note also that the approximations above are 
needed and used only at Qq in [28]. In the EKS98-parametrization [27] of the DGLAP-evolved results, 
it was observed that, to a good approximation, R^ y = R§ and R£ = Rg for all Q 2 . Further details 
of the EKRS analysis can be found in Ref. [28], here we summarize the constraints available in each 
x-region. Consider first quarks and antiquarks: 

• At x > 0.3 the DIS data constrains only the ratio R y : valence quarks dominate F^, so that 
i?^ 2 ps R v but Rg and Rq are left practically unconstrained. An EMC effect is, however, assumed 
also for the initial Rg and Rq since in the scale evolution the EMC effect of valence quarks is 
transmitted to gluons and, further on, to sea quarks [50]. In this way the nuclear modifications Rf 
also remain stable against the evolution. 

• At 0.04 < x < 0.3 both the DIS and DY data constrain Rg and R v but in different regions of 
Q 2 , as shown in Fig. 8. In addition, conservation of baryon number restricts Ry. The use of DY 
data [54] is essential in order to fix the relative magnitude of Ry and Rg since the DIS data alone 
cannot distinguish between them. As a result, no antishadowing appears in Rg(x, Qq). 

• At 5 • 10~ 3 < x < 0.04, only DIS data exist in the region Q > 1 GeV where the DGLAP analysis 
can be expected to apply. Once R v is fixed by the DIS and DY data at larger x, the magnitude of 
nuclear valence quark shadowing in the EKRS approach follows from the conservation of baryon 
number. As a result of these contraints, nuclear valence quarks in EKRS are less shadowed than 
the sea quarks, R v > Rg. 

• At x < 5 • 10~ 3 , as shown by Fig. 8, the DIS data for the ratio Rp 2 lie in the region Q < 1 GeV 
where the DGLAP equations are unlikely to be applicable in the nuclear case. Indirect constraints, 
however, can be obtained by noting the following: (i) At the smallest values of x (where Q ^ 
1 GeV), Rj, 2 depends only very weakly, if at all, on x [60,62]. (ii) d(F 2 Sn / 'F?)/dlogQ 2 > 
at x ~ 0.01 [47], indicating that 0R f a / d log Q 2 > 0. (hi) Negative log Q 2 -slopes of R^ 2 at 
x < 5 • 10 -3 have not been observed [60]. Based on these experimental facts, the EKRS approach 
assumes that the sign of the log Q 2 -slope of Rp 2 remains non-negative at x < 0.01 and therefore 
the DIS data [60, 62] in the nonperturbative region gives a lower bound on Rp 2 at Qq for very 
small x. The sea quarks dominate over the valence quarks at small values of x so that only Rg is 
constrained by the DIS data. R v remains restricted by baryon number conservation. 

Pinning down the nuclear gluon distributions is difficult in the absence of stringent direct constraints. 
What is available in different regions of x can be summarized as follows: 

• At x > 0.2, no experimental constraints are currently available on the gluons. Conservation of 
momentum is used as an overall constraint but this alone is not sufficient to determine whether 
there is an EMC effect for gluons. Only about 30 % of the gluon momentum comes from x > 0.2 
so that fairly sizable changes in R G (x, Qq) in this region can be compensated by smaller changes 



at x < 0.2 without violating the constraints discussed below. As mentioned above, in the EKRS 
approach an EMC effect is initially assumed for Rq- This guarantees a stable scale evolution and 
the EMC effect remains at larger Q 2 . 

• At 0.02 < x < 0.2 the Q 2 dependence of the ratio F 2 Sn /F 2 c measured by NMC [47] currently 
sets the most important constraint on R G [64]. In the small-x region where gluons dominate 
the DGLAP evolution, the Q 2 dependence of F2(x,Q 2 ) is dictated by the gluon distribution 
as dF2(x, Q 2 )/d\ogQ 2 ~ a s xg(2x,Q 2 ) [65]. As discussed in Refs. [28,52], this leads to 
dR^ 2 (x,Q 2 )/dlogQ 2 ~ a s [R£(2x,Q 2 ) - R£ 2 (x,Q 2 )}xg(2x,Q 2 )/F 2 D (x,Q 2 ). The log Q 2 
slopes of F 2 Sn /F 2 c; measured by NMC [47] therefore constrain R G . Especially, the positive log Q 2 - 
slope of F| n /F 2 c measured by NMC indicates that R g {2x,Q 2 q) > R^ 2 \x,Q 2 q) at x ~ 0.01 
[28, 52]. Thus, within the DGLAP framework, much stronger gluon than antiquark shadowing at 
x ~ 0.01, as suggested e.g. in Ref. [66], is not supported by the NMC data. The antishadowing 
in the EKRS gluons follows from the constraint i?^(0.03, Qq) 1 imposed by the NMC data 
(see also Ref. [64]) combined with the requirement of momentum conservation. The EKRS an- 
tishadowing is consistent with the E789 data [67] on D-meson production in pA collisions (note, 
however, the large error bar on the data). It also seems to be supported by J/ip production in DIS, 
measured by NMC [68]. 

• At x < 0.02, stringent experimental constraints do not yet exist for the nuclear gluons. It should 
be emphasized, however, that the initial Rq in this region is directly connected to the initial 
Rp 2 . As discussed above, related to quarks at small x, taking the DIS data on Rp 2 in the non- 
perturbative region as a lower limit for Rp 2 at Qq corresponds to 3R f a / d log Q 2 > so that 
R G {x < 0.02, Q 2 ) > Rj, 2 {x/2,Q 2 ) > R^(x/2,Q 2 < lGeV 2 ). Reference [28] observed that 
setting Rq(x, Qq) Rp 2 (x, Qq) at x < 0.01 fulfills this constraint. This approximation remains 
fairly good even after DGLAP evolution from Q ~ 1 GeV to Q ~ 100 GeV, see Ref. [27]. 

As explained in detail in Ref. [28], in the EKRS approach the initial ratios Ry(x, Qq), R a (x, Qq) 
and Rq(x, Qq) are constructed piecewise in different regions of x. Initial nPDFs are computed at Qq. 
Leading order (LO) DGLAP evolution to higher scales is performed and comparisons with the DIS and 
DY data are made. The parameters in the input ratios are iteratively changed until the best global fit to the 
data is achieved. The determination of the input parameters in EKRS has so far been done only manually 
with the best overall fit determined by eye. For the quality of the fit, see the detailed comparison with the 
data in Figs. 4-10 of Ref. [28]. The parametrization EKS98 [27] of the nuclear modifications R A (x, Q 2 ) 
was prepared on the basis of the results from Ref. [28]. It was also shown that when the free proton PDFs 
were changed from GRVLO [69] to CTEQ4L [70] (differing from each other considerably), the changes 
induced in Rf{x, Q 2 ) were within a few percent [27]. Therefore, accepting this range of uncertainty, the 
EKS98 parametrization can be used together with any (LO) set of the free proton PDFs. 

The HKM analysis. In principle, the definition of the nPDFs in the HKM analysis [29] differs slightly 
from that in EKRS: instead of the PDFs of the bound protons, HKM define the nPDFs as the average 
distributions of each flavour Una nucleus A: ff{x,Q 2 ) = (Z/A)ff A {x,Q 2 ) + (l- Z/A)f™ /A (x, Q 2 ). 
Correspondingly, the HKM nuclear modifications at the initial scale Qq = 1 GeV 2 are then defined 
through 

f t A (x, Q 2 ) = Wi (x, A, Z)[(Z/A)ff(x, Ql) + (1 - Z/A)ft(x, Q 2 )}. (11) 

In practice, however, since in the EKS98 parametrization R$ v = Rj y and R% = Rj, the EKS98 mod- 
ifications also represent average modifications. Flavour-symmetric sea quark distributions are assumed 
in HKM whereas the flavour asymmetry of the sea quarks in EKRS follows from that of the free proton. 

An improvement relative to EKRS is that the HKM method to extract the initial modifications 
wi(x,A,Z) at Qq is more automatic and more quantitative in the statistical analysis since the HKM 



analysis is strictly a minimum-x 2 fit. Also, with certain assumptions of a suitable form (see Ref. [29]) 
for the initial modifications wi, the number of parameters has been conveniently reduced to seven. The 
form used in the fit is 

where Hi(x) = b{X + CiX 2 . An analysis with a cubic polynomial is also performed, with similar results. 
Due to the flavour symmetry, the sea quark parameters are identical for all flavours. For valence quarks, 
conservation of charge Z (not used in EKRS) and baryon number A are required, fixing a Uv and aa v . 
In the case of non-isoscalar nuclei w Uv / Wd v - Also, momentum conservation is imposed fixing a g . 
Taking (3y, by and cy to be equal for uy and dy, and (3q = (3 g = 1, b g = — 2c g , the remaining seven 
parameters by, cy, j3y, a g b g , cq and c g are determined by a global DGLAP fit which minimizes \ 2 - 



The comparison. Irrespective of whether the best fit is found automatically or by eye, the basic proce- 
dure to determine the nPDFs in the EKRS and HKM analyses is the same. The results obtained for the 
nuclear effects are, however, quite different, as can be seen in the comparison at Q 2 = 2.25 GeV 2 shown 
in Fig. 9. The main reason for this is that different data sets are used: 




Fig. 9: The nuclear modifications for the uv, u and gluon distributions in EKS98 (solid) and in HKM (dot-dashed) at Q 2 — 
2.25 GeV . Upper panels: C/D. Lower panels: Pb/D with isoscalar Pb. 



• The HKM analysis [29] uses only the DIS data whereas EKRS include also the DY data from pA 
collisions. As explained above, the DY data are very important in the EKRS analysis to fix the rela- 
tive modifications of valence and sea quarks at intermediate x. Preliminary results reported in Ref. 
[71] show that, when the DY data are included in the HKM analysis, the antiquark modifications 
become more similar to those of EKRS. 

• The NMC data on F^/F^ [61], which imposes quite stringent constraints on the ^4-systematics 
in the EKRS analysis, is not used in HKM. As a result, Rp 2 has less shadowing in HKM than in 
EKRS (see also Fig. 1 of Ref. [52]), especially for heavy nuclei. 

• In addition to the recent DIS data sets, some older ones are used in the HKM analysis. The older 
sets are not included in EKRS. This, however, presumably does not cause any major differences 
between EKRS and HKM since the older data sets have larger error bars and will therefore typi- 
cally have less statistical weight in the \ 2 analysis. 



• The HKM analysis does not make use of the NMC data [47] on the Q 2 dependence of F 2 Sn /F^. 
As explained above, these data are the main experimental constraint on the nuclear gluons in the 
EKRS analysis. Figure 10 shows a comparison of the EKRS (solid) and the HKM (dot-dashed) 
results to the NMC data on F 2 Sn /F 2 c as a function of Q 2 . The HKM results do not reproduce 
the measured Q 2 dependence of F^/F® at the smallest values of x. This figure demonstrates 
explicitly that the nuclear modifications of gluon distributions at 0.02 < x < 0.1 can be pinned 
down with the help of these NMC data. 
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Fig. 10: The ratio i<f n /F 2 c as a function of Q 2 at different fixed values of x. The data are from NMC [47] and the curves are 
EKS98 [27] (solid) and HKM [29] (dot-dashed). 



4.13 Future prospects 

Constraints from pA collisions at LHC, RHIC and SPS. The hard probes in pA collisions, especially 
those at the LHC, will play a major role in probing the nPDFs at regions not accessible before. Figure 8 
shows the x and Q 2 regions of the phase space probed by certain hard processes at different center of 
mass energies 1 . Let us consider the figure in more detail: 

• The measurements of semileptonic decays of D and D mesons in pA will help to pin down 
the nuclear gluon distributions [72,73]. Borrowed from a LO analysis [73], the thick solid lines in 
Fig. 8 show how the average scale Q 2 = (m^) of open charm production is correlated with the average 
fractional momentum x = (22) of the incoming nuclear gluon in dimuon production in correlated DD 
decays at the LHC (computed for pPb at y/s = 5500 GeV and single muon acceptance 2.5 < y^ < 4.0, 
with no rapidity shifts), at RHIC (pAu at y/s = 200 GeV, 1.15 < < 2.44) and at the SPS (pPb at 
y/s = 17.3 GeV, < y M < 1). For each solid curve, the smallest Q 2 shown corresponds to a dimuon 
invariant mass M = 1.25 GeV, and the the largest Q 2 to M = 9.5 GeV (LHC), 5.8 GeV (RHIC) and 
4.75 GeV (SPS) (cf. Fig. 5 in Ref. [73]). We observe that data from the SPS, RHIC and the LHC will 
offer constraints in different regions of x. The data from NA60 at the SPS will probe the antishadowing 
of nuclear gluons, while the RHIC data probe the onset of gluon shadowing, and the LHC data constrain 
the nuclear gluons at very small x, deep in the shadowing region. 

• Similarly, provided that the DY dimuons can be distinguished from those coming from the D and 
B decays, the DY dimuon pA cross sections are expected to set further constraints on the nuclear effects 

'No counting rates have been considered for the figure. 



on the sea quark distributions. In 2 — > 1 kinematics, X2 = (M/ ^fs)e~ y where M and y are the invariant 
mass and rapidity of the lepton pair. The rapidity of the pair is always between those of the leptonic 
constituents of the pair. The shaded regions in Fig. 8 illustrate the regions in x = x A and Q 2 = M 2 
probed by DY dimuons in pA collisions at the LHC, RHIC and SPS. The parameters used for the LHC 
are ^fs = 5500 GeV and 2.5 < y < 4.0 with no rapidity shifts [74]. Note that from the (phase space) 
point of view of x A (but not of x\), Pb+Pb collisions at \/*PbPb = 5500 GeV at the LHC are equivalent 
to pPb collisions at \fs pPh = 8800 GeV: the decrease of X2 due to the increase in yfs is compensated 

by the rapidity shift y [74]: x\ A = (M/Vs pA ) e ~ (j/AA+S/o) = ( M /\^AA) e ~ VAA = X 2 A - For RHIC, 
the shaded region corresponds to yfs = 200 GeV and 1.2 < y < 2.2, the rapidity acceptance of the 
PHENIX forward muon arm [75]. For the SPS, we have again used yfs = 17.3 GeV and < y < 1. 
Single muon p^-cuts (see Refs. [75, 76]) have not been applied and the shaded regions all correspond 
to M 2 > 1 GeV 2 . At the highest scales shown, the dimuon regions probed at the SPS and RHIC are 
limited by the phase space. The SPS data may shed more light on the existence of an EMC effect for the 
sea quarks [77]. At the LHC the dimuons from decays of D and B mesons are expected to dominate the 
dilepton continuum up to the Z mass [78-80]. If, however, the DY contribution can be identified (see 
Sec. 6.1 for more discussion), DY dimuons at the LHC will probe the sea quark distributions deep in the 
shadowing region and also at high scales, up to Q 2 = M|. The RHIC region again lies conveniently 
between the SPS and the LHC, mainly probing the beginning of sea quark shadowing region. In Ap 
collisions at the LHC, DY dimuons with y = (y m i n + y m ax)/2 probe the regions indicated by the dashed 
line (the spread due to the y acceptance is not shown). 

• The dotted line shows the kinematical region probed by open heavy flavor production in pA 
(and also in AA) collisions, when both heavy quarks are at y = in the detector frame. In this case, 
again assuming 2^2 kinematics, x v 2 A = x AA = 2Q/y/s AA with Q = mf. To illustrate the effect 
of moving towards forward rapidities, the case yQ = y-q = 3 is also shown. Only the region where 
Q 2 > m 2 (above the upper horizontal line) is probed by bb production. Due to the same kinematics (in 
X2 and Q 2 ) as in open heavy quark production, the dotted lines also correspond to the regions probed by 
direct photon production (with Q = pr). 

• As demonstrated by Fig. 8, hard probes in pA collisions at RHIC and at the LHC in particular 
will provide us with very important constraints on the nPDFs at scales where the DGLAP evolution 
is expected to be applicable. Power corrections to the evolution [25,26] can, however, be expected to 
play an increasingly important role towards small values of x and Q 2 . The effects of the first of such 
corrections, the GLRMQ terms [25, 26] leading to nonlinear DGLAP+GLRMQ evolution equations, 
have been studied recently in light of the HERA data [82]. From the point of view of the DGLAP 
evolution of the (n)PDFs, gluon saturation occurs when their evolution becomes dominated by power 
corrections [81]. Figure 8 shows the saturation limits obtained for the free proton and for Pb in the 
DGLAP+GLRMQ approach [32] (solid curves, the dotted curves are extrapolations to guide the eye). 
The saturation limit obtained for the free proton in the DGLAP+GLRMQ analysis should be taken as an 
upper limit in Q 2 . It is constrained quite well by the HERA data (see Ref. [82]). The constraints from 
the Q 2 dependence of F^/F^, however, were not taken into account in obtaining the saturation limit 
shown for the Pb nucleus [32]. For a comparison of the saturation limits obtained in other models, see 
Ref. [32]. 

• Figure 8 also shows which hard probes of pA collisions can be expected to probe the nPDFs 
in the gluon saturation region. The probes directly sensitive to the gluon distributions are especially 
interesting from this point of view. Such probes include open cc production at small px and direct photon 
production at px ~ few GeV, both as far forward in rapidity as possible (see the dotted lines). Note that 
open bb production is already in the region where linear DGLAP evolution is applicable. In light of the 
saturation limits shown, the chances of measuring the effects of nonlinearities in the evolution through 
open cc production in pA at RHIC would seem marginal. At the LHC, however, measuring saturation 
effects in the nuclear gluon distributions through open cc in pA could be possible. 



Improvements of the DGLAP analyses. On the practical side, the global DGLAP fits of the nPDFs 
discussed above can be improved in obvious ways. The EKRS analysis should be made more automatic 
and a proper statistical treatment, such as in HKM, should be added. The automatization alone is, how- 
ever, not expected to change the nuclear modifications of the PDFs significantly from EKS98 but more 
quantitative estimates of the uncertainties and of their propagation would be obtained. This work is in 
progress. As also discussed above, more contraints from data should be added to the HKM analysis. 
More generally, all presently available data from hard processes in DIS and pA collisions have not yet 
been exhausted: for instance, the recent CCFR DIS data for and F3 Fe from uFe and uFe collisions 
[83] (not used in EKRS or in HKM), could help to pin down the valence quark modifications [84]. In the 
future, hard probes of pA collisions at the LHC, RHIC and the SPS will offer very important constraints 
on the nPDFs, especially for gluons and sea quarks. Eventually, the nPDF DGLAP analyses should be 
extended to NLO perturbative QCD. As discussed above, the effects of power corrections [25, 26] to the 
DGLAP equations and also to the cross sections [85] should be analysed in detail in the context of global 
fits to the nuclear data. 

4.2 Nonlinear corrections to the DGLAP equations; looking for the saturation limits 

Kari J. Eskola, Heli Honkanen, Vesa J. Kolhinen, Jianwei Qiu, and Carlos A. Salgado 

Free proton PDFs, fi(x,Q 2 ), are needed for the calculation of hard process cross sections in 
hadronic collisions. Once they are determined at a certain initial scale Qq, the DGLAP equations [43] 
well describe their evolution to large scales. Based on global fits to the available data, several different 
sets of PDFs have been obtained [19, 20,44, 86]. The older PDFs do not adequately describe the recent 
HERA data [87] on the structure function F2 at the perturbative scales Q 2 at small x. In the analysis 
of newer PDF sets, such as CTEQ6 [19] and MRST2001 [44], these data have been taken into account. 
However, difficulties arise when fitting both small and large scale data simultaneously. In the new MRST 
set, the entire HI data set [87] has been used in the analysis, leading to a good average fit at all scales 
but at the expense of a negative NLO gluon distribution at small x and Q 2 < 1 GeV 2 . In the CTEQ6 set 
only the large scale (Q 2 > 4 GeV 2 ) data have been included, giving a good fit at large Q 2 but resulting 
in a poorer fit at small-x and small Q 2 (Q 2 < 4 GeV 2 ). Moreover, the gluon distribution at small x with 
Q 2 < 1.69 GeV 2 has been set to zero. 

These problems are interesting since they could indicate a new QCD phenomenon. At small values 
of momentum fraction x and scales Q 2 , gluon recombination terms, which lead to nonlinear corrections 
to the evolution equations, can become significant. The first of these nonlinear terms have been calculated 
by Gribov, Levin and Ryskin [25], and Mueller and Qiu [26]. In the following, these corrections shall 
be referred to as the GLRMQ terms for short. With these modifications, the evolution equations become 
[26] 

|^ 2 G< 2 W) + ...G„t, (14) 

where the two-gluon density can be modelled as x 2 G^ 2 \x, Q 2 ) = ■^ I [xg(x, Q 2 )] 2 , with proton radius 
R = 1 fm. The higher dimensional gluon term, Ght [26], is here assumed to be zero. The effects of the 
nonlinear corrections on the DGLAP evolution of the free proton PDFs were studied in Ref. [82] in view 
of the recent HI data. The results are discussed below. 

4.21 Constraints from the HERA data 

The goal of the analysis in Ref. [82] was (1) to possibly improve the (LO) fit of the calculated ^(x, Q 2 ) 
to the HI data [87] at small Q 2 while (2) concurrently maintaining the good fit at large Q 2 and, finally, 
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Fig. 11: Left: F 2 (x, Q 2 ) calculated using CTEQ6L [19] (dotted-dashed) and the DGLAP+GLRMQ results with 
set 1 (solid) and set 2a (double dashed) [82], compared with the HI data [87]. Right: The Q 2 dependence of 
the gluon distribution function at fixed x, from set 1 evolved with DGLAP+GLRMQ (solid) and directly from 
CTEQ6L (dotted-dashed). 



(3) to study the interdependence of the initial distributions and the evolution. 

In CTEQ6L, a good fit to the HI data is obtained (see Fig. 11) with a fiat small-x gluon distribution 
at Q 2 ~ 1.4 GeV 2 . As can be seen from Eqs. (13)-(14), the GLRMQ corrections slow the scale evolution. 
Now one may ask whether the HI data can be reproduced equally well with different initial conditions 
(i.e. assuming larger initial gluon distributions) and the GLRMQ corrections included in the evolution. 
This question was addressed in Ref. [82] by generating three new sets of initial distributions using 
DGLAP + GLRMQ evolved CTEQ5 [86] and CTEQ6 distributions as guidelines. The initial scale was 
chosen to be Qq = 1.4 GeV 2 , slightly below the smallest scale of the data. The modified distributions 
at Qq were constructed piecewise from CTEQ5L and CTEQ6L distributions evolved down from Q 2 = 3 
and 10 GeV 2 (CTEQ5L) and Q 2 = 5 GeV 2 (CTEQ6L). A power law form was used in the small-x region 
to tune the initial distributions until good agreement with the HI data was found. 

The difference between the three sets in Ref. [82] is that in set 1 there is still a nonzero charm 
distribution at Qq = 1.4 GeV 2 , slightly below the charm mass threshold with m c = 1.3 GeV in CTEQ6. 
In set 2a and set 2b the charm distribution has been removed at the initial scale and the resulting deficit 
in F2 has been compensated by slightly increasing the other sea quarks at small x. Moreover, the effect 
of charm was studied by using different mass thresholds: m c = 1.3 GeV in set 2a and m c = a/L4 GeV 
in set 2b, i.e. charm begins to evolve immediately above the initial scale. 

The results of the DGLAP+GLRMQ evolution with the new initial distributions are shown in 
Fig. 11. The left panel shows the comparison between the HI data and the (LO) structure function 
i*2(x, Q 2 ) = J2i e i x [Qi{ x i Q 2 ) + <li( x , Q 2 )] calculated from set 1 (solid lines), set 2a (double dashed) 



and the CTEQ6L parametrization (dotted-dashed lines). As can be seen, the results are very similar, 
showing that with modified initial conditions and DGLAP+GLRMQ evolution, one obtains an as good 
or even a better fit to the HERA data (x 2 /N = 1.13, 1.17, 0.88 for the sets 1, 2a, 2b, respectively) as 
with the CTEQ6L distributions (x 2 /N = 1-32). 

The evolution of the gluon distribution functions in the DGLAP+GLRMQ and DGLAP cases is 
illustrated more explicitly in the right panel of Fig. 11, where the absolute distributions for fixed x 
are plotted as a function of Q 2 for set 1 and for CTEQ6L. The figure shows how large differences at 
the initial scale vanish during the evolution due to the GLRMQ effects. At scales Q 2 > 4 GeV 2 , the 
GLRMQ corrections fade out rapidly and the DGLAP terms dominate the evolution. 



4.22 Saturation limits: p and Pb 

The DGLAP+GLRMQ approach also offers a way to study the gluon saturation limits. For each x 
in the small-x region, the saturation scale Q 2 at can be defined as the value of the scale Q 2 where the 
DGLAP and GLRMQ terms in the nonlinear evolution equation become equal, 9 gi^Qi ^ \ q^=q 2 t ( x ) = 0- 
The region of applicability of the DGLAP+GLRMQ is Q 2 > <2 2 at (x) where the linear DGLAP part 
dominates the evolution. In the saturation region, at Q 2 < Ql^x), the GLRMQ terms dominate, and all 
nonlinear terms become important. 

In order to find the saturation scale Q 2 at (x), for the free proton, the initial distributions (set 1) 
at Qq = 1.4 GeV 2 have to be evolved downwards in scale using the DGLAP+GLRMQ equations. As 
discussed in Ref. [82], since only the first correction term has been taken into account, the gluon distribu- 
tion near the saturation region should be considered as an upper limit. Consequently, the saturation scale 
obtained is an upper limit as well. The result is shown by the asterisks in Fig. 12. The saturation line 
for the free proton from the geometric saturation model by Golec-Biernat and Wusfhoff (G-BW) [88] is 
also plotted (dashed line) for comparison. It is interesting to note that although the DGLAP+GLRMQ 
and G-BW approaches are very different, the slopes are very similar at the smallest values of x. 
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Fig. 12: The gluon saturation limits in the DGLAP+GLRMQ approach for proton (asterisks) and Pb (A = 208), 
with (crosses) and without (dots) nuclear modifications [82]. The saturation line for the proton from the geometric 
saturation model [88] (dashed line), and for Pb from [89] (dotted-dashed) are also plotted. 



Saturation scales for nuclei can be determined similarly. For a nucleus A, the two-gluon density 
can be modelled as x 2 G^ 2 \x, Q 2 ) = -jp-[xgA(x, Q 2 )] 2 , i.e. the effect of the correction is enhanced by 



a factor of A 1 / 3 . An initial estimate of the saturation limit can be obtained by starting the downwards 
evolution at high enough scales, Q 2 = WO... 200 GeV 2 , for the GLRMQ terms to be negligible. Like 
in the proton case, the result is an upper limit. It is shown for Pb in Fig. 12 (dots). The effect of the 
nuclear modifications was also studied by applying the EKS98 [27] parametrization at the high starting 
scale. As a result, the saturation scales Q 2 at (x) are somewhat reduced, as shown in Fig. 12 (crosses). 
The saturation limit obtained for a Pb nucleus by Armesto in a Glauberized geometric saturation model 
[89] is shown (dotted-dashed) for comparison. Again, despite the differences between the approaches, 
the slopes are strikingly similar. 

For further studies and for more accurate estimates of Qg at (x) in the DGLAP+GLRMQ approach, 
a full global analysis of the nuclear parton distribution functions should be performed, along the same 
lines as in EKRS [27, 28] and HKM [29]. 

4.3 Gluon distributions in nuclei at small x: guidance from different models 

Nestor Armesto and Carlos A. Salgado 

The difference between the structure functions measured in nucleons and nuclei [46] is a very 
important and well known feature of nuclear structure and nuclear collisions. At small values of x 
(< 0.01, shadowing region), the structure function F2 per nucleon turns out to be smaller in nuclei 
than in a free nucleon. This corresponds to shadowing of parton densities in nuclei. While at small x 
valence quarks are of little importance and the behaviour of the sea is expected to follow that of F^, 
the gluon distribution, which is not an observable quantity, is badly determined and represents one of 
the largest uncertainties for cross sections both at moderate and large Q 2 in collinear factorization [17]. 
For example, the uncertainty in the determination of the gluon distributions for Pb at Q 2 ~ 5 GeV 2 at 
LHC for y = 0, x ~ mj-/ 'y/s ~ 10~ 4 . . . 10~ 3 , is a factor ~ 3 (see Fig. 13), which may result in an 
uncertainty of up to a factor of ~ 9 for the corresponding cross sections in PbPb collisions. 

In this situation, while waiting for new data from lepton-ion [90-92] or pA colliders, guidance 
from different theoretical models is of utmost importance for controlled extrapolations from the region 
where experimental data exist to those interesting for the LHC. Two different approaches to the problem 
have been suggested: On one hand, some models try to explain the origin of shadowing, usually in terms 
of multiple scattering (in the frame where the nucleus is at rest) or parton interactions in the nuclear wave 
function (in the frame in which the nucleus is moving fast). On the other hand, other models parameterize 
parton densities inside the nucleus at some scale large enough for perturbative QCD to be applied 
reliably and then evolve these parton densities using the DGLAP [43] evolution equations. Then, the 
origin of the differences between parton densities in nucleons and nuclei is not addressed, but contained 
in the parameterization at Qq, obtained from a fit to data. 

4.31 Multiple scattering and saturation models 

The nature of shadowing is well understood qualitatively. In the rest frame of the nucleus, the incoming 
photon splits, at small enough x, into a qq pair long before reaching the nucleus, with a coherence length 
l c oc l/(m,Nx) with tun the nucleon mass. At small enough x, l c becomes on the order of or greater 
than the nuclear size. Thus the qq pair interacts coherently with the nucleus with a typical hadronic cross 
section, resulting in absorption [93-100]. (See Ref. [101] for a simple geometrical approach in this 
framework.) In this way nuclear shadowing is a consequence of multiple scattering and is thus related to 
diffraction (see e.g. Refs. [102, 103]). 

Multiple scattering is usually formulated in the dipole model [97, 104], equivalent to fc^-factorization 
[105] at leading order. In this framework, the 7*-nucleus cross section is expressed as the convolution 
of the probability of the transverse or longitudinal 7* to split into a qq pair of transverse dimension r 
times the cross section for dipole-nucleus scattering. It is this dipole-nucleus cross section at fixed im- 
pact parameter which saturates, i.e. reaches the maximum value allowed by unitarity. Most often this 



multiple scattering problem is modelled through the Glauber-Gribov approach [89, 103, 106, 107]. The 
dipole-target cross section is related through a Bessel-Fourier transform to the so-called unintegrated 
gluon distribution ipa(x, kr) in kx -factorization [89, 106, 108], which in turn can be related to the usual 
collinear gluon density by 



where A 2 is some infrared cut-off, if required. This identification is only true for Q 2 S> Q 2 at [109, 110], 
where Q 2 &t corresponds to the transverse momentum scale at which the saturation of the dipole-target 
cross section occurs. 

Other formulations of multiple scattering do not use the dipole formalism but relate shadowing to 
diffraction by Gribov theory. See recent applications in [102, 103] and the discussion in Sec. 4.4. 

An explanation equivalent to multiple scattering, in the frame in which the nucleus is moving fast, 
is gluon recombination due to the overlap of the gluon clouds from different nucleons. This makes the 
gluon density in a nucleus with mass number A smaller than A times that of a free nucleon [25, 26]. At 
small x, the interaction develops over longitudinal distances z ^ \/{m^x) which become of the order of 
or larger than the nuclear size, leading to the overlap of gluon clouds from different nucleons located in a 
transverse area ~ l/Q 2 . Much recent work in this field stems from the development of the semiclassical 
approximations in QCD and the appearance of non-linear evolution equations in x (see Refs. [Ill, 112], 
the contributions by Mueller, Iancu et al, Venugopalan, Kaidalov and Kharzeev in [113], and references 
therein), although saturation appears to be different from shadowing [109, 110] (i.e. the reduction in the 
number of gluons as defined in DIS is not apparent). In the semiclassical framework, the gluon field at 
saturation reaches a maximum value and becomes proportional to the inverse QCD coupling constant, 
gluon correlations are absent and the dipole-target cross section has a form which leads to geometrical 
scaling of the 7* -target cross section. Such scaling has been found in small x nucleon data [114] and 
also in analytical and numerical solutions of the non-linear equations [115-118] but the data indicate 
that apparently the region where this scaling should be seen in lepton-nucleus collisions has not yet been 
reached [119]. Apart from proposed eA colliders [90-92], the LHC will be the place to look for non- 
linear effects. However, the consequences of such high density configurations may be masked in AB 
collisions by other effects, such as final-state interactions. Thus pA collisions would be essential for 
these studies, see e.g. Ref. [120]. Moreover, they would be required to fix the baseline for other studies 
in AB collisions such as the search for and characterization of the QGR As a final comment, we note 
that the non-linear terms in Refs. [25, 26] are of a higher-twist nature. Then, in the low density limit, the 
DGLAP equations [43] are recovered. However, the non-linear evolution equations in x [111-113] do 
not correspond to any definite twist. Their linear limit corresponds to the BFKL equation [121, 122]. 

We now comment further on the difference between gluon shadowing and saturation. In saturation 
models, see Refs. [111-113], saturation, defined as either a maximum value of the gluon field or by the 
scattering becoming black, and gluon shadowing, defined as xGa/{xGn) < 1, are apparently different 
phenomena, i.e. saturation does not necessarily lead to gluon shadowing [109, 110]. Indeed, in numerical 
studies of the non-linear small x evolution equations, the unintegrated gluon distribution turns out to be 
a universal function of just one variable r = k^/Q 2 ^ [115-117], vanishing quickly for k\ > Q 2 at - 
(In Ref. [118] this universality is analytically shown to be fulfilled up to k\ much larger than Q 2 at ). 
This scaling also arises in the analysis of DIS data on nucleons [114]. It has also been searched for 
in nuclear data [119]. In nuclei, Q 2 &t increases with increasing nuclear size, centrality and energy. 
Thus, through the relation to the collinear gluon density given in Eq. (15), this scaling implies that 
the integral gives the same value (up to logarithmic corrections if a perturbative tail oc \jk\ exists 
for k\ 3> Q 2 a t), provided that the upper limit of the integration domain Q 2 ^> Q 2 at . Thus the ratio 
xGa/{xGn) remains equal to 1 (or approaching 1 as a ratio of logs if the perturbative tail exists) for 
Q 2 » Qs a f So saturation or non-linear evolution do not automatically guarantee the existence of gluon 
shadowing. Since shadowing is described in terms of multiple scattering in many models, this result 




(15) 



[109, 110] seems somewhat surprising. 



4.32 DGLAP evolution models 

A different approach is taken in Refs. [27-29, 123]. Parton densities inside the nucleus are parameterized 
at some scale Qq ~ 1 ... 5 GeV 2 and then evolved using the DGLAP [43] equations. Now all nuclear 
effects on the parton densities are included in the parameterization at Qq, obtained from a fit to experi- 
mental data. The differences between the approaches mainly come from the data used to constrain the 
parton distributions (e.g. including or excluding Drell-Yan data). In these calculations lack of data leaves 
the gluon badly constrained at very small x. Data on the (Revolution of give direct constraints [52] 
through DGLAP evolution but only for x > 0.01. See Sec. 4. 1 for more discussion of this framework. 

4.33 Comparison between different approaches 

Model predictions usually depend on additional semi-phenomenological assumptions. Since these as- 
sumptions are typically different, the predictions often contradict each other. For example, in Refs. 
[93-97] it is argued that large size qq configurations give the dominant contribution to the absorption, 
which does not correspond to any definite twist but to an admixture of all twists, resulting in essen- 
tially Q 2 -independent shadowing. On the other hand, in the gluon recombination approach [25, 26] the 
absorption is clearly a higher-twist effect dying out at large Q 2 . Finally, in the DGLAP-based models 
[27-29, 123], all Q 2 -dependence comes from QCD evolution and is thus of a logarithmic, leading-twist 
nature. Thus these approaches all lead to very different predictions at small x, particularly for the gluon 
density. 
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Fig. 13: Ratios of gluon distribution functions from different models at Q 2 = 5 GeV 2 ; HKM refers to the results from Ref. 
[29], Sarcevic, Ref. [107], EKS98, Refs. [27, 28], Frankfurt, Ref. [103], Armesto, Refs. [89, 106] and new HIJING, Ref. [66]. 
The bands represent the ranges of x = (Q/y r s)e y for processes with \y\ < 0.5, Q 2 = 5 GeV 2 at RHIC (y/5 = 200 GeV) and 
LHC (Vs = 5.5 TeV). 

We now compare available numerical results from different approaches. Recall that in dipole ap- 
proaches such as that of Ref. [89], there are difficulties in identifying the unintegrated gluon distribution 
with the ordinary gluon density at small and moderate Q 2 , see the discussion of Eq. (15). A comparison 



at Q 2 = 5 GeV 2 for the Pb/p gluon ratio can be found in Fig. 13. In the RHIC region, x ~ 10~ 2 , 
the results from Refs. [27-29, 89, 107] roughly coincide but lie above those from Refs. [66, 103]. At 
x ~ 10~ 4 . . . 10~ 5 (accessible at the LHC) the results of Ref. [89] are smaller than those of Refs. 
[27-29, 107], approaching those of Ref. [103] but remaining larger than those of Ref. [66]. Apart from 
the loose small x gluon constraints from DIS data on nuclei, in Refs. [27-29] saturation is imposed on 
the initial conditions for DGLAP evolution, while Refs. [89, 103, 107] are multiple scattering models. 
In Ref. [66], the Q 2 -independent gluon density is tuned at x ~ 10~ 2 to reproduce charged particle 
multiplicities in Au+Au collisions at RHIC. The strongest gluon shadowing is thus obtained. However, 
it seems to be in disagreement with existing DIS data [52]. Additional caution has to be taken when 
comparing results from multiple scattering models with those from DGLAP analysis [27-29]. The gluon 
ratios at some moderate, fixed Q 2 and very small x may become smaller (e.g. in Ref. [89]) than the F2 
ratios at the same x and Q 2 , which might lead to problems with leading-twist DGLAP evolution [52]. 

4.4 Predictions for the leading-twist shadowing at the LHC 

Leonid Frankfurt, Vadim Guzey, Martin McDermott, and Mark Strikman 

The basis of our present understanding of nuclear shadowing in high-energy scattering on nuclei 
was formulated in the seminal work by V. Gribov [124]. The key observation was that, within the 
approximation that the range of the strong interactions is much smaller than the average inter-nucleon 
distances in nuclei, there is a direct relationship between nuclear shadowing in the total hadron-nucleus 
cross section and the cross section for diffractive hadron-nucleon scattering. While the original derivation 
was presented for hadron-deuterium scattering, it can be straightforwardly generalized to lepton-nucleus 
DIS. 
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Fig. 14: Diagrams for hard diffraction in ep scattering and for the leading-twist nuclear shadowing. 



Reference [126] demonstrated how to use the Collins factorization theorem for hard diffraction in 
DIS [127] to generalize Gribov theory to calculate the leading-twist component of nuclear shadowing 
separately for each nuclear parton distribution. The correspondence between the two phenomena is 
illustrated in Fig. 14. 

As a result, the double scattering term of the shadowing correction to the nuclear parton distribu- 



tions, 5f^ A (x, Q 2 ), takes the form 
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Here j is a generic parton label (i.e. a gluon or a quark of a particular flavor), fP, N is the diffractive partem 
distribution function (DPDF) of the nucleon for parton j, p A (b, z) is the nucleon density normalized to 
unity, i] is the ratio of the real to imaginary parts of the nucleon diffractive amplitude, while xjp^ = 
0.1 for quarks and xjp : o = 0.03 for gluons. Since ffj N obeys the leading-twist DGLAP evolution 

equation, the Q 2 -evolution of Sfff A is also governed by DGLAP, i.e. it is by definition a leading-twist 
contribution. This explains why the approach of Ref. [126] can be legitimately called a leading-twist 
approach. However, in this approach, one can also take into account higher-twist effects in diffraction. 

The HERA diffractive data are consistent with the dominance of leading-twist diffraction in the 
process 7* + p — > Mx + p for Q 2 > 4 GeV 2 . These studies have determined NLO quark DPDFs 
directly and the gluon diffractive DPDF indirectly through scaling violations. More recent studies of 
charm production and diffractive dijet production in the direct photon kinematics measured the gluon 
DPDF directly. It was found to be consistent (within a 30% uncertainty primarily determined by the 
treatment of the NLO contribution to two jet production) with the inclusive diffraction in DIS. As a 
result, we were able to determine the strength of interaction in the quark and gluon channels, a 3 eS . This 
strength turns out to be large at Qq ~ 4 GeV 2 , where one can set up the boundary condition 2 for the 
subsequent QCD evolution. The interaction strength is typically of the order of 20 mb for the quark 
channel and 30-40 mb for the gluon channel at Q 2 ~ 4 GeV 2 , x < 10~ 3 . Consequently, it is necessary 
to include the interactions with N > 3 nucleons in the calculations. These interactions were modeled 
using the quasi-eikonal approximation. Corrections due to fluctuations in the cross section were also 
analyzed and found to be rather small for reasonable models of these effects [128]. Therefore, the master 
equation for the evaluation of nuclear shadowing takes the following form 



It is important to note that Eq. (17) defines the input nuclear PDFs at the initial, low scale Q = Qo = 
2 GeV for subsequent QCD evolution to higher scales Q 2 . In order to have consistent results after 
evolution, Eq. (17), applicable for 10~ 5 < x < 0.05, should be supplemented by a model of nuclear 
PDFs for larger x, x > 0.05. In our analysis, this is done by assuming no enhancement for quarks and 
some enhancement (antishadowing) for gluons (see Fig. 15). 

We have recently performed a detailed analysis of nuclear shadowing within the leading-twist 
approach using the current HERA diffractive data [103]. We have analyzed the uncertainties originat- 
ing from the input diffractive parton distribution functions, related to the uncertainties in the data, and 
found that they are less than < 20 . . . 30% for x < 10~ 3 . The biggest uncertainty in the gluon chan- 
nel originates from the t-dependence of the gluon DPDF, not directly measured. An example of our 
calculations [129] is presented in Fig. 15. 

2 For Q 2 < 4 GeV 2 , higher-twist effects in ep diffraction appear to be significant. This leads to 30 . . . 50% contributions 
of higher-twist effects in nuclear shadowing at Q 2 ~ 1 ... 2 GeV 2 , which may make the use of the small-a; NMC data for the 
extraction of nuclear PDFs problematic. 
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Fig. 15: Predictions of quark and gluon shadowing for Q = 2 GeV (solid curves), Q — 10 GeV (dashed curves) and Q = 100 
GeV (dot-dashed curves). 



We also calculated the higher-twist shadowing, which we estimated in the eikonal approximation, 
and demonstrated that the leading-twist contribution dominates down to very small x at Q 2 > 4 GeV 2 . 
(Numerical results are available from V. Guzey upon request.) Also, we found large differences between 
the LO and NLO calculations of the nuclear structure function 

One can see from Fig. 15 that the leading-twist shadowing effects are expected to be large for 
small x for both the quark and gluon channels at Q ~ few GeV. These effects are expected to survive up 
to scales characteristic of Z and W production. 

Our formalism can be readily applied to evaluate nuclear shadowing in nPDFs at all impact pa- 
rameters. The nuclear shadowing correction to the impact-parameter-dependent nPDFs, Sfj /a{ x , Q 2 \ b), 
can be found from Eq. (17) by simply differentiating with respect to the impact parameter b [103]. The 
effect of nuclear shadowing increases with decreasing impact parameter. This is illustrated in Fig. 16, 
where we plot the ratios fj/A{x, Q 2 , b)/(ATA(b)fj/N(x, Q 2 )), where Ta(6) = f dzpA(b, z), for quarks 
and gluons at b = (solid curves) and b = 6 fm (dashed curves). For comparison, the impact-parameter- 
averaged ratios fj/A(x, Q 2 ) / (Afj/ N (x, Q 2 )) are given by the dot-dashed curves. Figure 16 demonstrates 
that nuclear shadowing significantly depends on the impact parameter b. Hence, comparing peripheral 
and central pA collisions will provide an additional measure of the nuclear shadowing. 

It is worth emphasizing that the diagrams corresponding to leading-twist shadowing are usually 
neglected in the BFKL-type approaches to nuclear scattering. For example, these diagrams are neglected 
in the model of Balitsky [130] and Kovchegov [131]. Note also that higher-twist effects become impor- 
tant at sufficiently small x and Q < Q m i n (x). Rough estimates indicate that, in the gluon sector with 
A ~ 200, Qmin(10~ 3 ) ~ 2 ... 5 GeV, depending on how the NLO effects are handled. Also, Q m \ n (x) 
should increase with decreasing x. 

In conclusion, leading-twist shadowing effects are expected to lead to significant modifications of 




Fig. 16: Predictions of quark and gluon shadowing as a function of the impact parameter: b = (solid curves); b = 6 rm 
(dashed curves); fr-averaged ratios (dot-dashed curves). For all curves, Q = 2 GeV. 



the hard process cross sections at the LHC in pA and AA collisions. Further progress in the diffraction 
studies at HERA will significantly reduce the uncertainties in the predictions. 



4.5 Structure functions are not parton distributions 

Stanley J. Brodsky 

Ever since the earliest days of the parton model, it has been assumed that the leading-twist struc- 
ture functions Fi(x, Q 2 ) measured in deep inelastic lepton scattering are determined by the probability 
distributions of quarks and gluons as determined by the light-cone (LC) wavefunctions of the target. For 
example, the quark distribution is 
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The identification of structure functions with the square of light-cone wavefunctions is usually made in 
the LC gauge, n • A = A + = 0, where the path-ordered exponential in the operator product for the 
forward virtual Compton amplitude apparently reduces to unity. Thus the deep inelastic lepton scatter- 
ing cross section (DIS) appears to be fully determined by the probability distribution of partons in the 
target. However, Paul Hoyer, Nils Marchal, Stephane Peigne, Francesco Sannino, and I have recently 
shown that the leading-twist contribution to DIS is affected by diffractive rescattering of a quark in the 
target, a coherent effect which is not included in the light-cone wavefunctions, even in light-cone gauge 
[132-134] . The distinction between structure functions and parton probabilities is already implied by 
the Glauber-Gribov picture of nuclear shadowing [125, 135-137]. In this framework shadowing arises 
from interference between complex rescattering amplitudes involving on-shell intermediate states, as in 



Fig. 17. In contrast, the wavefunction of a stable target is strictly real since it does not have on energy- 
shell configurations. A probabilistic interpretation of the DIS cross section is thus precluded. 

It is well-known that in the Feynman and other covariant gauges one has to evaluate the correc- 
tions to the "handbag" diagram due to the final-state interactions of the struck quark (the line carrying 
momentum p\ in Fig. 18) with the gauge field of the target. In light-cone gauge, this effect also involves 
rescattering of a spectator quark, the p2 line in Fig. 18. The light-cone gauge is singular - in particu- 
lar, the gluon propagator d l ^ c (k) = k >2 % +ie [— g^ u + ^^n-k^™" ] nas a pole at k + = which requires 
an analytic prescription. In final-state scattering involving on-shell intermediate states, the exchanged 
momentum k + is of O (1/V) in the target rest frame, which enhances the second term in the propagator. 
This enhancement allows rescattering to contribute at leading twist even in LC gauge. 




Fig. 17: Glauber-Gribov shadowing involves interference between rescattering amplitudes. 



The issues involving final-state interactions even occur in the simple framework of abelian gauge 
theory with scalar quarks. Consider a frame with q + < 0. We can then distinguish FSI from ISI using 
LC time-ordered perturbation theory, LCPTH [138]. Figure 18 illustrates two LCPTH diagrams which 
contribute to the forward j*T — ► j*T amplitude, where the target T is taken to be a single quark. In 
the aligned jet kinematics the virtual photon fluctuates into a qqpair with limited transverse momentum, 
and the (struck) quark takes nearly all the longitudinal momentum of the photon. The initial q and q 
momenta are denoted p\ and p2 — h\, respectively. 
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Fig. 18: Two types of final state interactions, (a) Scattering of the antiquark (p2 line), which in the aligned jet kinematics is part 
of the target dynamics, (b) Scattering of the current quark (pi line). For each LC time-ordered diagram, the potentially on-shell 
intermediate states - corresponding to the zeroes of the denominators D a , Db, D c - are denoted by dashed lines. 



The calculation of the rescattering effects on DIS in Feynman and light-cone gauge through three 
loops is given in detail in Ref. [132]. The result can be resummed and is most easily expressed in 
eikonal form in terms of transverse distances tt, Rt conjugate to p2T> for- The DIS cross section can be 



expressed as 
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is the resummed result. The Born amplitude is 

A(p 2 ,r T ,R T ) = 2eg 2 MQp 2 V (m\\r T )W (f T , R T ) 
where mf, = p 2 Mx B + rn 2 and 



f d 2 p T e i?T 'VT 1 

V{mr T ) = / 77^2^— 5 = ^-^o(mr T ). 

7 (2it) z pL + 27r 



(18) 



(19) 



(20) 



(21) 



The rescattering effect of the dipole of the qq is controlled by 



VF(r r ,i? T )= / 77T e l 



(2tt) 2 Jfc2 2vr S I R T 



(22) 



The fact that the coefficient of A in Eq. (19) is less than unity for all fx, Rt shows that the rescattering 
corrections reduce the cross section. It is the analog of nuclear shadowing in our model. 

We have also found the same result for the DIS cross sections in light-cone gauge. Three prescrip- 
tions for defining the propagator pole at k + = have been used in the literature: 
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(23) 



the principal-value (PV), Kovchegov (K) [139], and Mandelstam-Leibbrandt (ML) [140] prescriptions. 
The 'sign function' is denoted e(x) = @(x) — 6(— x). With the PV prescription we have I v = 
f dk 2 [k^],^ = 0. Since an individual diagram may contain pole terms ~ ^/kf , its value can depend 

on the prescription used for light-cone gauge. However, the kf = poles cancel when all diagrams 
are added. The net is thus prescription-independent and agrees with the Feynman gauge result. It is 
interesting to note that the diagrams involving rescattering of the struck quark p\ do not contribute to the 
leading-twist structure functions if we use the Kovchegov prescription to define the light-cone gauge. In 
other prescriptions for light-cone gauge the rescattering of the struck quark line p\ leads to an infrared 
divergent phase factor exp(z</>), where 
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(24) 



where A is an infrared regulator, and I v = 1 in the K prescription. The phase is exactly compensated 
by an equal and opposite phase from FSI of line pi. This irrelevant change of phase can be understood 
by the fact that the different prescriptions are related by a residual gauge transformation proportional to 
5(k + ) which leaves the light-cone gauge A + = condition unaffected. 

Diffractive contributions which leave the target intact thus contribute at leading twist to deep 
inelastic scattering. These contributions do not resolve the quark structure of the target, and thus they are 
contributions to structure functions which are not partem probabilities. More generally, the rescattering 
contributions shadow and modify the observed inelastic contributions to DIS. 



Our analysis in the K prescription for light-cone gauge resembles the "covariant parton model" 
of Landshoff, Polkinghorne and Short [141, 142] when interpreted in the target rest frame. In this de- 
scription of small x DIS, the virtual photon with positive q + first splits into the pair p\ and p2- The 
aligned quark p\ has no final state interactions. However, the antiquark line p2 can interact in the target 
with an effective energy s oc k\jx while staying close to mass shell. Thus at small x and large s, the 
antiquark p2 line can first multiply scatter in the target via pomeron and Reggeon exchange, and then it 
can finally scatter inelastically or be annihilated. The DIS cross section can thus be written as an integral 
of the (Tqp->x cross section over the p2 virtuality. In this way, the shadowing of the antiquark in the 
nucleus g^a-^x cross section yields the nuclear shadowing of DIS [136]. Our analysis, when interpreted 
in frames with q + > 0, also supports the color dipole description of deep inelastic lepton scattering at 
small x. Even in the case of the aligned jet configurations, one can understand DIS as due to the coherent 
color gauge interactions of the incoming quark-pair state of the photon interacting first coherently and 
finally incoherently in the target. For further discussion see Refs. [133, 143, 144]. The same final-state 
interactions which produce leading-twist diffraction and shadowing in DIS also lead to Bjorken-scaling 
single-spin asymmetries in semi-inclusive deep inelastic, see Refs. [145-147]. 

The analysis presented here has important implications for the interpretation of the nuclear struc- 
ture functions measured in deep inelastic lepton scattering. Since leading-twist nuclear shadowing is 
due to the destructive interference of diffractive processes arising from final-state interactions (in the 
q + < frame), the physics of shadowing is not contained in the wavefunctions of the target alone. For 
example, the light-front wavefunctions of stable states computed in light-cone gauge (PV prescription) 
are real, and they only sum the interactions within the bound-state which occur up to the light-front time 
r = when the current interacts. Thus the shadowing of nuclear structure functions is due to the mutual 
interactions of the virtual photon and the target nucleon, not the nucleus in isolation [132]. 

4.6 Nuclear parton distribution functions in the saturation regime 

Yuri V. Kovchegov 

4.61 Introduction 

At the high collision energies reached at the LHC, parton densities in the nuclear wavefunctions will 
become very large, resulting in the onset of parton recombination and saturation of parton distribution 
functions [23,25,26]. In the saturation regime, the growth of partonic structure functions with energy 
is expected to slow down sufficiently to unitarize the total nuclear cross sections. Gluonic fields in the 
saturated nuclear wavefunction would become very strong, reaching the maximum possible magnitude 
in QCD of ~ 1/g [139, 148, 149]. The transition to the saturation region can be characterized by 
the saturation scale Q 2 (s, A), which is related to the typical two dimensional density of the parton color 
charge in the infinite momentum frame of the hadronic or nuclear wavefunctions [23, 139, 148, 149]. The 
saturation scale is an increasing function of energy yfs and of the atomic number A [116, 118, 131, 150- 
154]. It is expected that it (roughly) scales as 

Q 2 s {s,A) ~ A a s 5 (25) 

where 8 » 0.2 - 0.3 based on HERA data [88, 155] and a > 1/3 [116, 118, 152-154]. Various 
estimates predict the saturation scale for lead nuclei at 14 TeV to be Q 2 S « 10 — 20 GeV 2 . This high 
value of the saturation scale would ensure that the strong coupling constant a s (Q 2 ) is small, making the 
vast majority of partons in the nuclear wave functions at the LHC perturbative [109, 110, 149, 156, 157]! 
The presence of an intrinsic large momentum scale Q s justifies the use of a perturbative QCD expansion 
even for such traditionally nonperturbative observables as total hadronic cross sections [23, 150]. 

Below we first discuss the predictions for gluon and quark distribution functions of a large nucleus, 
given by the quasi-classical model of McLerran and Venugopalan [23, 109, 139, 148, 150]. We then 
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Fig. 19: The non-Abelian Weizsacker- Williams field of a nucleus. 



discuss how quantum corrections should be included to obtain an equation for the realistic F2 structure 
function of the nucleus [130, 131, 151, 158, 159]. We will conclude by deriving the gluon production 
cross section for pA scattering in the quasi-classical approximation [109] and generalizing it to include 
quantum corrections [160]. 



4.62 Structure functions in the quasi-classical approximation 
The gluon distribution 

As suggested originally by McLerran and Venugopalan [23], the gluon field of a large ultrarelativistic 
nucleus is classical and is given by the solution of the Yang-Mills equations of motion with the nucleus 
providing the source current. The classical solution in the light cone gauge of the nucleus, known as 
the non-Abelian Weizsacker- Williams field (WW), has been found in Refs. [139, 148, 149]. It includes 
the effects of all multiple rescatterings of the small-x gluon on the nucleons in the nucleus located at 
the same impact parameter as the gluon. This resummation parametrically corresponds to resumming 
powers of Q 2 s /q 2 (higher twists), where q is the gluon transverse momentum. Diagrammatically, the 
field is shown in Fig. 19. The correlator of two WW fields gives rise to the unintegrated nuclear gluon 
distribution given by [109, 149] 



cIxGa(x, q 2 



dq 2 (2 



^ J d 2 ze-^- f d 2 b_Tr(A ww (0)A ww (z)) 
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where b is the gluon impact parameter and 



Q 2 s (z) = ^ 2 a f[ c P xG N {x, l/z 2 ) T A {b) (27) 



is the saturation scale with atomic number density p = A/[(4/3)ttR\] and = 2y R 2 A — b 2 for 

spherical nuclei A considered here. The gluon distribution in a single nucleon, xGn(x, Q 2 ), is assumed 
to be a slowly varying function of Q 2 . In the two-gluon exchange approximation [109] xGn{x, Q 2 ) = 
' 2j ^r L 111 O 2 \ 2 . ' 

The gluon distribution function of Eq. (26) falls off perturbatively as xGa ~ Q 2 s /q 2 f° r large 
transverse momenta q, but is much less singular in the infrared region, scaling as xGa ~ InQs/q 
for small q [149]. We conclude that multiple rescattering softens the infrared singular behavior of the 
traditional perturbative distribution. Note also that the number of gluons given by Eq. (26) is of the order 
of l/a s , corresponding to strong gluonic fields of the order of A^ ~ 1/g. 



Fig. 20: DIS on a nucleus in the quasi-classical approximation. 



The quark distribution 

To calculate the quark distribution function in the quasi-classical approximation, one has to consider 
deep inelastic scattering in the rest frame of the nucleus including all multiple rescatterings, as shown in 
Fig. 20. The F 2 structure function (and the quark distribution) was first calculated by Mueller [150] in 
this approximation, yielding 

F 2 (x,Q 2 ) = Y,4[xq(x,Q 2 )+xq(x,Q 2 )] = ^—J2 e l J __<F*-^(r,z) d 2 b N (r, b, Y), 

(28) 

where $ 7 ~~* qq (r, z) is the wavef unction of a virtual photon in DIS, splitting into a qq pair with transverse 
separation r. The fraction of photon longitudinal momentum carried by the quark is z. The quantity 
N(r,b, Y) is the forward scattering amplitude of a dipole with transverse size r at impact parameter b 
with rapidity Y on a target nucleus. In the quasi-classical approximation of Refs. [23, 150], N(r, b,Y = 
0) is given by 

N(r,b,Y = 0) = l - e - r2Q li^ /4 (29) 

where Q 2 quark is the quark saturation scale which can be obtained from the gluon saturation scale in Eq. 
(27) by replacing N c by Cp in the numerator [150, 160]. 



4.63 Inclusion of quantum corrections 

At high energy, the quasi-classical approximation previously considered ceases to be valid. The energy 
dependence of the cross sections comes in through the quantum corrections which resum logarithms of 
energy [122, 161-165]. The evolution equation for N, resumming logarithms of energy, closes only in 
the large-iV c limit of QCD and reads [104, 130, 131, 151, 158, 159, 166-168] 

N(x ol ,b,Y) = N(x ol ,b,0)e * 'H r ) Y + dye - ln i - ) {Y y) 

n Jo 



L 



x 2 1 11 

d 2 x 2 2 01 2 [2N(x 02 ,b + -x 12 ,y) - N(x 02 ,b + -x 12 ,y) N(x 12 ,b + -x 02 ,y)], (30) 

x 02 a: '12 ^ Z Z 



with initial condition in Eq. (29) and x^ = x j — Xy In the usual Feynman diagram formalism, it resums 
multiple pomeron exchanges (fan diagrams). Together with Eq. (28), it provides the quark distribution 
function in a nucleus at high energies. The equation has been solved both by approximate analytical 
methods and numerically in Refs. [116, 118, 152-154], providing unitarization of the total DIS cross 
section and saturation of the quark distribution functions in the nucleus. 

Using Eq. (30) together with Eq. (26), one may conjecture the following expression for the unin- 
tegrated gluon distribution function including quantum corrections [160] 

dxG A (x,q 2 ) 2 f , 2j , 2 _,-„.„ C F 



dq 2 



j d 2 bd 2 ze-^-^N G {z_,b_,Y = lnl/x) (31) 




with Nq the forward amplitude of a gluon dipole scattering on the nucleus which, in the large-iV c limit, 
is Nq = 2N — N 2 [160], where N is the forward amplitude of the qq dipole scattering on the nucleus 
used in Eq. (30). Equation (31) must still be proven. 

4.64 Gluon production in pA collisions 

The pA program at the LHC will test the above predictions for parton distribution functions by measuring 
particle production cross sections which depend on these quantities. Below we discuss gluon production 
in pA. 

Quasi-classical result 

Gluon production in pA interactions has been calculated in the quasi-classical approximation in Ref. 
[109]. The relevant diagrams in the light cone gauge of the proton are shown in Fig. 21. The cross 
section is [109] 

4^ = f dHd'xd'y^f^^e^y) (l-e-^^-e-y-^ + e-^ 2 ^ 
d A k dy J tt 2 {2ir) z xfy z \ 

(32) 

Gluon production including quantum corrections 

To include quantum corrections in Eq. (32), it must be rewritten in the factorized form [160] and gen- 
eralized to include the linear evolution of the proton structure function and nonlinear evolution of the 
nuclear structure function. The result, see also Ref. [160], is 

^-2 —,.2 

<larA ( '' r d*zn G (z,Y-y) V * e ** V * N G (z,y), (33) 
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where Y is the total rapidity interval between the proton and the nucleus, Nq{z_, y) = J d 2 bNc{z, b, y) 

-*2 

and V 2 is the transverse gradient squared. The function nc(z_,Y — y) in Eq. (33) is the solution of 
the linear part of Eq. (30) (BFKL) integrated over all b with the two-gluon exchange amplitude between 
the proton and a gluon dipole as initial conditions. In the no-evolution limit, Eq. (33) reduces to Eq. 
(32). Equation (33) gives the saturation physics prediction for particle production at midrapidity in pA 
collisions at the LHC. 



5. THE BENCHMARK TESTS IN pA 



5.1 W, Z, and Drell-Yan cross sections 

Ramona Vogt 

Studies of W ± and Z° gauge bosons as well as high mass Drell-Yan production in pA interactions 
at the LHC will provide information on the nuclear parton densities at much higher Q 2 than previously 
available. The high Q 2 measurements will probe the x range around 0.015 at y = 0. This range of x has 
been probed before by fixed target nuclear deep-inelastic scattering so that a measurement at the LHC 
will constrain the evolution of the nuclear quark distributions over a large lever arm in Q 2 . Knowledge 
of the nuclear quark distributions in this region will aid in studies of jet quenching with a Z° opposite 
a jet in AA collisions. At least one of the initial partons will be a quark or antiquark in the reactions 
qq -> Z°g and q(q)g -> q(q)Z°. 

Since the W and Z° studies will open a new window in the (x, Q 2 ) plane for nuclear parton distri- 
bution functions, these measurements are obvious benchmarks. The Z° measurement probes symmetric 
qq annihilation while measurements such as the W + — W~ asymmetry via high p F decay leptons are 
better probes of the nuclear valence quark distributions. Drell-Yan production at high Q 2 is a somewhat 
different story because, to be an effective benchmark, the Drell-Yan signal should dominate the dilepton 
mass spectrum. In AA collisions, the signal dilepton mass distribution is dominated by bb decays until 
M > Mz [169]. Thus only in this large mass region is Drell-Yan a useful benchmark. 

We very briefly outline the next-to-leading order, NLO, calculation here. TheNNLO cross sections 
are available [170] but are only a small addition (a few percent) compared to the NLO correction. The 
NLO cross section for production of a vector boson, V, with mass M at scale Q in a pp interaction is 
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where is proportional to the LO partonic cross section, ij — ► V, and the sum over Q includes u, 
d, s and c. The matrices C n and C lf contain information on the coupling of the various quark flavors 
to boson V. The parton densities in the proton, ff, are evaluated at momentum fraction x and scale 
Q 2 . The convolutions of the parton densities, including the proton and neutron (isospin) content of the 
nucleus, with shadowing parameterizations formulated assuming ff(x, Q 2 ) = ff(x, Q 2 )Rf(x, Q 2 )/A 
are given in the appendix of Ref. [171] 3 . All the results are calculated with the MRST HO densities [172] 
at Q = My and the EKS98 shadowing parameterization [27,28]. In the Drell-Yan process, the mass 
distribution can also be calculated by adding a dM in the denominator of the left-hand side of Eq. (34) 
and the delta function S(M — My ) on the right-hand side. The pref actors Uj- are rather simple [170], 



z o 8ir G F , 2 ; 2 Mf w ± 2vr G f 



where G F = 1.16639 x 10" 5 GeV 2 , M z = 91.187 GeV, and M w = 80.41 GeV. In the case of the Drell- 
Yan process, there are three contributions to : virtual photon exchange, Z° exchange, and 7* — Z° 

3 InRef. [171], S\x,Q 2 ,A) = Rf(x,Q 2 ). 
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where sin 2 O w = 1 - M^/M|, T z = 2.495 GeV, and T z ^ l+l - = 3.367%. The functions A^O, Q 2 ) 
in Eq. (34) are universal for gauge bosons and Drell-Yan production [170]. We work in the MS scheme. 
The NLO correction to the qq channel includes contributions from soft and virtual gluons as well as hard 
gluons from the process qq — > Vg. The quark-gluon contribution only appears at NLO through the real 
correction qg — > qV. At NLO, a s (Q 2 ) is calculated to two loops with nf = 5. 

As we show in Fig. 22, the isospin of the nucleus is most important for (/(/-dominated processes. 
Gauge boson and Drell-Yan production are the best examples of this domination in hadroproduction. 
We compare the pPb and Pbp W + , W~ and Z° distributions per nucleon with and without the EKS98 
shadowing parameterization to pp production at 5.5 TeV. The pp distributions closely follow the pPb 
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Fig. 22: The W + , W~ and Z° rapidity distributions in pp, pPb and Pbp collisions at 5.5 TeV/nucleon evaluated at Q = My- 
The solid and dashed curves show the results without and with shadowing respectively in Pbp collisions while the dotted and 
dot-dashed curves give the results without and with shadowing for pPb collisions. The dot-dot-dot-dashed curve is the pp result. 



distributions without shadowing except for y < 2. Since the W + distribution is most sensitive to the 
fu(xi)f^(x2) contribution, the pp and pPb distributions peak away from y = 0. The Pbp distribution, 
on the other hand, decreases with increasing y because the neutron content of lead washes out the peak. 
The opposite effect is seen for W~ production where, in Pbp, the neutrons produce a peak at forward 
rapidity since f% = fu- The Z° distribution is almost independent of isospin, as are the high mass 



Drell-Yan distributions. We see antishadowing in Pbp interactions with large nuclear x\ while we have 
small nuclear X2 in the pPb calculations, reducing the y distributions. The W + pA distributions are all 
reduced compared to the pp distributions, the W~ pA results are typically enhanced relative to pp while 
the Z° pA/pp ratios do not change much with y. 

The gauge boson cross sections in pp, Pbp, and pPb collisions at 5.5 TeV/nucleon as well as Pbp 
and pPb collisions at 8.8 TeV/nucleon are given in Table 5. The Drell-Yan cross sections in the mass 



Table 5: Gauge boson cross sections per nucleon. No decay channel is specified. 



Ay 


ons (nb) 


OEKS98 (nb) 


cjns (nb) 


OEKS98 (nb) 


<tns (nb) 


CEKS98 (nb) 






7^ 


W + 












pp 5.5 TeV 








< y < 2.4 


8.10 


- 


15.13 


- 


11.41 


- 


< y < 1 


3.51 


- 


6.13 


- 


5.16 


- 


2.4 < y < 4 


2.01 


- 


5.68 


- 


2.42 


- 








Pbp 5.5 TeV 








< y < 2.4 


8.23 


8.37 


13.12 


13.34 


13.37 


13.63 


< y < 1 


3.53 


3.53 


5.71 


5.67 


5.56 


5.54 


2.4 < y < 4 


2.14 


2.06 


3.72 


3.62 


4.36 


4.24 








pPb 5.5 TeV 








< y < 2.4 


8.12 


7.39 


14.87 


13.43 


11.62 


10.57 


< y < 1 


3.52 


3.33 


5.95 


5.60 


5.31 


5.01 


2.4 < y < 4 


2.01 


1.68 


5.67 


4.69 


2.43 


2.02 








Pbp 8.8 TeV 








< y < 2.4 


12.37 


12.47 


19.62 


19.71 


19.87 


20.01 


< y < 1 


5.22 


5.11 


8.37 


8.16 


8.22 


8.02 


2.4 < y < 4 


5.07 


5.03 


8.30 


8.31 


9.46 


9.47 








pPb 8.8 TeV 








< y < 2.4 


12.26 


10.93 


21.41 


18.96 


18.08 


16.09 


< y < 1 


5.20 


4.82 


8.60 


7.92 


7.98 


7.35 


2.4 < y < 4 


4.83 


3.95 


11.80 


9.58 


5.98 


4.87 



bins 60 < M < 80 GeV, 80 < M < 100 GeV, and 100 < M < 120 GeV are given in Table 6 for pp, 
Pbp, and pPb collisions at 5.5 TeV/nucleon. The cross sections are presented in three rapidity intervals, 
< y < 2.4, corresponding to the CMS barrel+endcaps, < y < 1, the central part of ALICE, and 
2.4 < y < 4, the ALICE muon arm. The difference in the Pbp and pPb cross sections without shadowing 
is entirely due to isospin. Thus the pA cross section in \y\ < 2.4 is not a simple factor of two larger than 
the < y < 2.4 cross section since the results are not symmetric around y = 0, as can be seen in Fig. 22. 

We have checked the dependence of the gauge boson results on scale, parton density, and shadow- 
ing parameterization. The scale dependence enters in the PDFs, a s , and the scale-dependent logarithms 
ln(Q 2 /M£) [170]. At these high Q 2 values, M v /2 < Q < 2M V , the change in a s is rs ±10%, small 
compared to the variation at lower scales. When Q 2 = My, the scale dependent terms drop out. If 
Q 2 > M v , the logs are positive, enhancing the 0(a s ) component, but if Q 2 < M 2 , the scale-dependent 
logs change sign, reducing the NLO contribution. Evolution with Q 2 increases the low x density of the 
sea quarks and gluons while depleting the high x component and generally reducing the valence distri- 
butions. Since the x values do not change when Q 2 is varied, the higher scales also tend to enhance the 
cross sections. The gauge boson cross sections are decreased ~ 30% when Q 2 is decreased and increased 
by 20% with increased Q 2 . The scale change is even smaller for high mass Drell-Yan production, a few 
percent in the 60-80 GeV bin and 1% or less for higher masses. 



Table 6: Drell-Yan cross sections in the indicated mass intervals at 5.5 TeV/nucleon in pp, Pbp and pPb collisions. 



Ay 


cr(nb) 


<7ns (nb) 0-EKS98 (nb) 


<7ns (nb) 


CEKS98 (nb) 




pp 


Pbp 




pPb 






60 < M < 80 GeV 






< y < 2.4 


0.012 


0.012 0.012 


0.012 


0.011 


< y < 1 


0.0048 


0.0048 0.0047 


0.0048 


0.0044 


2.4 < y < 4 


0.0038 


0.0036 0.0035 


0.0038 


0.0031 






80 < M < 100 GeV 






< y < 2.4 


0.26 


0.27 0.27 


0.26 


0.24 


< y < 1 


0.11 


0.11 0.11 


0.11 


0.10 


2.4 < y < 4 


0.055 


0.059 0.057 


0.055 


0.046 






100 < M < 120 GeV 






< y < 2.4 


0.0088 


0.0088 0.0090 


0.0088 


0.0080 


< y < 1 


0.0037 


0.0037 0.0037 


0.0037 


0.0035 


2.4 < y < 4 


0.0015 


0.0015 0.0014 


0.0015 


0.0012 



The variation with parton density can be as large as the scale dependence. The MRST and CTEQ 
[86] rapidity distributions are similar but the CTEQ cross sections are 7% higher. The GRV 98 [45] 
results are stronger functions of rapidity than the MRST. The total GRV 98 cross sections, however, are 
about 13% smaller. 

We have also compared the pA/pp ratios to the HPC [173] and HKM [29] parameterizations. 
The HPC parameterization does not distinguish between the sea and valence effects and includes no Q 2 
evolution. HKM actually parameterizes the nuclear parton distributions themselves and focuses on the 
distinction between the up and down nuclear distributions. It is only useful to LO. However, shadowing 
at LO and NLO is nearly identical [174]. 

In Pbp collisions, the EKS98 and HKM ratios are similar for y < 1.5. At higher rapidities, the 
HKM ratio is larger because its sea distribution increases more with x\ but its valence shadowing is 
weaker. The HPC ratio is lower than the EKS98 ratio at y ~ 0, ~ 0.8 relative to ~ 0.97. However 
the two ratios are very similar for y > 2. In pPb collisions, x 2 is in the shadowing region and all three 
ratios decrease with rapidity. Again, near y ~ 0, the HKM and EKS98 ratios are similar but the HKM 
sea quark shadowing is nearly independent of x 2 while the EKS98 ratio drops from 0.97 to 0.8 between 
y = and 4. The HPC ratio is parallel to but lower than the EKS98 ratio. 

Finally, we note that, without keeping the energy fixed, extracting the nuclear quark distributions 
from pp and pA and applying them to AA is difficult 4 since the isospin is an additional, nonnegligible 
uncertainty. Therefore, it is desirable to make pp and pPb runs at 5.5 TeV for a more precise correlation 
with the Pb+Pb data at 5.5 TeV/nucleon. Indeed, for best results, both pPb and Pbp runs should be 
done since the ALICE muon detector is not symmetric around y = 0. Its large rapidity coverage would 
extend the nuclear x coverage to much higher (Pbp) and lower (pPb) values since y ~ 4 is the maximum 
rapidity for Z° production at 5.5 TeV/nucleon, the upper limit of the muon coverage, where x\ ~ 1 and 
x 2 ~ 0.0003. 

Systematic studies with nuclear beams could be very useful for fully mapping out the A depen- 
dence of the nuclear parton densities. The LHC Ar beam would link to the NMC Ca data while the LHC 
O beam could connect back to previous C data [59]. However, these A systematics are best tuned to the 
AA collision energies, different for O, Ar, and Pb. Comparing pA and pp at the same energies does not 

4 For a discussion of shadowing effects on gauge boson production in Pb+Pb collisions at 5.5 TeV/nucleon and the possible 
impact parameter dependence of the effect, see Ref. [171]. 



shift x while a comparison to pp at 14 TeV rather than 5.5 TeV would shift the x values by a factor of 
~ 3 between Pb+Pb and pp. Energy systematics of gauge boson production between RHIC and LHC 
are more difficult because W and Z° production in pA at RHIC is quite close to threshold. However, at 
RHIC, larger x values can be probed at y ~ 0, x ~ 0.45, where the cross sections are quite sensitive to 
the valence distributions. High mass Drell-Yan is likely to be beyond the reach of RHIC. 



5.2 Z° and W transverse momentum spectra 

Xiaofei Zhang and George Fai 

At LHC energies, perturbative QCD (pQCD) provides a powerful calculational tool. For Z° and 
W transverse momentum spectra, pQCD theory agrees with the CDF [175] and DO [176] data very well 
at Tevatron energies [177]. The LHC pp program will test pQCD at an unprecedented energy. The 
heavy-ion program at the LHC will make it possible for the first time to observe the full pt spectra 
of heavy vector bosons in nuclear collisions and will provide a testing ground for pQCD resummation 
theory [178]. 

In nuclear collisions, power corrections will be enhanced by initial and final state multiple scatter- 
ing. As we will show, higher-twist effects are small at the LHC for heavy boson production. The only 
important nuclear effect is the nuclear modification of the PDFs (shadowing). Because of their large 
masses, the W and Z° will tell us about the nPDFs at large scales. Since contributions from different 
scales need to be resummed, the p? spectra of heavy vector bosons will provide information about the 
evolution of nPDFs from small scales to large scales. Since it is more difficult to detect W ± than Z°, we 
will concentrate on the Z° here (the results for W ± production are very similar [179]). 

Resummation of large logarithms in QCD can be carried out either directly in p^-space, or in the 
so-called "impact parameter", 6-space, the Fourier conjugate of the p^-space. Using the renormalization 
group equation technique, Collins, Soper, and Sterman (CSS) derived a formalism for the transverse 
momentum distribution of vector boson production in hadronic collisions [178]. In the CSS formalism, 
nonperturbative input is needed for the large b region. The dependence of the pQCD results on the 
nonperturbative input is not weak if the original extrapolation proposed by CSS is used. Recently, a new 
extrapolation scheme was introduced, based on solving the renormalization group equation including 
power corrections [177]. Using the new extrapolation formula, the dependence of the pQCD results on 
the nonperturbative input was significantly reduced. 

For vector boson (V) production in a hadron collision, the CSS resummation formalism yields 

[178] 

MhA l£%£ + X) = I ^ni,M,, A , XB) + Y^,M M , XB) , (37, 

where xa = e y M/^/s and xb = e~ y M/y/s, with rapidity y and collision energy y/s. In Eq. (37), 
the W term dominates the pr distributions when pr <C M, and the Y term gives corrections that are 
negligible for small px but become important when px ~ M. 

The function W(b, M, xa,xb) can be calculated perturbatively for small b, but an extrapolation 
to the large b region, requiring nonperturbative input, is necessary in order to complete the Fourier 
transform in Eq. (37). In order to improve the situation, a new form was proposed [177] by solving 
the renormalization equation including power corrections. In the new formalism, W(b, M, xa,xb) = 
W pcrt {b, M, xa, xb), when b < 6 max , with 

^P crt (6, M, x A , x B ) = e 5 ^ M ) W(b, c/b, x A , x B ) , (38) 

where all large logarithms from ln(l/6 2 ) to ln(M 2 ) have been completely resummed into the exponential 
factor S(b, M), and c is a constant of order unity [178]. For b > b max , 

W(b, M, x A , x B ) = VF pcrt (6max, M, x A , x B )F NP (b; 6 max ) , (39) 



where the nonperturbative function F NP is given by 



F = exp{- ln(M 2 6Lx/c 2 ) [gi ((b 2 ) a - (&Lx) a j +92 [b 2 - b 2 max )\ - g 2 (b 2 - b 2 ^)}. (40) 

Here, b max is a parameter that separates the perturbatively calculated part from the nonperturbative input. 
Unlike in the original CSS formalism, W(b, M, xa, xb) is independent of the nonperturbative param- 
eters when b < b max . In addition, the fe-dependence in Eq. (40) is separated according to different 
physics origins. The (b ) a -dependence mimics the summation of the perturbatively calculable leading- 
power contributions to the renormalization group equations to all orders in the running coupling constant 
a s (/j,). The 6 2 -dependence of the g 2 term is a direct consequence of dynamical power corrections to the 
renormalization group equations and has an explicit dependence on M. The g 2 term represents the effect 
of the non-vanishing intrinsic parton transverse momentum. 

A remarkable feature of the 6-space resummation formalism is that the resummed exponential 
factor exp[S'(6, M)\ suppresses the 6-integral when b is larger than 1/M. It can be shown using the 
saddle point method that, for a large enough M, QCD perturbation theory is valid even at px = [178]. 
As discussed in Refs. [177, 179], the value of the saddle point strongly depends on the collision energy 
y/s, in addition to its well-known M 2 dependence. The predictive power of the 6-space resummation 
formalism should be even better at the LHC than at the Tevatron. 

In Z° production, since final-state interactions are negligible, power corrections can arise only 
from initial state multiple scattering. Equations (39) and (40) represent the most general form of W and 
thus, apart from isospin and shadowing effects, which will be discussed later. The only way nuclear 
modifications associated with scale evolution enter the W term is through the coefficient g 2 . 

The parameters g\ and a of Eq. (40) are fixed by the requirement of continuity of the func- 
tion W(b) and its derivative at b = 6 max - (The results are insensitive to changes of 6 max € [0.3 
GeV _1 ,0.7 GeV -1 ]. We use 6 max = 0.5 GeV _1 .) A fit to low energy Drell-Yan data gives g 2 = 
0.25 ± 0.05 GeV 2 and g 2 = 0.01 ± 0.005 GeV 2 . As the b dependence of the g 2 and g 2 terms in Eq. (40) 
is identical, it is convenient to combine these terms and define G 2 = ln(M 2 b 2 aax /c 2 )g 2 + g 2 . Using 
the values of the parameters listed above, we get G 2 = 0.33 ± 0.07 GeV 2 for Z° production in pp col- 
lisions. The parameter G 2 can be considered the only free parameter of the nonperturbative input in Eq. 
(40), arising from the power corrections in the renormalization group equations. An impression of the 
importance of the power corrections can be obtained by comparing results with the above value of G 2 to 
those with power corrections turned off (G 2 = 0). We therefore define the ratio 

dpT I dpT 

The cross sections in the above equation and in the results presented here have been integrated over 
rapidity (—2.4 < y < 2.4) and M 2 in the narrow width approximation. For the PDFs, we use the 
CTEQ5M set [86]. 

Figure 23 displays the differential cross sections and the corresponding Rg 2 ratio (with the limit- 
ing values of G 2 = 0.26 GeV 2 (dashed) and G 2 = 0.40 GeV 2 (solid)) for Z° production as functions of 
Pt in pp collisions at yfs = 14 TeV The deviation of Rg 2 from unity decreases rapidly as px increases. 
It is smaller than one percent for both yfs = 5.5 TeV (not shown) and y/s = 14 TeV in pp collisions, 
even when pt = 0. In other words, the effect of power corrections is very small at the LHC over the 
whole pt region. 

Without nuclear effects on the hard collision, the production of heavy vector bosons in nucleus- 
nucleus (AB) collisions should scale with the number of binary collisions, and gab oc AB. However, 
there are several additional nuclear effects on the hard collision in a heavy-ion reaction. First of all, the 
isospin effects, which come from the difference between the neutron PDFs and the proton PDFs, is about 
2% at LHC. This is expected, since at the LHC x ~ 0.02 where the u — d asymmetry is very small. 
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Fig. 23: (a) Cross section da/dpr for Z production in pp collisions at the LHC with 
Eq. (41) with G 2 = 0.26 GeV 2 (dashed) and 0.40 GeV 2 (solid), (c) Cross section in pPb collisions at ^/s = 8.8 TeV; (d) Ra 2 
in pPb with G 2 = 0.8GeV 2 . 

The dynamical power corrections entering the parameter g 2 should be enhanced by the nuclear 
size, i.e. proportional to A 1 / 3 . Taking into account the A-dependence, we obtain G 2 = 1.15 ± 0.35 
GeV 2 for PbPb reactions. We find that with this larger value of G 2 , the effects of power corrections 
appear to be enhanced by a factor of about three from pp to PbPb collisions at the LHC. Thus, even the 
enhanced power corrections remain under 1 % when 3 GeV < pt < 80 GeV. This small effect is taken 
into account in the following nuclear calculations. 

Next we turn to the phenomenon of shadowing, expected to be a function of x, the scale fi, and 
of the position in the nucleus. The latter dependence means that in heavy-ion collisions, shadowing 
should be impact parameter (b) dependent. Here we concentrate on impact-parameter integrated results, 
where the effect of the 6-dependence of shadowing is relatively unimportant [180], and focus on the scale 
dependence. We therefore use the EKS98 parametrization [27]. We define 



da^ (p T , Z A /A, Z B /B) I da(p 7 



dpi 



dpi 



(42) 



where Za and Zb are the atomic numbers and A and B are the mass numbers of the colliding nuclei. 
The cross section da^ (px, Za/A, Zb/B)/ dpr has been averaged over AB, while da(pr) / dpr is the 
pp cross section. We have seen above that shadowing remains the only significant effect responsible for 
nuclear modifications. Figure 23(c) shows da^ jdpT and (d) the corresponding Rq 2 for pPb collisions 
at = 8.8 TeV. 

In Fig. 24(a), the result for i? s h (solid line) is surprising because even at px = 90 GeV, x ~ 0.05, 
and we are still in the "strict shadowing" region. Therefore, the fact that > 1 for 20 GeV < pr < 
70 GeV is not "anti-shadowing". To better understand the shape of the ratio as a function of px, we also 
show i? s h with the scale fixed to 5 GeV (dashed line) and 90 GeV (dotted), respectively, in Fig. 24(a). 
The nuclear modification to the PDFs is only a function of x and flavor in the calculations represented 
by the dashed and dotted lines. These two curves are similar in shape, but rather different from the solid 
line. In b space, W(b, M, xa,xb) is almost equally suppressed in the whole b region if the fixed scale 
shadowing is used. However, with scale-dependent shadowing, the suppression increases as b increases 
since the scale \i ~ 1/6 in the nPDFs. Thus the scale dependence redistributes the shadowing effect. 
Here the redistribution brings above unity for 20 GeV < px < 70 GeV. When px increases further, 
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Fig. 24: (a) Cross section ratios for Z production in PbPb collisions at yfs — 5.5 TeV: i? s h of Eq. (42) (solid line), and R B h 
with the scale fixed at 5 GeV (dashed) and 90 GeV (dotted); (b) R Bh in pPb at y/s = 8.8 TeV. 

the contribution from the Y term starts to be important and i? s h dips back below one to match the fixed 
order pQCD result. Figure 24(b) shows in pPb at y/s = 8.8 TeV. 

We see from Fig. 24 that the shadowing effects on the pr distribution of Z° bosons at the LHC are 
intimately related to the scale dependence of the nPDFs where the data are very limited [27]. Theoretical 
studies such as EKS98 are based on the assumption that the nPDFs differ from the parton distributions 
in the free proton but obey the same DGLAP evolution [27]. Therefore, the transverse momentum 
distribution of heavy bosons at the LHC in PbPb and pPb collisions can further test our understanding of 
the nPDFs. 

In summary, higher-twist nuclear effects appear to be negligible in Z° production at LHC energies. 
We have demonstrated that the scale dependence of shadowing may lead to unexpected phenomenology 
of shadowing at these energies. Overall, the Z° transverse momentum distributions can be used as 
a precision test for leading-twist pQCD in the TeV energy region for both proton-proton and nuclear 
collisions. We propose that measurements of the Z° spectra be of very high priority at the LHC. 



5.3 Jet and dijet rates in pA collisions 

Alberto Accardi and Nestor Armesto 

The study of jet production has become one of the precision tests of QCD. Next-to-leading order 
(NLO) computations have been successfully confronted with experimental data. For example, jet pro- 
duction at the Tevatron [181] provides a very stringent test of available QCD Monte Carlo codes and 
gives very valuable information on the behaviour of the gluon distribution in the proton at large x. 

On the other hand, few studies on jets in a nuclear environment have been performed in fixed 
target experiments [182, 183] and at colliders. There are prospects to measure high-p^ particles at RHIC 
but the opportunity of detecting and studying jets with Et ~ 100 GeV, where Et is the total transverse 
energy of all jets, will appear only at the LHC. Concretely, such measurements in pA collisions would 
settle some uncertainties such as the validity of collinear factorization [17, 184, 185], the modifications 
of nPDFs inside nuclei [46, 186], and the existence and amount of energy loss of fast partons inside cold 
nuclear matter [187]. Such issues have to be well understood before any claim, based on jet studies, of 




the existence of new physics in nucleus-nucleus collisions (hopefully QGP formation), can be considered 
conclusive. Furthermore, the study of high-pr particles and jets can help not only to disentangle the 
existence of such new state of matter but also to characterize its properties, see e.g. [188-190]. 

For our computations we will use a NLO Monte Carlo code [191-193], adapted to include isospin 
effects and modifications of the PDFs inside nuclei. This code is based on the subtraction method to 
cancel the infrared singularities between real and virtual contributions. A full explanation, a list of 
available codes and detailed discussions on theoretical uncertainties in nucleon-nucleon collisions can 
be found in Ref. [194]. The accuracy of our computations, limited by CPU time, is estimated to be: 

• 2 % for the lowest and 15 % for the highest £/r-bins of the transverse energy distributions. 

• 3 % for the pseudorapidity distributions. 

• 20 % for the least populated and 3 % for the most populated bins of the dijet distributions of the 
angle between the two hardest jets. 

The results will be presented in the LHC lab frame, so that for pPb collisions we have a 7 TeV proton 
beam on a 2.75 TeV Pb beam, see Sec. 2.1 on the experimental parameters in pA collisions. All the 
energies will be given per nucleon and, in order to compare with the pp case, all cross sections will be 
presented per nucleon-nucleon pair, i.e. divided by AB. 

Unless explicitly stated, we will use the MRST98 central gluon distribution [172] with the EKS98 
parameterization [27,28] and identical factorization and renormalization scales n = \xp = [1r = Et/2. 
We will employ the kp -clustering algorithm for jet reconstruction [195, 196] with D = 1 which is more 
sound on theoretical grounds than the cone algorithm [197-199]. The kinematical regions we consider 
are: 

• |?7i| < 2.5, where rji is the pseudorapidity of the jet, corresponding to the acceptance of the central 
part of the CMS detector. 

• Eti > 20 GeV for a single jet in the pseudorapidity distributions, to ensure the validity of pertur- 
bative QCD. 

• Eti > 20 GeV and Et2 > 15 GeV for the hardest and next-to-hardest jets entering the CMS 
acceptance in the (^-distributions, where 4> is the angle between these two jets. 

Please note that, even in the absence of nuclear effects, the 7^ -distributions may be asymmetric with 
respect to rji = since we present our results in the lab frame, not the center-of-mass frame. 

Asymmetric Exi-cuts have been imposed on the 4> distributions to avoid the generation of large 
logarithms at certain points in phase space, see Ref. [193]. Also, the results near <j> = it are unreliable 
[193], because they require all-order resummation, not available in the code we are using. 

We caution that a comparison of our results with data is not straightforward. Comparison with 
AB data will be even more problematic, see the jet Chapter for details. Several physical effects are 
not included in our computations. For example, the so-called underlying event (i.e. the soft particle 
production that coexists with the hard process) is not considered here. It may cause difficulties with jet 
reconstruction and increase the uncertainties in jet-definition algorithms. This effect has been considered 
in antiproton-proton collisions at the Tevatron [200, 201] but its estimate is model dependent, relying 
on our limited knowledge of soft multiparticle production. In collisions involving nuclei, the situation 
is even worse since our knowledge is more limited. The only way to take this into account would 
be to use available Monte Carlo simulations for the full pA event [202,203], modified to include the 
NLO contributions. Other processes which might need to be taken into account more carefully in our 
computations are multiple hard parton collisions [204-206]. They may be divided into two classes: 

1) Disconnected collisions: more than one independent parton-parton collisions in the same event, 
each one producing a pair of high-p^ jets. This has been studied at the Tevatron [207, 208] for 
single pp collisions, but its quantitative explanation and the extension to nuclear collisions is model 
dependent. Simple estimates of the influence of disconnected collisions on jet production in pA 



show it to be negligible 5 . This is not the case for AB collisions, see the jet Chapter of this report. 

2) Rescatterings: a high-p^ partem may undergo several hard collisions before hadronizing into a jet, 
mimicking a single higher-order process. Quantitative descriptions of the modification of the jet 
Exi spectrum are model dependent. However, at the LHC the effect is expected to be very small, 
almost negligible at the transverse energies considered here (see Sec. 6.5 on the Cronin effect in 
pA collisions). 

Finally, no centrality dependence has been studied. It is not clear how to implement such dependence in 
theoretical computations because both the number of nucleon-nucleon collisions and the modification of 
the nPDF may be centrality dependent. 

5.31 da/ cIEt and da / dr/ for jets at large Et 
Uncertainties 

A number of theoretical uncertainties in pp and pPb collisions are examined in Figs. 25 and 26. Varying 
the scale between £r/4 and Et gives differences of order ±10 %, the smaller scale giving the larger 
result. We have estimated the uncertainty due to the PDF by using CTEQ5M [86]. It gives a ~ 5 % 
larger results than the default MRST98 central gluon [172]. The effect of isospin alone (obtained from 
the comparison of pp and pPb without any modification of nPDF at the same energy per nucleon) is very 
small, as one would expect from gg dominance. The effect of nPDFs, estimated using EKS98 [27,28], 
is also small. However it produces an extra asymmetry (an excess in the region m < 0) in pA collisions 
at reduced energy (e.g. pPb at 5.5 TeV), which disappears at maximum energy (e.g. pPb at 8.8 TeV), 
see Fig. 26 (bottom-left) and Fig. 28 (bottom-right), due to the asymmetry of the projectile and target 
momentum in the LHC lab frame. Finally, there is some uncertainty due to the choice of jet-finding 
algorithm. The cone algorithm [197-199] with R = 0.7 slightly reduces the cross section ~ 1 ... 2 % 
compared to the fer-clustering algorithm [195, 196] with D = 1. The cone with R = 1 gives results 
~ 15 % larger than our default choice. 

Finally, we discuss the ratio of the NLO to LO cross sections, the iv'-factor. In Fig. 27 the iv'-factor 
is calculated for different collision systems and energies, varying the inputs. For our default options, the 
ratio is quite constant, ~ 1.2, over all En and i]i for all energies. The ^-dependence results in K ~ 1.35 
for ii = Et and ~ 1.15 for p = £t/4. The cone algorithm with R = 0.7 gives results very similar 
to those obtained with the /cp-clustering algorithm with D = 1, while a cone with R = 1 produces 
K ~ 1.45. (Note that at LO the choice of jet-finding algorithm has no influence on the results.) Finally, 
the choice of PDF, isospin and modifications of the nPDFs do not sizably influence the iv'-factor. 

Results 

The expected LHC luminosities for different collisions are shown in Table 7, along with the one month 
(10 6 s) integrated luminosity expressed in units of fib/(AB). Using this table and the cross sections for 
inclusive single jet and dijet production in the figures, it is possible to extract the number of expected jets 

5 We make a simple estimate of the average number of binary nucleon-nucleon collisions (n) associated with jet production 
for En > Eto- In the Glauber model [209, 210], one obtains: 

(n)(b,E T0 ) ^-(Eto) = AT A (b)a pp (E T0 ) , 

where b is the impact parameter and the nuclear profile function Ta(&) is normalized to unity. The cross sections are a(Ero) = 
o- pp (Eto) = J^° T dEn dcr pp (ET)/dETi, the jet cross section for jets with En > Eto in pp collisions and a p A(b, Eto) = 
da p A(ETo)/d 2 b — 1— [1 — TA(b)a pp (ETo)] A , the inelastic jet production cross section for Eti > Eto in pA collisions as 
a function of b. In both cases the cross sections are integrated over the pseudorapidity acceptance. Taking ct(Eto) = 527, 1.13 
and 0.046 /Jo as representative values in pPb collisions at 8.8 TeV for Eto = 20, 100 and 200 GeV respectively, see Fig. 28, 
(n) = 1.0 for all Eto, both for minimum bias collisions (i.e. integrating Ta(6) and da p A(b, ETo)/d 2 b over b between b = 
and oo) and for central collisions (e.g. integrating < b < 1 fm). Thus the contribution of multiple hard scattering coming 
from different nucleon-nucleon collisions seems to be negligible. 
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Fig. 25: Left-hand side: Scale dependence of jet cross sections [we use /i = Et/2 (solid), /i = £t /4 (dashed) and n = £t 
(dotted) lines] as a function of En for r/i < 2.5 (upper plot) and rji for En > 20 GeV (lower plot). Right-hand side: PDF 
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(lower lines, scaled by 0.1 in the upper plots and by 0.25 in the lower plots) are shown. Unless otherwise stated, default options 
are used, see text. 
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From top to bottom: scale dependence (p = St/2: solid; [i — £V/4: dashed; [i = Et'- dotted); dependence on the jet 
reconstruction algorithm (fcT-algorithm with D = 1: solid; cone algorithm with R = 0.7: dashed; cone algorithm with 
R = 1: dotted); nucleon PDF dependence (MRST98 central gluon: solid; CTEQ5M: dashed lines); and isospin and nuclear 
PDF dependence (pp: solid; pA without: dashed; pA with EKS98: dotted). Six upper plots: pp at 14, 6.3 (scaled by 0.5) and 
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Table 7: Luminosities L (cm _2 s _1 ) and L x 10 6 s (fib/AB) for different collision systems at the LHC. The numbers of 
expected jets and dijets in a certain kinematical range are obtained by multiplying the last column by the cross sections given 
in Figs. 25, 26, 28 (jets) and 29, 30 (dijets). 
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(Figs. 25, 26 and 28) or dijets (Figs. 29 and 30). For example, from the solid curve on the upper left-hand 
side of Fig. 28 and C = 10 34 cm _2 s _1 , we expect 10 10 jets with Et% ~ 70 GeV from a cross section of 
1 (jb/(AB), and 10 6 jets with E Ti ~ 380 GeV from a cross section of 10~ 4 jjbl(AB), m our rji range. 

Looking at the results of Fig. 28, it is evident that samples of ~ 10 3 jets are feasible, up to 
~ 325 GeV in one month run at the given luminosity. Indeed, from Table 7, 10 3 jets in pPb at 8.8 
TeV would correspond to a cross section of 4.8 • 10~ 5 /jb/(AB) which intersects the upper right solid 
curve in Fig. 28 at E Ti ~ 325 GeV. 

Some conclusions can now be drawn. First, since nuclear effects on the PDFs are not large (at most 
10 % at Et > 20 GeV for Pb according to EKS98), an extensive systematic study of the ^-dependence of 
the large-T^r jet cross sections does not seem to be the most urgent need. If, however, the nuclear effects 
in these cross sections turn out to be significant, the measurements of the A systematics, especially 
towards the smallest values of Et, would be necessary in order to understand the new underlying QCD 
collision dynamics. Second, noting the asymmetry in the pseudorapidity distributions in pPb collisions 
at the reduced energy of 5.5 TeV, and to a lesser extent in pAr at 6.3 TeV, when EKS98 corrections are 
implemented, a pA run at reduced proton energies might test the nuclear PDF effects; a study of En- 
distributions in different r/j-slices going from the backward (A) to the forward (p) region would scan 
different regions in x of the ratio of the nuclear gluon distribution over that of a free nucleon, which is 
so far not well constrained. Nevertheless, this would require a detailed comparison with pp results at 
the same energy in order to disentangle possible biases (e.g. detector effects). If these pp results were 
not available experimentally, the uncertainty due to the scale dependence, which in our computations at 
reduced energies is as large as the asymmetry, in the extrapolation from 14 TeV to the pA energy should 
be reduced as much as possible for this study to be useful. 

5.32 High-ET dijet momentum imbalance 

Dijet distributions as a function of (f>, the angle between the two hardest jets, can test the perturbative 
QCD expansion. At LO in collinear factorization [17, 184, 185], the two jets are produced back-to-back, 
so that any deviation from a peak at (f> = it is a signal of NLO corrections. However, recall that the 
results near = 7r are unreliable [193] since the NLO corrections become large and negative, requiring 
an all-order resummation in this region. 

Uncertainties 

In Fig. 29, the same uncertainties examined for ETi and rji are studied in the dijet cross sections as 
a function of 4>. First, the smaller the scale, the larger the results, so that the results for fi = Et are 
up to ~ 50 % smaller than for /i = £r/4. Second, the PDF choice has less than a 10 % effect when 
comparing CTEQ5M with MRST98 central gluon. Third, the variation due to isospin, obtained from the 
comparison of pp and pPb without nuclear modifications of PDFs, is again negligible, and the effect of 
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Fig. 28: Energy dependence of jet cross sections as a function of En for \r]i\ < 2.5 (upper plots) and r\i for En > 20 GeV 
(lower plots). Left-hand side: pp collisions at 14 (solid), 6.3 (dashed) and 5.5 TeV (dotted). Center. pAr collisions at 9.4 (solid) 
and 6.3 TeV (dashed). Right-hand side: pPb collisions at 8.8 (solid) and 5.5 TeV (dashed). Default options are used, see text. 



including EKS98 [27, 28] is also small, ~ 3 %. Finally, the jet-finding algorithm gives an uncertainty of 
order 20 % for a cone with R = 1, R = 0.7 and the fer-clustering algorithm with D = 1. 

Results 

Our results for pp and pA collisions at different energies are presented in Fig. 30. It is clear that the 
dijet momentum imbalance should be measurable provided jet reconstruction is possible in the nuclear 
environment and no other physics contribution, such the underlying event or multiple hard parton scat- 
tering, interferes, spoiling the comparison with experimental data. Only extensive studies using Monte 
Carlo simulations including full event reconstruction will be able to clarify whether such measurements 
are feasible, see the Section on Jet Detection at CMS in the jet Chapter. 

5.4 Direct photons 

For direct photons, we refer to the photon Chapter by P. Aurenche et al in this report. 
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Fig. 29: Uncertainties in dijet cross sections as a function of <f> for Eti > 20 GeV, Et2 > 15 GeV and |»ji|, \r/2\ < 2.5 inpp 
collisions at 14 TeV. Upper left: Scale dependence; Upper right plot: dependence on jet reconstruction algorithm; Lower left: 
PDF dependence in pp collisions at 14 TeV. Lower right plot shows the isospin and nuclear PDF dependence in pp and pPb 
collisions at 5.5 TeV. Unless otherwise stated, default options are used, see text. 
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Fig. 30: Energy dependence of dijet cross sections as a function of 4> for En > 20 GeV, Et2 > 15 GeV and |»ji|, |»j2| < 2.5. 
Upper: pp collisions; Center: pAr collision; Lower: pPb collisions. Default options are used, see text. 



5.5 Quarkonium transverse momentum distributions 

Ramona Vogt 

Since a major experimental focus at the LHC will be quarkonium suppression in heavy ion colli- 
sions, an important plasma probe [222] already at the SPS [223], it is important that systematic studies 
also be undertaken at the LHC. These systematics are especially valuable for the T family where high 
statistics measurements of T production will be available in pp, pA, and AA interactions for the first 
time at the LHC. 

Studies of quarkonium production in pA interactions are particularly relevant for LHC because 
the energy dependence of cold nuclear matter effects such as nucleon absorption are not fully understood 
and need to be checked with measurements. The existence of nuclear effects beyond the modifications 
of the parton densities in the nucleus argues that quarkonium production is not a good reaction by which 
to study the nuclear gluon distribution. However, the nuclear gluon distribution can be obtained through 
careful measurements of heavy flavor production in pp and pA interactions and used as an input to 
quarkonium production. Thus quarkonium is still a benchmark process to test cold matter effects in pA 
interactions at the LHC. It is important to study the shape of the quarkonium px distributions in pA 
interactions because plasma effects on the px distributions are expected to persist to high px [224, 225]. 
Detailed comparisons of the px dependencies of different quarkonium states in AA collisions will reveal 
whether or not the ratios of e.g. ip'/J/ip and T'/T and T"/T are independent of px, as predicted 
by the color evaporation model [226-228] or vary over the px range, as expected from nonrelativistic 
QCD [229]. Indeed quarkonium studies in pA may provide a crucial test of the quarkonium production 
process because it is likely to be possible to separate the P and S states only in pA interactions where 
the accompanying decay photons can be found. 

Two approaches have been successful in describing quarkonium production phenomenologically — 
the color evaporation model (CEM) [226, 227] and nonrelativistic QCD (NRQCD) [229]. Both are briefly 
discussed. 

In the CEM, any quarkonium production cross section is some fraction Fc of all QQ pairs below 
the HH threshold where H is the lowest mass heavy hadron containing Q. Thus the CEM cross section is 
simply the QQ production cross section with a mass cut imposed but without any contraints on the color 
or spin of the final state. The produced QQ pair then neutralizes its color by interaction with the collision- 
induced color field — "color evaporation". The Q and the Q either combine with light quarks to produce 
heavy-flavored hadrons or bind with each other in a quarkonium state. The additional energy needed 
to produce heavy-flavored hadrons is obtained nonperturbatively from the color field in the interaction 
region. The yield of all quarkonium states may be only a small fraction of the total QQ cross section 
below the heavy hadron threshold, 2m u- At leading order, the production cross section of quarkonium 
state C is 

a CEM = F C ^2 H d§ dx ^ f!( x ^Q 2 ) fj( x ^ Q 2 ) M*) s ( s - X ^ S ) , (43) 

ij m Q J 

where ij = qq or gg and aij(s) is the ij — > QQ subprocess cross section. The fraction Fc must be 
universal so that, once it is fixed by data, the quarkonium production ratios should be constant as a 
function of ^/s, xf and px- The actual value of Fc depends on the heavy quark mass, tuq, the scale, 
Q 2 , the parton densities and the order of the calculation. 

Of course the leading order calculation in Eq. (43) is insufficient to describe high px quarkonium 
production since the QQ pair px is zero at LO. Therefore, the CEM was taken to NLO [228, 230] using 
the exclusive QQ hadroproduction code of Ref. [214]. At NLO in the CEM, the process gg — > gQQ 
is included, providing a good description of the quarkonium px distributions at the Tevatron [230]. The 
values of Fc for the individual charmonium and bottomonium states have been calculated from fits 
to the total J/ip and inclusive T data combined with relative cross sections and branching ratios, see 



Refs. [224,231] for details. The pp total cross sections at 5.5 TeV are shown in Tables 13 and 14 in 
Section 6.6, along with the fit parameters for the partem densities, quark masses and scales given for QQ 
production in Table 9 in Section 6.4. 

NRQCD was motivated by the success of potential models in determining the quarkonium mass 
spectra. It was first applied to high px quarkonium production when it was clear that the leading color 
singlet contribution severely underestimated direct J/ip and tp' production at the Tevatron [232,233]. 
NRQCD describes quarkonium production as an expansion in powers of v, the relative Q-Q velocity. 
Thus NRQCD goes beyond the leading color singlet state to include color octet production. Now the 
initial QQ quantum numbers do not have to be the same as in the final quarkonium state because an 
arbitrary number of soft gluons with E < itlqv can be emitted before bound state formation and change 
the color and spin of the QQ. These soft gluon emissions are included in the hadronization process. 

The cross section of quarkonium state C in NRQCD is 

^nrqcd = f 1 dxidx2 ff( Xl ,Q^ff( X2 ,Q^C^ [n] {Q 2 ){0 C n ) , (44) 

i,j n J 

where the partonic cross section is the product of perturbative expansion coefficients, (Q 2 ), and 

nonperturbative parameters describing the hadronization, (C?„ ). The expansion coeffients can be factor- 
ized into the perturbative part and the nonperturbative parameters because the production time, oc I/tuq, 
is well separated from the timescale of the transition from the QQ to the bound state, oc 1/rriQV. If 
1/rriQV 3> 1/rriQ, the bound state formation is insensitive to the creation process. In principle, this 
argument should work even better for bottomonium because mg is larger. The color singlet model result 
is recovered when v — > [234]. The NRQCD parameters were obtained from comparison [235-238] to 
unpolarized high px quarkonium production at the Tevatron [232,233]. The calculations describe this 
data rather well but fail to explain the polarization data [239]. 

Both the CEM and NRQCD can explain the Tevatron quarkonium data because both are dominated 
by gg gg* and gg -► gQQ at high p T . In the CEM g* -► QQ while in NRQCD, g* -► QQ[ 3 s{ 8) ] 
but QQ -> QQ^sj^Zp^] [234]. Direct J/ip and V' production at the Tevatron are predominantly 
by 3 s[ 8 ^ states followed by 1 5q 8 ' ) + 3 -Pj 8 ^ states and, while the \c singlet contribution does not under- 
predict the yield as badly, a 3 s[ 8 ^ component is needed for the calculations to agree with the data. The 
CEM includes the octet components automatically but the relative strengths of the octet contributions 
are all determined perturbatively. The two approaches differ in their predictions of quarkonium polar- 
ization. NRQCD predicts transverse polarization at high px [234] while the CEM predicts unpolarized 
quarkonium production. Neither approach can successfully explain the polarization data [239]. 

The bottomonium results at the Tevatron [240] extend to low px where the comparison to NRQCD 
is somewhat inconclusive. Some intrinsic kx is needed to describe the data below pt ~ 5 GeV. 

In this section, we only present the CEM px distributions. The NRQCD approach has not yet been 
fully calculated beyond LO [234]. We include intrinsic transverse momentum broadening, described 
further in Section 6.7, although it is a rather small effect at high px- The direct J/ij) and T px distribu- 
tions above 5 GeV in pp interactions at 5.5 TeV are shown in Fig. 31 for all the fits to the quarkonium 
data, described in Section 6.6. They all agree rather well with each other. However, the GRV98 result is 
slightly larger than the other J/ip curves at high px- Thus changing the quark mass and scale does not 
have a large effect on the px dependence. In general, the J/ ijj px distributions are more steeply falling 
than the T distributions. 

Figure 32 shows the px distributions in pp and pPb interactions at 5.5 and 8.8 TeV along with pp 
interactions at 14 TeV. The shadowing effect is negligible at high px so that the two distributions are 
indistinguishable. 
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Fig. 3 1 : The pr distributions for J/ip (left) and T (right) production assuming (fcy ) = 1 GeV 2 . On the left hand side, the solid 
curve employs the MRST HO distributions with m c = Q/2 = 1.2 GeV, the dashed, MRST HO with m c = Q = 1.4 GeV, the 
dot-dashed, CTEQ 5M with m c = Q/2 = 1.2 GeV, and the dotted, GRV 98 HO with m c = Q = 1.3 GeV. On the right hand 
side, the solid curve employs the MRST HO distributions with m b = Q = 4.75 GeV, the dashed, m b = Q/2 = 4.5 GeV, the 
dot-dashed, m b = 2Q = 5 GeV, and the dotted, GRV 98 HO with rn b = Q = 4.75 GeV. 
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Fig. 32: ThepT distributions for J/ip (left) and T (right) production assuming (k^) pp = 1 GeV 2 forpp interactions at 5.5 
(lower), 8.8 (center) and 14 (upper) TeV On the left hand side, we use the MRST HO distributions with m c = Q/2 = 1.2 
GeV while on the right hand side, we use the MRST HO distributions with rn b = Q = 4.75 GeV. 



6. PROCESSES WITH POTENTIALLY LARGE NUCLEAR EFFECTS 



6.1 Drell-Yan cross sections at low M 2 : rapidity distributions 

Ramona Vogt 

Low mass Drell-Yan production could provide important information on the quark and antiquark 
distributions at low x but still rather high Q 2 relative to the deep-inelastic scattering measurements. 
When M = 4, 11,20 and 40 GeV respectively, x ~ 0.00073,0.002,0.0036 and 0.0073 at y = and 
\fs = 5.5 TeV. Such low values of x have been measured before in nuclear deep-inelastic scattering but 
always at Q 2 < 1 GeV 2 , out of reach of perturbative analyses. Thus, these measurements could be a 
benchmark for analysis of the nuclear quark and antiquark distributions. 

As previously explained in Sec. 5.1, in AA collisions, the inclusive Drell-Yan contribution is 
unlikely to be extracted from the dilepton continuum which is dominated by correlated and uncorrected 
QQ (DD, BB) pair decays [79, 80]. However, in pA, the low mass Drell-Yan may be competitive with 
the QQ decay contributions to the dilepton spectra [79]. The number of uncorrected QQ pairs is also 
substantially reduced from the AA case [73]. Further calculations of the heavy quark contribution to 
dileptons are needed at higher masses to determine whether the Drell-Yan signal ever dominates the 
continuum or not. Perhaps the Drell-Yan and QQ contributions can be separated using their different 
angular dependencies, isotropic for QQ and oc (1 + cos 2 0) for Drell-Yan, as proposed by NA60 [241] 
and discussed in Sec. 6.2. When M is small and px is large relative to M, the QQ background may be 
reduced and the DY signal observable. This is discussed in the next section. 

The calculation of the Drell-Yan cross section is described in Section 5.1. Since the nuclear effects 
are larger at low Q 2 , we will focus on the variation of the results with scale and shadowing parameter- 
ization. We again make our calculations with the MRST HO [172] partem densities with Q 2 = M 2 
and the EKS98 parameterization [27, 28] in the MS scheme. The Drell-Yan convolutions with nuclear 
parton densities including isospin and shadowing effects are the same as for the Z° in the appendix of 
Ref. [171]. Only the pref actor changes from to Bf- . We present results for pp, Vbp and pPb col- 
lisions at 5.5 TeV/nucleon, the same energy as the Pb+Pb center of mass for better comparison without 
changing x. This is most desirable from the theoretical point of view because the nuclear parton densities 
are most easily extracted from the data. However, as was not touched upon in Section 5.1, the desirability 
of this approach ultimately depends upon the machine parameters. The highest pPb and Pbp luminosities 
can be achieved at 8.8 TeV/nucleon. Running the machine at lower energies will reduce the maximum 
luminosity, up to a factor of 20 if only injection optics are used [242], which could reduce the statistical 
significance of the measurement. This issue remains to be settled. 

Figure 33 shows the rapidity distributions in pPb and Pbp collisions at 5.5 TeV/nucleon in four 
mass bins: 4 < M < 9 GeV, 11 < M < 20 GeV, 20 < M < 40 GeV and 40 < M < 60 GeV. The 
shadowing effect decreases as the mass increases, from a ~ 25% effect at y = in the lowest mass bin 
to a ~ 10% effect in the highest bin. Isospin is a smaller effect on low mass Drell-Yan, as is clear from 
Fig. 34, where the pPb/pp and Vbplpp ratios are given. 

The pVblpp ratio without shadowing is indistinguishable from unity over the whole rapidity range. 
The Vbplpp ratio without shadowing only differs from unity for y < 2. The dashed curves in Fig. 34 
represent the Vbplpp ratios with EKS98 shadowing and shows the pattern of increasing x\ through the 
'antishadowing' region into the EMC region. The ratio does not increase above unity because we present 
the Ap/pp ratio rather than a ratio of Ap with shadowing to Ap without shadowing. This difference can 
be easily seen by comparison to Fig. 33 where the dotted curves are equivalent to the pp distributions. 
The dashed curves in Fig. 34 are the ratios of the dashed to solid curves of Fig. 33. Clearly, at the peak of 
the rapidity distributions, the dashed curves are higher than the solid curves, signifying antishadowing, 
yet do not cross the dotted curves. The isospin of the Pb nucleus thus causes the Vbplpp ratio to remain 
less than unity over all rapidity. Note also that the peaks of the ratios occur at smaller values of x\, 
from x\ ~ 0.11 in the lowest bin to ~ 0.062 in the highest. The evolution of the antishadowing region 
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Fig. 33: The Drell-Yan rapidity distributions in pPb and Pbp collisions at 5.5 TeV evaluated at Q — M for 4 < M < 9 GeV, 
11 < M < 20 GeV, 20 < M < 40 GeV and 40 < M < 60 GeV. The solid and dashed curves show the results without 
and with shadowing respectively in Pbp collisions while the dotted and dot-dashed curves give the results without and with 
shadowing for pPb collisions. 
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Fig. 34: The ratios of pPb and Pbp collisions to pp collisions at 5.5 TeV evaluated at Q — M for 4 < M < 9 GeV, 
11 < M < 20 GeV, 20 < M < 40 GeV and 40 < M < 60 GeV, as a function of rapidity. The solid and dashed curves 
show the Pbp/pp ratios without and with shadowing respectively while the dotted and dot-dashed curves are the pPb/pp ratios 
without and with shadowing. 



broadens and the average x decreases somewhat with increasing Q 2 . On the other hand, the pVblpp ratios 
are rather constant, decreasing only slowly with rapidity although the slope of the ratio increases with 
M. The pPb calculations are always at low X2 where the shadowing parameterization does not change 
much with X2- These results are reflected in the integrated cross sections in Table 8. 

Table 8: Rapidity integrated Drell-Yan cross sections at 5.5 TeV/nucleon in the indicated mass intervals. The calculations are 
to NLO with the MRST HO parton densities and Q 2 = M 2 . 
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In Fig. 35 we show the effects of scale (upper two plots) and shadowing parameterization (lower 
two plots) on the results. 

We choose the lowest mass bin to emphasize the overall effect. The scale variation, discussed in 
Section 5.1, shows that the cross sections are larger for scale 2M and smallest for scale M/2. This scale 
dependence is opposite to typical expectations since lower scales generally lead to larger cross sections. 
At lower energies and other processes, this is true but at high energies, the low x evolution exceeds the 
effect of the scale dependent logarithms [170]. The scale variation is a 50% effect in this bin, decreasing 
to 10% for 40 < M < 60 GeV. The HPC [173] and HKM [29] shadowing calculations have the same 
general trends as the EKS98 parameterization. The HKM model has the weakest effect on sea quarks at 
low x and produces the largest ratios for pPb/pp and in Vbplpp for y < 3. The HPC parameterization has 
the strongest shadowing without evolution and is generally lower than the other two. 

If the low Q 2 Drell-Yan cross section can be extracted from the other contributions to the contin- 
uum in pA, it could provide useful information on the low x, moderate Q 2 regime not covered in previous 
experiments. Comparison could be made to the only previous available pA Drell-Yan measurements as 
a function of xp by the E772 [54] and E866 [55] collaborations at much higher x. 

RHIC pA data should also be available by the time the LHC starts up which could increase the 
^fs systematics, providing, with both colliders, a collective measurement over a wide range of x values 
at similar Q 2 . 
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Fig. 35: The upper two plots present the scale dependence of the Drell-Yan cross section in the mass interval 4 < M < 9 
GeV in pPb and Pbp collisions at 5.5 TeV respectively. The EKS98 parameterization is used in all cases. The solid curve is 
calculated with Q = M, the dashed curve with Q = M/2 and the dot-dashed with Q = 2M. The lower plots present the 
dependence on the shadowing parameterization in the mass range 4 < M < 9 GeV for the ratios pPhlpp and Pbp/pp at 5.5 
TeV. The solid curve is the EKS98 parameterization, the dashed curve is the HPC parameterization and the dot-dashed curve is 
the HKM parameterization. 



6.2 Low mass Drell-Yan production at LHC energies 

George Fai, Jianwei Qiu, and Xiaofei Zhang 

6.21 Introduction 

Massive dilepton production in hadronic collisions is an excellent laboratory for theoretical and exper- 
imental investigations of strong interaction dynamics, and it is a channel for discovery of quarkonium 
states and a clean process for the study of the PDF. In the Drell-Yan process, the massive lepton pair is 
produced via the decay of an intermediate virtual photon 7*. If both the physically measured dilepton 
mass M and the transverse momentum px are large, the cross section in a collision between hadrons 
(or nuclei) A and B, A(Pa) + B(P B ) — > 7*(^ 11) + X, can be factorized systematically in QCD 
perturbation theory and expressed as [17] 

^dM^dpl = (3^) £ f dxifa/A(xi^) I dx 2 f b/B ( X2 ^)^^( Xl , X2 ,M,p T ,y-,p). 

1 a,b 1 

(45) 

The sum J2 a b runs over au> P^ton flavors, f a / A and f b / B axe normal parton distributions and fi represents 
the renormalization and factorization scales. The partonic cross section da^^* x / ' dp\dy in Eq. (45), 
the short-distance probability for partons of flavors a and b to produce a virtual photon of invariant mass 
M, is perturbatively calculable in terms of a power series in a s (//). The scale p, is of the order of the 

energy exchange in the reaction, jx ~ J M 2 + p\. 



When M ~ Mz, high mass dilepton production in heavy ion collisions at LHC energies is dom- 
inated by the Z° channel and is an excellent hard probe of QCD dynamics [243]. In this contribution, 
we demonstrate that the transverse momentum distribution of low mass (Aqcd *C M <C Mz) dilepton 
production at LHC energies is a reliable probe of both hard and semihard physics at LHC energies and 
is an advantageous source of constraints on gluon distributions in the proton and in nuclei [244]. In 
addition, it provides a significant contribution to the total dilepton spectra at the LHC, also an important 
channel for quarkonium and heavy flavour decays. 

The transverse momentum (pr) distribution of the dileptons can be divided into three regions: low 
Pt (-C M), intermediate px (~ M), and high pt (3> M). When both the physically measured M and 
Pt are large and are of the same order, the short-distance partonic part da^^* x / 'dp T dy in Eq. (45) can 
be calculated reliably in conventional fixed-order QCD perturbation theory in terms of a power series 
in a s (p). However, when pt is very different from M, the calculation of Drell-Yan production in both 
low and high px regions becomes a two-scale problem in perturbative QCD and the calculated partonic 
parts include potentially large logarithmic terms proportional to a power of Iu(M/pt). As a result, the 
higher-order corrections in powers of a s are not necessary small. The ratio a NLO /a LO oc a s x (large 
logarithms) can be of order 1 and convergence of the conventional perturbative expansion in powers of 
a s is possibly impaired. 

In the low pt region, there are two powers of ln(M 2 /p T ) for each additional power of a s and the 
Drell-Yan p? distribution calculated in fixed-order QCD perturbation theory is known to be unreliable. 
Only after all-order resummation of the large cc™ In +1 (M 2 /p T ) logarithms do predictions for the pt 
distributions become consistent with data [177, 245]. We demonstrate in Sec. 6.22 that low mass Drell- 
Yan production at pt as low as Aqcd at LHC energies can be calculated reliably in perturbative QCD 
with all order resummation. 

When pt > M/2, the lowest-order virtual photon "Compton" subprocess g + q — ► 7* + q 
dominates the pt distribution, and the high-order contributions including all-order resummation of 
a" ln n ~ 1 (M 2 /p T ) preserve the fact that the pt distributions of low mass Drell-Yan pairs are domi- 
nated by gluon initiated partonic subprocesses [246]. We show in Sec. 6.23 that the pt distribution of 
low mass Drell-Yan pairs can be a good probe of the gluon distribution and its nuclear dependence. We 
give our conclusions in Sec. 6.24. 

6.22 Low mass Drell-Yan production at low transverse momentum 

Resummation of large logarithmic terms at low pt can be carried out in either pt or impact parameter 
(b) space, the Fourier conjugate of pr space. All else being equal, the b space approach has the advantage 
of explicit transverse momentum conservation. Using renormalization group techniques, Collins, Soper, 
and Sterman (CSS) [178] devised a b space resummation formalism that resums all logarithmic terms as 
singular as (l/p T ) \n m (M 2 /p T ) when px — ► 0. This formalism has been widely used for computations 
of vector boson pt distributions in hadron reactions [247]. 

At low mass M and pr, Drell-Yan transverse momentum distributions calculated in the CSS b- 
space resummation formalism strongly depend on the nonperturbative parameters determined at fixed 
target energies. However, it was pointed out recently that the predictive power of pQCD resummation 
improves with y/s and, when the energy is high enough, pQCD should have good predictive power even 
for low mass Drell-Yan production [177]. The LHC will provide us a chance to study low mass Drell-Yan 
production at unprecedented energies. 

In the CSS resummation formalism, the differential cross section for Drell-Yan production in 
Eq. (45) is reorganized as the sum 



d° 'ab->u(m)x 
dM 2 dy dp T 



7 (resum) , (Y) 

AB—tll(M)X aC7 AB^ll(M)X 

dM 2 dy dp\ + dM 2 dydp\ ' 



(46) 



The all-orders resummed term is a Fourier transform from the 6-space, 



, (resum) 

dM^dpl = W?J d2 ~ h ^ T:b W(b,M, XAl x B ) = -j dbJ (p T ~b)~bW(b,M,x A ,X B ) , 

(47) 

where Jo is a Bessel function, x A = e y M/yfs and xb = e~ y M/ yfs, with rapidity y and collision 
energy sfs. In Eq. (46), the cr ( rcsum ) term dominates the pr distributions when px <C M, while the 
term gives negligible corrections for small px but becomes important when pr ~ M. 

The function W(b, M, x A ,x B ) resums to all orders in QCD perturbation theory the singular terms 
that would otherwise behave as 5 2 (pt) and (1/p^) ln m (M 2 /p|.) in transverse momentum space for all 
m > 0. It can be calculated perturbatively for small b, 

W(b, M, x A , x B ) = e- s ^ W(b, c/b, x A , x B ) = W pcrt (b, M, x A , x B ) ■ (48) 

All large logarithms from ln(c 2 /6 2 ) to ln(M 2 ) have been completely resummed into the exponential 
factor S(b, M) = fj?l 2 d^ 2 /^ [ln(M 2 /^, 2 )„4(a s (^)) + B(a s (n))] where the functions A and B are 
given in Ref. [178] and c = 2e~ lE with Euler's constant 7^ « 0.577. With only one large momentum 
scale 1/6, the function W(b, c/b, x A ,x B ) in Eq. (48) is perturbatively calculable and is factorized as 



W(b, c/b,x A ,x B ) = ^2 

a,b 



/a/A(6,V6)®C (^ 



h/B{&M~b)®c b ^f) 

?2 . 



do, (49) 



where the sum J2 a b runs over au parton flavors, the <g) represents the convolution over parton momen- 
tum fraction, and cr is the lowest order partonic cross section for the Drell-Yan process [178]. The 
functions C = ^2 n C^(a s /n) n in Eq. (49) are perturbatively calculable and given in Ref. [178]. Since 
the perturbatively resummed W pert (b, M,x A ,x B ) in Eq. (48) is only reliable for the small b region, an 
extrapolation to the nonperturbative large b region is necessary in order to complete the Fourier transform 
in Eq. (47). 

In the original CSS formalism, a variable 6* and a nonperturbative function FQg S (b, x A ,x B ) 
were introduced to extrapolate the perturbatively calculated J^ pcrt into the large b region. The full b- 
space distribution was of the form 

W css (b,M,x A ,x B ) = W pcvt (k,M,x A ,XB)F^i s (b,M, XA ,x B ), (50) 



where 6* = b/\Jl + (fe/6 max ) 2 , with 6 max = 0.5 GeV -1 . This construction ensures that 6* < fe max for 
all values of b. 

In terms of the 6* formalism, a number of functional forms for the Fq$ s have been proposed. A 
simple Gaussian form in b was first proposed by Davies, Webber and Stirling (DWS) [248], 

F^ s (b,M,x A ,x B ) =exp{-( 5l + 52 ln(M/2M ))6 2 } , (51) 

with the parameters Mq = 2 GeV, g\ = 0.15 GeV 2 , and 52 = 0.4 GeV 2 . In order to take into account the 
apparent dependence on collision energies, Ladinsky and Yuan (LY) introduced a new functional form 
[249], 

F^(b,M,x A ,x B ) = exp{-( gi +g 2 ln(M/2M ))b 2 -g ig3 \n(W0x AX B)b} , (52) 

with M = 1.6 GeV, gi = 0.11+°;^ GeV 2 , g 2 = 0.58^^ GeV 2 , and g 3 = -1.5+°;} GeV" 1 . Recently, 
Landry, Brook, Nadolsky, and Yuan proposed a modified Gaussian form [250], 

Fblny(&, M, x a ,x b ) = exp {- [ 9l + g 2 ln(M/2M ) + gi 93 HW0x A x B )} 6 2 } , (53) 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

b(1/GeV) 



Fig. 36: The &-space resummed functions bW(b) in Eq. (47) 
for Drell-Yan production of dilepton mass M = 5 GeV at 
= 5.5 TeV with the QZ (solid), DWS (dashed), LY (dot- 
ted), and BNLY (dot-dashed) formalism of nonperturbative 
extrapolation. 
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Fig. 37: Differential cross section da/dM dydp T for pro- 
duction of Drell-Yan pairs of M = 5 GeV in pp collisions 
at the LHC with y = 0, and yfs = 5.5 TeV with the NNLL 
(solid), NLL (dashed), and LL (dotted) accuracy. 



with M = 1.6 GeV, g x = 0.21+°;^ GeV 2 , g 2 = 0.68+°;^ GeV 2 , and g 3 = -0.6+°;°^. All these param- 
eters were obtained by fitting low energy Drell-Yan and high energy W and Z data. Note, however that 
the 6* formalism introduces a modification to the perturbative calculation. The size of the modifications 
strongly depends on the nonperturbative parameters in i ?NP (6, M, xa,xb), M, and yfs [247]. 

A remarkable feature of the 6-space resummation formalism is that the resummed exponential 
factor exp[— S(b, M)] suppresses the 6-integral when b is larger than \ jM. It can be shown using the 
saddle point method that, for a large enough M, QCD perturbation theory is valid even at px = 
[178]. For high energy heavy boson (W, Z, and Higgs) production, the integrand of 6-integration in 
Eq. (47) at pr = is proportional to bW(b, Q, xa,xb), with a saddle point b sp well within the pertur- 
bative region (b sp < 6 max )- Therefore, the 6-integration in Eq. (47) is dominated by the perturbatively 
resummed calculation. The uncertainties from the large- b region have very little effect on the calculated 
Pt distributions and the resummation formalism has good predictive power. 

On the other hand, in low energy Drell-Yan production, there is no saddle point in the perturbative 
region for the integrand in Eq. (47). Therefore the dependence of the final result on the nonperturbative 
input is strong [177, 245]. However, as discussed in Refs. [177, 179, 247], the value of the saddle point 
strongly depends on ^fs in addition to its well-known M 2 dependence. 

Figure 36 shows the integrand of the 6-integration in Eq. (47) at pt = for production of Drell- 
Yan pairs of mass M = 5 GeV in pp collisions at yfs = 5.5 TeV and 6 max = 0.5 GeV -1 . The curves 
represent different extrapolations to the large- b region. The three curves (dashed, dotted, and dot-dashed) 
are evaluated using the b* formalism with the DWS, LY, and BLNY nonperturbative functions, respec- 
tively. Although these three nonperturbative functions give similar 5-space distributions for heavy boson 
production at Tevatron energies, they predict very different 6-space distributions for low mass Drell- 
Yan production at LHC energies even within the perturbative small-6 region. Since the 6-distribution 
in Fig. 36 completely determines the resummed pt distribution through the ^-integration weighted by 
the Bessel function Joiprb), we need to worry about the uncertainties of the resummed low mass pt 
distributions calculated with different nonperturbative functions. 



In order to improve the situation, a new formalism of extrapolation (QZ) was proposed [177], 

Wfh M ,f W^(b,M,x A ,x B ) b<b max 

w { o,M,x A ,x B) I W pcrt (b max , M, xa,xb) F NP (b, M, xa, %b', b max , a) b > 6 max ' ^ V 

where the nonperturbative function F NP is given by 

F np = exp |_ i n (M 2 6Lx/c 2 ) [91 ((b 2 ) a - (6Lx) a ) +92 (b 2 - b 2 max )] - g 2 (b 2 - b 2 max ) } . (55) 

Here, 6 max is a parameter to separate the perturbatively calculated part from the nonperturbative input, 
and its role is similar to the 6 max in the 6* formalism. The term proportional to g\ in Eq. (55) represents 
a direct extrapolation of the resummed leading power contribution to the large b region. The parameters 
gi and a are determined by the continuity of the function W(b\ M, xa,xb) at 6 max . On the other hand, 
the values of g 2 and g 2 represent the size of nonperturbative power corrections. Therefore, sensitivity to 
g 2 and g~2 in this formalism clearly indicates the precision of the calculated px distributions. 

The solid line in Fig. 36 is the result of the QZ parameterization with 6 max = 0.5 GeV^ 1 and 
92 = 92 = 0. Unlike in the 6* formalism, the solid line represents the full perturbative calculation and 
is independent of the nonperturbative parameters for b < 6 max . The difference between the solid line 
and the other curves in the small b region, which can be as large as 40%, indicates the uncertainties 
introduced by the 6* formalism. 

It is clear from the solid line in Fig. 36 that there is a saddle point in the perturbative region even 
for dilepton masses as low as M = 5 GeV in Drell-Yan production at yfs = 5.5 TeV. At that energy, 
xa, xb ~ 0.0045. For such small values of x, the PDFs have very strong scaling violations, leading to 
large parton showers. It is the large partem shower at the small x that strongly suppresses the function 
W(b, c/b, xa,xb) in Eq. (48) as b increases. Therefore, for 

frmax ~(a few GeV) \ the predictive power 
of the 6-space resummation formalism depends on the relative size of contributions from the small-6 
(b < 6 max ) and large-& (b > S max ) regions of the 6-integration in Eq. (47). With a narrow b distribution 
peaked within the perturbative region for the integrand, the ^-integration in Eq. (48) is dominated by the 
small-6 region and, therefore, we expect pQCD to have good predictive power even for low M Drell-Yan 
production at LHC energies. 

Figure 37 presents our prediction of the fully differential cross section da / dM 2 dydp 2 ^ for Drell- 
Yan production in pp collisions at =5.5 TeV and y = [244]. Three curves represent the different 
order of contributions in a s to the perturbatively calculated functions A(a s ), B(a s ), and C(a s ) in the 
resummation formalism [251]. The solid line represents a next-to-next-to-leading-logarithmic (NNLL) 
accuracy corresponding to keeping the functions, A(a s ), B(a s ), and C(a s ) to order af, a 2 , and a], 
respectively. The dashed line has next-to-leading-logarithmic (NLL) accuracy with the functions, A(a s ), 
B(a s ), and C{a s ) at a 2 , a s , and respectively, while the dotted line has the lowest leading-logarithmic 
(LL) accuracy with the functions, A(a s ), B(a s ), and C{a s ) at the a s , a° s , and a° s , respectively. Similar to 
what was seen in the fixed order calculation, the resummed p? distribution has a if-factor about 1.4-1.6 
around the peak due to the inclusion of the coefficient 

In Eq. (55), in addition to the g\ term from the leading power contribution of soft gluon showers, 
the #2 term corresponds to the first power correction from soft gluon showers and the g~2 term is from 
the intrinsic transverse momentun of the incident parton. The numerical values of g 2 and g 2 have to be 
obtained from fits to the data. From fits of low energy Drell-Yan data and heavy gauge boson data at 
the Tevatron, we found that the intrinsic transverse momentum term dominates the power corrections. It 
has a weak energy dependence. For convenience, we combine the parameters of the b 2 term as G 2 = 
g 2 ln(M 2 6^ ax /c 2 ) + g 2 . For M = 5 GeV and y = 0, we use G 2 ~ 0.25 in the discussion here [177]. 
To test the G 2 dependence of our calculation, we define 

„ , s _ d(T AB^ll(M)X^ I da AB ^U(M)x(PT) 

UgAPt)- dM 2 dydp 2 / dM^dydp 2 ' { ' 
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Fig. 38: The ratio Rg 2 defined in Eq. (56) with G2 = 0.25 Fig. 39: Rs h as a function of p T using EKS shadowing in 

GeV 2 for the differential cross section shown in Fig. 37. pPb and PbPb collisions at = 5.5 TeV. 



where the numerator represents the result with finite G2, and the denominator contains no power correc- 
tions, G2 = 0. 

The result for Rg 2 is shown in Fig. 38. It deviates from unity by less than 1%. The dependence 
of our result on the nonperturbative input is indeed very weak. 

Since the G2 terms represent the power corrections from soft gluon showers and partem intrinsic 
transverse momenta, the smallness of the deviation of Rg 2 from unity also means that leading power 
contributions from gluon showers dominate the dynamics of low mass Drell-Yan production at LHC 
energies. Even though the power corrections will be enhanced in nuclear collisions, we expect them to 
be still less than several percent [179]. The isospin effects are also small here because x A and xb are 
small. 

Since the leading power contributions from initial-state parton showers dominate the production 
dynamics, the important nuclear effect is the modification of the parton distributions. Because xa and xb 
are small for low mass Drell-Yan production at LHC energies, shadowing is the only dominant nuclear 
effect. In order to study the shadowing effects, we define [179] 

-+ll(M)X (Pt) 

sh[PT) = dM^dydp\ / dM*dydp\ ' { ° 

We plot the ratio R s h as a function of px in PbPb collisions at ^fs = 5.5 TeV for M = 5 GeV and y = 
in Fig. 39. The EKS parameterization of nuclear parton distributions [27,28] was used to evaluate the 
cross sections in Eq. (57). It is clear from Fig. 39 that low mass Drell-Yan production can be a good 
probe of nuclear shadowing [244], both in PbPb and in pPb collisions at the LHC. 



6.23 Low mass Drell-Yan production at high transverse momentum 

The gluon distribution plays a central role in calculating many important signatures at hadron colliders 
because of the dominance of gluon-initiated subprocesses. A precise knowledge of the gluon distribution 
as well as its nuclear dependence is absolutely vital for understanding both hard and semihard probes at 
LHC energies. 

It was pointed out recently that the transverse momentum distribution of massive lepton pairs 
produced in hadronic collisions is an advantageous source of constraints on the gluon distribution [252], 
free from the experimental and theoretical complications of photon isolation that beset studies of prompt 



photon production [253,254]. Other than the difference between a virtual and a real photon, the Drell- 
Yan process and prompt photon production share the same partonic subprocesses. Similar to prompt 
photon production, the lowest-order virtual photon "Compton" subprocess g + q — > 7* + q dominates 
the pt distribution when px > M/2. The next-to-leading order contributions preserve the fact that the 
Pt distributions are dominated by gluon-initiated partonic subprocesses [252]. 

There is a phase space penalty associated with the finite mass of the virtual photon, and the Drell- 
Yan factor a cm /(37rM 2 ) < 10~ 3 /M 2 in Eq. (45) renders the production rates for massive lepton pairs 
small at large values of M and px- In order to enhance the Drell-Yan cross section while keeping 
the dominance of the gluon-initiated subprocesses, it is useful to study lepton pairs with low invariant 
mass and relatively large transverse momenta [246]. With the large px setting the hard scale of the 
collision, the invariant mass of the virtual photon M can be small, as long as the process can be identified 
experimentally and M » Aqcd- For example, the cross section for Drell-Yan production was measured 
by the CERN UA1 Collaboration [255] for virtual photon mass M G [2m M , 2.5] GeV. 

When p\ » M 2 , the perturbatively calculated short-distance partonic parts, da a b->-y*x /dp T dy 
in Eq. (45), receive one power of the logarithm ln(p T /M 2 ) at every order of a s beyond leading order. 
At sufficiently large pt, the coefficients of the perturbative expansion in a s will have large logarithmic 
terms, so that these higher-order corrections may not be small. In order to derive reliable QCD predic- 
tions, resummation of the logarithmic terms ln m (p 2 n /M 2 ) must be considered. It was recently shown 
[246] that the large h\ m {p\/M 2 ) terms in the low mass Drell-Yan cross sections can be systematically 
resummed into a set of perturbatively calculable virtual photon fragmentation functions [256]. Similar to 
Eq. (46), the differential cross section for low mass Drell-Yan production at large pr can be reorganized 
as 

ao 'AB^fl(M)X _ UU AB-*ll(M)X AB^ll(M)X 

dM 2 dydp 2 T ~ dM 2 dydp\ dM 2 dydp 2 T ' ( } 

where o-( rcsum ) includes the large logarithms and can be factorized as [246] 

(rcsum) 

f f 

Yl / dx lfa/A(xi,fJ-) / dx 2 f b/B (x 2 ,fi) 

a,b,c 



<ln \li >ll(M)X ( n, 



dM 2 dydp\ \3irM 2 J 

X / ^ ^§ {PCT =PT/Z) D ^x(z,» 2 F ;M 2 ) , (59) 

with the factorization scale [i and fragmentation scale [ip, and virtual photon fragmentation functions 
-D c ^ 7 * (z, fj, F ; Q 2 ). The cr( Dir ) term plays the same role as term in Eq. (46), dominating the cross 
section when px — ► M. 

Figure 40 presents the fully resummed transverse momentum spectra of low mass Drell-Yan pro- 
duction in pp collisions with M = 5 GeV at y = and y/s = 5.5 TeV (solid). For comparison, we 
also plot the leading order spectra calculated in conventional fixed order pQCD. The fully resummed 
distribution is larger in the large p? region and smoothly convergent as pr — ► 0. In addition, as dis- 
cussed in Ref. [246], the resummed differential cross section is much less sensitive to the changes of 
renormalization, factorization, and fragmentation scales and should thus be more reliable than the fixed 
order calculations. 

To demonstrate the relative size of gluon initiated contributions, we define the ratio 

R _ doAB-^f (M)x (gluon-initiated) / da AB ^ T{M)x 
9 dp\ dy I dp\ dy 



The numerator includes the contributions from all partonic subprocesses with at least one initial-state 
gluon while the denominator includes all subprocesses. 
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Fig. 40: Differential cross section do j dM 1 dydp% for pro- Fig. 41: Ratio of gluonic over total contributions to Drell- 

duction of Drell-Yan pairs of M = 5 GeV in pp collisions at Yan production at the LHC, R g , defined in Eq. (60) with 

yfs = 5.5 TeV with low and high pr resummation (solid), M — 5 GeV at ^J~s = 5.5 TeV 
in comparison to conventional lowest result (dashed). 

In Fig. 41, we show R g as a function of px in pp collisions at y = and y/s = 5.5 TeV with 
M = 5 GeV. It is clear from Fig. 41 that gluon-initiated subprocesses dominate the low mass Drell-Yan 
cross section and that low mass Drell-Yan lepton-pair production at large transverse momentum is an 
excellent source of information about the gluon distribution [246]. The slow falloff of R g at large px is 
related to the reduction of phase space and the fact that cross sections are evaluated at larger values of x. 



6.24 Conclusions 

In summary, we present the fully differential cross section of low mass Drell-Yan production calculated 
in QCD perturbation theory with all-order resummation. For px <C M, we use CSS 6-space resumma- 
tion formalism to resum large logarithmic contributions as singular as (l/p T ) ln m (M 2 /p T ) to all orders 
in a s . We show that the resummed px distribution of low mass Drell-Yan pairs at LHC energies is domi- 
nated by the perturbatively calculable small 6-region and thus reliable for px as small as Aqcd- Because 
of the dominance of small x PDFs, the low mass Drell-Yan cross section is a good probe of the nuclear 
dependence of partem distributions. For px S> M, we use a newly derived QCD factorization formalism 
[246] to resum all orders of \n m {p T /M 2 ) type logarithms. We show that almost 90% of the low mass 
Drell-Yan cross sections at LHC energies is from gluon initiated partonic subprocesses. Therefore, the 
low mass Drell-Yan cross section at px > M is an advantageous source of information on the gluon 
distribution and its nuclear dependence — shadowing. Unlike other probes of gluon distributions, low 
mass Drell-Yan does not have the problem of isolation cuts associated with direct photon production at 
collider energies and does not have the hadronization uncertainties of J/ tp and charm production. More- 
over, the precise information on dilepton production in the Drell-Yan channel is critical for studying 
charm production at LHC energies. 



6.3 Angular distribution of Drell-Yan pairs in pA at the LHC 

Rainer J. Fries 



6.31 Introduction 

The Drell-Yan cross section can be described by a contraction of the lepton and hadron tensors 

( '° ~ rl. in ,\V". (61) 



dM 2 dp 2 T dydVl 64vr 3 sM 4 ^ 

Here we parametrize the lepton pair by the invariant mass M of the virtual photon and its transverse 
momentum pt and rapidity y in the center of mass frame of the colliding hadrons. In addition we give 
the direction of one of the leptons, say the positively charged one, in the photon rest frame using polar 
and azimuthal angles <p and 9: dS7 = d<j)d cos 9. Now the cross section can be understood in terms of 
four helicity amplitudes [257-259] 



dM 2 dpi dy dQ 64vr 3 sM 2 



Wtl(1 + cos 2 9) + W L (l/2 - 3/2 cos 2 9) 

(62) 

+Wa sin 29 cos <fi + Waa sin 2 9 cos 2<f> 



These are defined as contractions W a y = e^(a)W^ u e* u (a') of the hadron tensor with polarization 
vectors of the virtual photon for polarizations a = 0, ±1. Only four out of all possible contractions are 
independent, the others can be related by symmetries of the hadron tensor. The usual choice is to pick 
the longitudinal Wl = Wo,o, the helicity flip W A = (W^o + W 7 o,i)/v / 2 and the double helicity-flip 
amplitude W A A = Wi _i together with the trace W T l = W T + W L /2 = -Wtf/2 as a basis. Note that 
integration over the angles 9 and 4> leaves only contributions from the trace 

dM 2 dp\dy = QaJ7m 2 ~ W ^ = 2aJ7m 2 ^~ 9iiv)W ^ ' (63) 

On the other hand, if we are only interested in relative angular distributions — i.e. in the ratio 

16-7T da j da 

~ dM 2 dp\ dy dn I dM 2 dp\ dy 

— we can make use of angular coefficients. Two different sets can be found in the literature [260]. One 
set consists of the coefficients 

^."k, Al = i^, .4 2 = H!^, (65) 

Wtl Wtl " Wtl 



the other one is defined by 



2-3Ao 2A 1 2A 2 

^ = : — ! M= v= ir- (66) 

2 + 2 + A 2 + A 



The helicity amplitudes are frame dependent. In principle we allow all frames where the photon 
is at rest, i.e. qf 1 = (M, 0,0, 0). However there are some frames with particular properties studied in 
[258]. Here we only use the Collins-Soper (CS) frame. It is characterized by two properties. First, the 
y-axis is perpendicular to the plane spanned by the two hadron momenta Pi and P2 (which are no longer 
collinear in the photon rest frame as long as pr / 0, as in our kinematic domain) and second, the z-axis 
cuts the angle between Pi and — P2 into two equal halves, see Fig. 42. 
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Fig. 42: The Collins-Soper frame: the z-axis cuts the angle between Pi and P 2 into halves (the half angle is called the 
Collins-Soper angle 70s) while the x-axis is perpendicular to Pi and P2. The direction of one lepton momentum h can then 
be given by the angles 9 and <j>. 



6.32 Leading twist 

The hadron tensor for the Drell-Yan process in a leading-twist (twist-2) calculation is given by the well 
known factorization formula as a convolution of two partem distributions with a perturbative parton cross 
section. We are interested here in the kinematic region characterized by intermediate photon mass M of 
a few GeV and intermediate transverse momentum pt ~ M. For these values, the dominant contribution 
to twist-2 is given by the NLO perturbative diagrams like the one in Fig. 43 (left) and we can safely omit 
logarithmic corrections of type ln 2 (M 2 /p|,). It has been shown that a leading-twist calculation up to 
NLO respects the so-called Lam-Tung sum rule, Wl = 2Waa [261]. In terms of angular coefficients 
this can be rewritten as Aq = A2 or 2u = 1 — A. Furthermore, the spin-flip amplitude Wa has to vanish 
for a symmetric colliding system like pp. For pA we expect small contributions for Wa due to lost 
isospin symmetry and nuclear corrections to the parton distributions. Results for pp at yfs = 5.5 TeV 
have already been presented elsewhere [262]. 

6.33 Nuclear-enhanced twist-4 

For large nuclei, corrections to the leading-twist calculation, induced by multiple scattering, play an im- 
portant role. The formalism for including these nuclear-enhanced higher-twist contributions was worked 
out by Luo, Qiu and Sterman [36, 263, 264]. The leading nuclear corrections (twist-4 or double scatter- 
ing) have already been calculated for some observables. For Drell-Yan this was first done by Guo [37] 
and later generalized to the Drell-Yan angular distribution [260, 265]. 

Figure 43 (right) shows an example for a diagram contributing at twist-4 level. Now two partons 
a, b from the nucleus and one (c) from the single hadron are involved. As long as pt ~ M, twist-4 is 
dominated by two different contributions. The double-hard (DH) process where each parton from the 
nucleus has a finite momentum fraction and the soft-hard (SH) process where one parton has vanishing 



Fig. 43: Examples for diagrams contributing to the dilepton production in hadron (h) nucleus (A) scattering at twist-2 (left) and 
twist-4 (right) level. 
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Fig. 44: Rapidity dependence of iV CT Wtl (left) and N a Wa (right) for pPb at sfs = 8.8 TeV, M = 5 GeV and p T = 4 GeV: 
twist-2 NLO (short dashed), double hard with EKS98 modifications (long dashed), soft hard (dot-dashed, scaled up by a factor 
of 10) and the sum of all contributions (bold solid line). The double hard contribution calculated without EKS98 modifications 
is also shown (dotted line). Note that the incoming nucleus has positive rapidity. 

momentum fraction. The factorization formulas are given by 



for double-hard and soft-hard scattering respectively [37,260,265]. The function D XaiXb (q,x c ) is a 
second order differential operator in x a and x b while and are new twist-4 matrix elements 
which encode nonperturbative correlations between the partons a and b. Since we are still missing solid 
experimental information about these new quantities, they are usually modeled in a simple way through 



parton distributions. We use T^ b n (x a ,x b ) = CA 4 / 3 f a (x a )f b (x b ) and T^ b u (x a ) = X 2 A 4 / 3 f a (x a ) where 



C and A 2 are normalization constants. The key feature, the nuclear enhancement, is their scaling with 
the nuclear size. 

It has been shown in Refs. [260, 265] that the DH process shows a trivial angular pattern in the 
sense that it is similar to the LO annihilation of on-shell quarks. The only difference is that one of the 
quarks now carries finite pr- In this spirit, it is no surprise that the DH contribution respects the Lam- 
Tung relation. On the other hand, SH scattering is more complicated and violates the Lam-Tung sum 
rule. We also expect that the spin-flip amplitude, Wa, can receive large contributions from the twist-4 
calculation. 

6.34 Numerical results 

In this section, we present some numerical results obtained for pPb collisions at the LHC energy i/i = 
8.8 TeV. Results for RHIC energies can be found elsewhere [260]. We use the CTEQ5L parton distribu- 
tions [86] combined with EKS98 [27] nuclear modifications both for the nuclear parton distributions and 
for the models of the twist-4 matrix elements. In some plots we also give results for double-hard con- 
tributions without nuclear modification. Since we do not know anything about the correct x-dependence 
of the higher-twist matrix elements, this gives an impression of the theoretical uncertainty we may as- 
sume for the higher-twist calculation. The normalization constants for the twist-4 matrix elements are 
chosen to be A 2 = 0.01 GeV 2 and C = 0.005 GeV 2 . In order to enable convenient comparison with 
cross sections, we show all helicity amplitudes multiplied by the prefactor N a = Q 2 m /(64-7r 3 sM 2 ) of 




c 



Tab ( x a,Xb) H^ b+C (q, X a , X b , X c ) f c {x c ) 



D x a ,x b (q, x c )T ab (x a ) H% b+C (q, x a , x b = 0, x c ) f c {x c ) 



(67) 



(68) 



Eq. (62). 
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Fig. 45: The angular coefficients A (long dashed), \i (short dashed) and v (dot dashed) for pPb at ^fs = 8.8 TeV and M = 5 
GeV as functions of rapidity (left, at pr =4 GeV) and transverse momentum (right, at y — 0). Thin lines represent pure 
twist-2 calculations, thick lines show the results including twist-4. The solid line gives the violation of the Lam-Tung relation 
2v — (1 — A) scaled up by a factor of 10. 



In Fig. 44 we give results for the helicity amplitudes Wtl and Wa as functions of rapidity. We 
observe that DH scattering gives a large contribution at negative rapidities (the direction of the proton) 
which can easily balance the suppression of the twist-2 contribution by shadowing. Soft-hard scattering 
is strongly suppressed at these energies. Wl (not shown) has qualitatively the same behavior as Wtl- 
The helicity-flip amplitude Wa picks up only a small contribution from twist-2 and is entirely dominated 
by double-hard scattering, as already expected. This would be a good observable to pin down nuclear 
effects in pA collisions. Figure 45 gives the full set of angular coefficients A, \x and v as functions 
of y and pr- Both the twist-2 results and the modifications by twist-4 are given. The violation of the 
Lam-Tung sum rule is numerically almost negligible since the soft-hard contribution is so small. Note 
here that earlier experiments have discovered a large violation of the Lam-Tung relation in it A collisions 
[267]. This issue is still not fully resolved. 

Thus pA collisions offer the unique opportunity to study nuclear effects directly via the rapidity 
dependence of observables. Helicity amplitudes and angular coefficients can help to pin down the role 
of nuclear-enhanced higher twist. The helicity-flip amplitude Wa and the coefficient /i, which both 
vanish for pp, are particularly promising. Several different species of nuclei are advisable to address the 
important question of A scaling. 

In addition, the angular distribution of Drell-Yan lepton pairs might play a crucial role in the 
extraction of the Drell-Yan process from dilepton data. Dilepton spectra at the LHC have a large contri- 
bution from the decay of correlated and uncorrected heavy quark pairs [79, 80]. However the angular 
pattern of these dileptons is isotropic rather than modulated as in the Drell-Yan case. This provides 
additional means to discriminate between Drell-Yan dileptons and those from heavy quark decays. 



6.4 Heavy flavor total cross sections 

Ramona Vogt 

Heavy quark production will be an important part of the LHC program, with both open heavy 
flavors and quarkonium studied extensively. Heavy quark production will probe the small x region of the 
gluon distribution in the proton and in the nucleus: at y = 0, x ~ 4 x 10 -4 for charm and x ~ 0.0011 for 
bottom. Direct information on the gluon distribution in the proton is still lacking, let alone the nucleus. 
Open heavy flavor production is a better direct test of the nuclear gluon distribution than quarkonium 
which is subject to other cold matter effects such as nuclear absorption. 



Heavy quark studies in pA interactions provide important benchmarks for further measurements 
in AA. More than one QQ pair (bottom as well as charm) will be produced in each central AA collision. 
Therefore heavy quark decays will dominate the lepton pair continuum between the J/ ijj and Z° peaks 
[169]. Although heavy quark energy loss in the medium is currently expected to be small [211], the 
shape of the dilepton continuum could be significantly altered due to modifications of the heavy quark 
distributions by energy loss [80,212]. An understanding of heavy quark production in pA is also a 
necessity for the study of quarkonium production in AA collisions. Models predicting a quarkonium 
enhancement relative to heavy flavor production in nuclear collisions require a reliable measure of the 
heavy flavor production rates in AA collisions with shadowing included to either verify or disprove these 
predictions. 

In this section, we first discuss our estimates of the QQ total cross sections and some of their 
associated uncertainties. We then describe a more detailed measurement of the nuclear gluon distribution 
using dileptons from heavy quark decays. 

At LO heavy quarks are produced by gg fusion and qq annihilation while at NLO qg and qg scatter- 
ing is also included. To any order, the partonic cross section may be expressed in terms of dimensionless 

(k I) 

scaling functions f>- ' that depend only on the variable 77 [213], 

Vij(s,m 2 Q ,Q 2 ) = (4WQ 2 )) fc £ f%%) In' , (69) 

III \ I I I 

fc=o i=o \ Q/ 

where s is the partonic center of mass energy squared, tuq is the heavy quark mass, Q is the scale 
and rj = s/Arriq — 1. The cross section is calculated as an expansion in powers of a s with k = 
corresponding to the Born cross section at order 0(ar s ). The first correction, k = 1, corresponds to the 
NLO cross section at 0(a z s ). It is only at this order and above that the dependence on renormalization 
scale, Qr, enters the calculation since when k = 1 and / = 1 the logarithm ln(Q 2 /nig) appears, as 
in gauge boson and Drell-Yan production at 0(a s ), described in Section 5.1. The dependence on the 
factorization scale, Qp, in the argument of a s , appears already at LO. We assume that Qr = Qf = 
Q. The next-to-next-to-leading order (NNLO) corrections to next-to-next-to-leading logarithm (NNLL) 
have been calculated near threshold [213] but the complete calculation only exists to NLO. 

The total hadronic cross section is obtained by convoluting the total partonic cross section with 
the partem distribution functions (PDFs) of the initial hadrons, 

a pp (s,m 2 Q ) = V L m 2 —5{x lX2 - r) ^ ^llR llA I d i: j{T,m 2 Q ,Q 2 ) , (70) 

* — / Q T X\ X2 

i,j=q,q,g 

where the sum i is over all massless partons and x\ and x 2 are fractional momenta. The partem densities 
are evaluated at scale Q. All our calculations are fully NLO, applying NLO partem distribution functions 
and the two-loop a s to both the 0(a 2 ) and 0(a%) contributions, as is typically done [213,214]. 

To obtain the pp cross sections at LHC, we compare the NLO cross sections to the available cc and 
bb production data by varying the mass, rag, and scale, Q, to obtain the 'best' agreement with the data 
for several combinations of tuq, Q, and PDF. We use the recent MRST HO central gluon [172], CTEQ 
5M [86], and GRV 98 HO [45] distributions. The charm mass is varied between 1.2 and 1.8 GeV with 
scales between m c and 2m c . The bottom mass is varied over 4.25 to 5 GeV with 7X1^/2 < Q < 2m{,. The 
best agreement for Q = m c is with m c = 1.4 GeV and m c = 1.2 GeV is the best choice using Q = 2m c 
for the MRST HO and CTEQ 5M distributions. The best agreement with GRV 98 HO is Q = m c = 1.3 
GeV while the results with Q = 2m c lie below the data for all m c . All five results agree very well with 
each other for pp — > cc. We find reasonable agreement with all three PDFs for = Q = 4.75 GeV, 
m b = Q/2 = 4.5 GeV and m b = 2Q = 5 GeV. 

Before calculating the QQ cross sections at nuclear colliders, some comments need to be made 
about the validity of the procedure. Since the cc calculations can only be made to agree with the data 



when somewhat lower than average quark masses are used and the pp — ► bb data also need smaller m& 
to agree at NLO, it is reasonable to expect that corrections beyond NLO could be large. Indeed, the 
preliminary HERA-B cross section agrees with the NNLO-NNLL cross section in Ref. [213], suggesting 
that the next-order correction could be nearly a factor of two. Thus the NNLO correction could be nearly 
as large as the NLO cross section. 

Unfortunately, the NNLO-NNLL calculation is not a full result and is valid only near threshold. 
The pp data at higher energies, while not total cross sections, also show a large discrepancy between the 
perturbative NLO result and the data, nearly a factor of three [215]. This difference could be accounted 
for using unintegrated partem densities [216] although the unintegrated distributions vary widely. The 
problem remains how to connect the regimes where near-threshold corrections are applicable and where 
high-energy, small- x physics dominates. The problem is increased for charm where, even at low energies, 
we are far away from threshold. 

We take the most straightforward approach-comparing the NLO calculations to the pp and n~p 
data without higher-order corrections. The cc data are extensive while the bb data are sparse. The cc data 
tend to favor lighter charm masses but the bb data are less clear. A value of mj = 4.75 GeV underpredicts 
the Tevatron results [215] at NLO but agrees reasonably with the average of the ir~p data. However, for 
the HERA-B measurement [217] to be compatible with a NLO evaluation, the b quark mass would have 
to be reduced to 4.25 GeV, perhaps more compatible with the Tevatron results. Therefore, if the NNLO 
cross section could be calculated fully, it would likely be more compatible with a larger quark mass. 
A quantitative statement is not possible, particularly for charm which is far enough above threshold at 
fixed-target energies for the approximation to be questionable [218], see however Ref. [219]. 

If we assume that the energy dependence of the NNLO and higher order cross sections is not 
substantially different from that at LO and NLO, then we should be in the right range for collider energies. 
The theoretical K factors have a relatively weak yfs dependence, < 50% for 20 GeV < yfs < 14 TeV 
[220]. The heavy quark distributions might be slightly affected since the shapes are somewhat sensitive 
to the quark mass but the LHC energy is far enough above threshold for the effect to be small. 

Table 9: Charm and bottom total cross sections per nucleon at 5.5, 8.8 and 14 TeV. 

5.5 TeV 8.8 TeV 14 TeV 

pp pPb pp pPb pp 
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The total pp and pA cross sections at 5.5, 8.8 and 14 TeV, calculated with a code originally from J. 
Smith, are given in Table 9. Shadowing is the only nuclear effect included in the calculation. The cross 



sections do not reflect any possible difference between the pA and Ap rates because they are integrated 
over all X\ and x 2 . However, the different x regions probed by flipping the beams could be effectively 
studied through the dilepton decay channel. Thus we now turn to a calculation of the nuclear gluon 
distribution in pA interactions [73]. We show that the dilepton continuum can be used to study nuclear 
shadowing and reproduces the input shadowing function well, in this case, the EKS98 parameterization 
[27, 28]. To simplify notation, we refer to generic heavy quarks, Q, and heavy-flavored mesons, H. The 
lepton pair production cross section in pA interactions is 



da pA->a+x 

dM a dy a 



J dPpvFpi J d 3 pHd 3 p w 5(M fl -M(p l ,pj))5(y fl -y(p l ,pj)) 

dT H ^ l+X {p H ) dF^+X^ da pA^HH+X 

d 3 pi d 3 P[ d 3 pnd 3 pjj 

X0(Vmin <yi,yj< 2/max)0(</>min < <Pl,<kf < 0max) (71) 

where M(pi,pj) and y(pi,Pj) are the invariant mass and rapidity of the lepton pair. The decay rate, 
dF H ^ l+x (pn) /d 3 pi, is the probability that meson H with momentum pn decays to a lepton I with 
momentum p\. The 9 functions define single lepton rapidity and azimuthal angle cuts used to simulate 
detector acceptances. 

Using a fragmentation function Dq to describe quark fragmentation to mesons, the HH produc- 
tion cross section can be written as 

da pA^HH+X fd 3 p Q d 3 PQ d(T pA^QQ+X r 

x 5® (p H - Zl p Q ) 5 ( V (p^j - Z2 p^) . (72) 

Our calculations proved to be independent of Dq. The hadronic heavy quark production cross section 
per nucleon in Eq. (72) is 



1 da pA^QQ+X p da ij ^ QQ 



(73) 



where the parton densities in the nucleus are related to those in the proton by f A = R A where Rf is 
the EKS98 parameterization [27, 28]. The partonic cross section in Eq. (73) is the differential of Eq. (69) 
at k = 0. Note that the total lepton pair production cross section is equal to the total QQ cross section 
multiplied by the square of the lepton branching ratio. 

We compare the ratios of lepton pair cross sections with the input R A in Fig. 46. The RHIC and 
LHC results are integrated over the rapidity intervals appropriate to the PHENIX and ALICE dilepton 
coverages. The ratio is similar to but smaller than R A at all energies. The higher the energy, the better 
the agreement. The ratio always lies below R A because we include qq annihilation, which decreases with 
energy, but quark and gluon shadowing are different and the phase space integration smears the shadow- 
ing effect relative to R A ((x2), (Q))- The calculated ratio deviates slightly more from Rf{(x 2 ),(Q)) for 
e + e~ than for because the curvature of R A with x is stronger at larger x and, due to the differences 

in rapidity coverage, the average values of x 2 are larger for e + e~. 

The average Q 2 increases with energy and quark mass while the average x 2 decreases with energy. 
At the SPS, 0.14 < (x 2 ) < 0.32, Rf is decreasing. At RHIC, 0.003 < (x 2 ) < 0.012, Rf is increasing 
quite rapidly. Finally, at the LHC, 3 x 10~ 5 < (^2) < 2 x 10~ 4 , Rf is almost independent of x. The 
values of (x 2 ) are typically larger for electron pairs at collider energies because the electron coverage 
is more central than the muon coverage. We have only considered pPb collisions but in Pbp collisions 
X\ is an order of magnitude or more greater than x 2 at forward rapidities such as in the ALICE muon 
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Fig. 46: The ratios of lepton pairs from correlated DD and BB decays in pA to pp collisions at the same energies (solid curves) 
compared to the input Rf at the average X2 and Q (dashed)/ yj (Q 2 ) (dot-dashed) of each M bin. From Ref. [73]. 



spectrometer. Thus doing both pA and Ap collisions would greatly increase the x range over which the 
nuclear gluon distribution could be determined. By assessing cc and bb pairs together, it is possible to 
investigate the evolution of the nuclear gluon distribution with Q 2 , albeit with slightly shifted values of 
x due to the mass difference. 

The importance of such a study as a function of A and y/s is clear, especially when there is cur- 
rently little direct information on the nuclear gluon distribution. It is desirable for the pp, pPb, and Pbp 
measurements to be performed at the same energy as the Pb+Pb collisions at 5.5 TeV even though the 
relative proton/neutron number is negligible in gluon-dominated processes. The effect of shadowing is 
easier to determine if the x value does not shift with energy, particularly because the dilepton measure- 
ment does not provide a precise determination of the heavy quark momenta, only an average. 

As shown in Fig. 46, energy systematics would provide a much clearer picture of the nuclear gluon 
distribution over the entire x range. The PHENIX collaboration has already demonstrated that charm can 
be measured using single leptons [221]. Perhaps it might be possible to obtain the bottom cross section 
the same way as well at higher px- Another handle on the QQ cross sections may be obtained by e/j, 
correlations, since the DY and quarkonium decays will not contribute. At the LHC, uncorrected QQ 
pairs can make a large contribution to the continuum. These pairs, with a larger rapidity gap between the 
leptons, are usually of higher mass than the correlated pairs. However, the larger rapidity gap reduces 
their acceptance in a finite detector. Like-sign subtraction also reduces the uncorrected background. For 
further discussion, see Refs. [72, 73, 79]. 




6.5 Cronin effect in proton-nucleus collisions: a survey of theoretical models 

Alberto Accardi 



6.51 Introduction 

In this short note I will compare available theoretical models for the description of the so-called Cronin 
effect [268] in inclusive hadron spectra in pA collisions. The analysis will be limited to references 
containing quantitative predictions for pA collisions [180,205,269-271]. The observable we study is 
the Cronin ratio, R, of the inclusive differential cross sections for proton scattering on two different 
targets, normalized to the respective atomic numbers A and B: 



R{pt) 



B d(T pA /d 2 p T 
A da pB /d?p T 



In absence of nuclear effects one would expect R(pr)=l, but for A > B a suppression is measured at 
small pt, and an enhancement at moderate pt with R{pr) —> 1 as px — > oo. This behaviour may be 
characterized by the value of three parameters: the transverse momentum p x at which R crosses unity 
and the transverse momentum pj\j at which R reaches its maximum value Rm=R(pm), see Fig. 47. 
These Cronin parameters will be studied in Sec. 6.53. 




Px Pm p T 

Fig. 47: Definition of p x ,Pm, Rm- 



The Cronin effect has received renewed interest after the experimental discovery at RHIC of an 
unexpectedly small R < 1 in Au+Au collisions at ^/s=l30 GeV compared to pp collisions at the same 
energy [272,273]. This fact has been proposed as an experimental signature of large jet quenching, 
suggesting that a Quark-Gluon Plasma was created during the collision [274,275]. However, the ex- 
trapolation to RHIC and LHC energies of the known Cronin effect at lower energies is haunted by large 
theoretical uncertainties which may make the inteipretation such signals unreliable. More light has been 
recently shed on this problem by the RHIC data on dA collisions at y/s = 200 GeV [6-8, 276]. 

Since no hot and dense medium is expected to be created in pA collisions, a pA run at the same 
nucleon-nucleon energy as in AA collisions would be of major importance to test the theoretical models 
and to have reliable baseline spectra for the extraction of novel physical effects. This has been the key in 
establishing the relevance of final state effects (jet quenching) as the the cause of the suppression of the 
Cronin ratio R in Au+Au collision at RHIC [276, 277]. This measurement will also be essential at LHC, 
where at the accessible Bjorken x < 10~ 3 new physics such as gluon saturation may (or may not) come 
into play. 

A further advantage of pA collisions is the relatively small multiplicity compared to AA collisions. 
For this reason, at ALICE, minijets may be observed at px > 5 GeV, see Section 2. As I will discuss in 
Sec. 6.53, the Cronin effect on minijet production may then be a further check of the models. 



6.52 The models 

Soon after the discovery of the Cronin effect [268], it was realized that the observed nuclear enhancement 
of the p^-spectra could be explained in terms of multiple interactions [278,279]. The models may be 
classified according to the physical object which is rescattering (the projectile hadron or its partons), and 
to the "hardness" of the rescattering processes taken into account. Note that a parton is commonly said 
to undergo a hard scattering if the exchanged momentum is of the order or greater than approximately 
1 GeV. However, physically there is no sharp distinction between soft and hard momentum transfer. 
Therefore, I prefer to refer to the so-called two component models of hadron transverse spectra, and call 
a scattering which is described by a power-law differential cross section at large p? hard, and a scattering 
whose cross section is decreasing faster than any inverse power of the transverse momentum at large pt 
soft. In Tables 10 and 11,1 provide a quick comparison of the hadronic and partonic rescattering models. 



Soft hadronic rescattering models [180, 269] 

These models are based on the pQCD collinearly factorized cross section for inclusive particle production 
in pp collisions. In order to better describe the hadron transverse momentum spectra at yfs ~ 20 GeV, 
and in the semi-hard region pt ~ 1 . . . 5 GeV at higher cms energies, one has to also include an intrinsic 
transverse momentum [280, 281] for the colliding parton. The pr spectrum in a pA collision is obtained 
from a Glauber-type model 6 : 

-J^- = KY J F i / P ® Fj/A ® - kl) D h k , (74) 

i,j,k,l 

where the proton and nucleus parton distribution functions are, respectively, 



e -k} T /(k 2 T ) pA {b) e -kj T /(k' T ) Ap {b) 
Fi/P = fi/ P (xi,Q 2 ) A — and F j/A = T A (b)f j/A (x j ,Q 2 ) — -^y- . (75) 



In Eq. (74), da/dt(ij — > kl) is the pQCD parton-parton cross section, the variables with a hat are the 
partonic Mandelstam variables, and a sum over incoming and outgoing parton flavours is performed. 
The proton is considered point-like compared to the target nucleus, and the scattering occurs at impact 
parameter b. The nucleus is described by the Woods-Saxon nuclear thickness function T A (b). In Eq. (75), 
fi/p(A) ( x j Q 2 ) w& the PDFs of the proton (nucleus). Isospin imbalance is taken into account and nuclear 
shadowing is included by the HIJING parametrization [282]. Partons are assumed to have an intrinsic 
transverse momentum with average squared value (k^) pA ^ Ap ) and a Gaussian distribution. Due to the 
hi and kj integrations a regulator mass m=0.8 GeV has been used in the pQCD cross section. Finally, 
D%(z, Q' 2 ) are the fragmentation functions of a parton k into a hadron h with a fraction z of the parton 
momentum. 

Soft proton-nucleon interactions are assumed to excite the projectile proton wavefunction, so that 
when the proton interacts with the next target nucleon its partons have a broadened intrinsic momentum. 
Each proton rescattering is assumed to contribute to the intrinsic momentum broadening in the same way 
so that: 

(k 2 T ) pA {b, y/s) = (4) pp + 5 h A (b, yfi) , (76) 

where {k?f) pp is the parton intrinsic momentum needed to describe hadron transverse spectra in pp colli- 
sions, 5 is the average momentum squared acquired in each rescattering, and 

M^) = { M6 '^ )_1 { ' VA -\l n HD 

In if v A — 1 > n 



integrations are schematically indicated with crossed circle symbols. For details see the original references. 



Table 10: Parameters of the soft hadronic rescattering models [180]. 



model 


hard scales (GeV) 


K 


regul. 


proton intrinsic k T (GeV 2 ) 


n 


average fc T -kick (GeV 2 ) 


nPDF 


Soft 


Q = Q' =pt 
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{k 2 T ) pp = 1.2+0.2a s (? 2 )g 2t 
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z (n \ _ n nrr ln^Q/GcV) 
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0.8 GeV 


{k 2 T ) PP = F{y/s) * 


4 


8 = 0.4 


HIJING 



q 2 = 2siu/(s + t + u) 2 . Parametrization chosen to best reproduce pp data. 

No explicit parametrization is given. Values of (kfypp extracted from a "best fit" to pp data, see Fig. 15 of Ref. [180]. 
These scales used in the computations of Table 12. In Ref. [180] Q = Pt/2 and Q' = pr/2z c . 



represent the average number of collisions which are effective in broadening the intrinsic momentum. 
Both models assume h A to be a function of the number of proton-nucleon collisions u a (b,-\/s) = 
a pp {y/s)T A {b), with a pp the nucleon-nucleon inelastic cross section. However, Ref. [269] takes n=oo, 
while Ref. [180] assumes an upper limit n=4 justified in terms of a proton dissociation mechanism: after 
a given number of interactions, the proton is so excited that it can no longer interact as a whole with the 
next nucleon. I will call the first model simply soft and the second soft-saturated. In both models, target 
nucleons do not rescatter, so that 

(kT)Ap(b, y/s) = {k^jpp ■ 

Further differences between the models are related to the choices of the hard scales, of the If -factor 
which simulates NLO contributions to the parton cross section, and to the parametrizations of {k?f) pp 
and 5 (see Table 10). 



Soft partonic rescatterings: the colour dipole model [270] 

In this model, the particle production mechanism is controlled by the coherence length l c =\fs / '(rn,Nkr), 
where is the nucleon mass and kr the transverse momentum of the parton which fragments in the 
observed hadron. Depending on the value of l c , three different calculational schemes are considered. 
(a) In fixed target experiments at low energy (e.g., at SPS), where l c <C R A , the projectile partons 
interact incoherently with target nucleons. High-py hadrons are assumed to originate mainly from a 
projectile parton which experienced a hard interaction and whose transverse momentum was broadened 
by soft parton rescatterings. The parton is then put on-shell by a single semihard scattering computed 
in factorized pQCD. This scheme I will discuss in detail below, (b) At LHC, where the cm. energy is 
very large and l c 3> Ra, the partons interact coherently with the target nucleons and high-p^ hadrons 
are assumed to originate from radiated gluons. Parton scatterings and gluon radiation are computed in 
the light-cone dipole formalism in the target rest frame, (c) At intermediate energies, like RHIC, an 
interpolation is made between the results of the low- and high-energy regimes discussed above. All the 
phenomenological parameters needed in this model are fixed in reactions other than from pA collisions, 
and in this sense the model is said to be parameter-free. 

In the short coherence length scheme, pQCD factorization is assumed to be valid and Eq. (74) is 
used with parton masses m g = 0.8 GeV and m q = 0.2 GeV for, viz., gluons and quarks. Moreover, 

Fi/v = fi/ P {xi + ^Q 2 )-§^{xM and F j/A = T A {b) f j/A ( Xj , Q 2 ) ^(x, b) . 

Parton rescatterings are computed in terms of the propagation of a qq pair through the target nucleus. 
The final parton transverse momentum distribution dNi/d 2 kiT is written as: 



d 2 h 



iT 
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(78) 



The first bracket in Eq. (78) represents the contribution of the proton intrinsic momentum. The sec- 
ond bracket is the contribution of soft parton rescatterings on target nucleons, expressed through the 



Table 11: Parameters of the soft partonic rescattering model (at short-/ c ) and of the hard partonic rescattering models. 



model 


hard scales 


K 


regulators (GeV) 


intr. kr 


dE/dz 


dipole cross-sect. 


nPDF 


Col. dip. 


Q — Q' — pr 


® 


m 9 =0.8, m q =0.2 


as Soft mod. 


-2.5 GeV/fm 


tTo=23 mb, A=0.288, zo=3-lCr 4 


EKS98 T 


Hard AT 


Q = Q'=P* T 


2 


/i free param. 


no 


no 


computed from pQCD 


no 


HardGV 


Q = Q' = pt 




fi=0A2 b 


as Soft mod. 


no 




EKS98 



* Q=n in Ref. [205]. f Only at large Xj (EMC effect). 18 Factors out in the Cronin ratio. 
b /i determines only the typical momentum transfer in elastic rescatterings. 



phenomenological dipole cross section. For a quark, cr^(r<r, x) = do [l — exp ( — \ r\ Q 2 (x))] with 
Q s =l GeV(x/xo) x / 2 (the parameters ar given in Table 11). For a gluon, a^ g =9/4a^ q is used. The ex- 
pansion of Eq. (78) to zeroth order in a^ q gives the intrinsic kx distribution, dN^> /d 2 kx, of the nucleon 
partons. The first order term represents the contribution of one-rescattering processes, and so on. Partons 
from the target nucleus are assumed not to rescatter because of the small size of the projectile. Energy 
loss of the projectile partons is taken into account by a shift of their fractional momentum proportional to 
the energy of the radiated gluons, given by the product of the average mean path length AL and the en- 
ergy loss per unit length dE/dz [283, 284]. Since nuclear shadowing effects are computed in the dipole 
formalism, see Eq. (78), parton distribution functions in the target are modified only to take into account 
antishadowing at large x according to the EKS98 parametrization [27]. 

By Fourier transforming the dipole cross section we find that the transverse momentum distribu- 
tion of single parton-nucleon scattering is Gaussian, which justifies the classification of this model as 
"soft". However, the single distributions are not just convoluted to obtain a broadening proportional 
to the average number of rescatterings. Indeed, in the second bracket, the rescattering processes have a 
Poisson probability distribution. As a result, the nuclear broadening of the intrinsic momentum is smaller 
than the product of the average number of rescatterings and the single scattering broadening; perhaps a 
dynamical explanation of the assumption that n < oo in Eq. (77), used in the Soft-saturated model [180]. 



Hard partonic rescattering model [205, 271] 

The model of Ref. [205], hereafter labeled "hard AT", assumes that parton rescatterings are responsible 
for the Cronin enhancement. Up to now, it includes only semihard scatterings, i.e., scatterings described 
by the pQCD parton-parton cross section. It is the generalization to an arbitrary number of hard parton 
rescatterings of the early models of Refs. [278, 279, 285] and of the more recent Refs. [286, 287], which 
are limited to a single hard rescattering. As shown in Ref. [279], considering only one rescattering may 
be a reasonable assumption at low energy to describe the gross features of the Cronin effect, but already 
at RHIC energies this might be insufficient to determine the Cronin peak Rm [205]. The AT model 
assumes that the S-matrix for a collision of n partons from the projectile on m partons from the target is 
factorizable in terms of S -matrices of parton-parton elastic-scattering. It also assumes generalized pQCD 
factorization [18, 288]. The result is a unitarized cross section, as discussed in Refs. [205, 289]: 

- E A* ® #f ° D ? + E h,A Ta ^ ® D) . (79) 

* 3 

The first term accounts for multiple semihard scatterings of partons in the proton on the nucleus. In the 
second term, partons of the nucleus are assumed to undergo a single scattering, with dNj/ p /d 2 kT = 
Yli fi/p ® a f/N'- Nuclear effects are included in dN^/^/dPkT, the average transverse momentum dis- 
tribution of a proton parton that suffered at least one semihard scattering. In impact parameter space it 
reads 



dNi l A t^ _ f d2r „-ik 



d 2 k T (b) = / IT 



<T-r 



(80) 



where unitarity is explicitly implemented at the nuclear level, as discussed in Ref. [289]. In Eq. (80), 
a% d (r) = KY,jJd 2 p [l - e'^] ^(ij -» ij) ® f j/p . Moreover, w< = cr^T^fe) is identified with 
the target opacity to the partem propagation. Note that o^f^(r) oc r 2 as r — > and cr^^r) — ► cjm 1 as 
r — > oo. This, together with the similarity of Eqs. (78) and (80), suggests the interpretation of <r^ d (r) 
as a Ziara? dipole cross section, accounting for hard parton rescatterings on target nucleons analogous to 
what does for soft parton rescatterings. Note that no nuclear effects on PDF's are included. However, 
shadowing is partly taken into account by the multiple scattering processes. 

To regularize the IR divergences of the pQCD cross section, a small mass regulator, fi, is intro- 
duced in the parton propagators. It is considered a free parameter of the model that sets the scale at 
which pQCD computations break down. As a consequence of the unitarization of the interaction, due to 
the inclusion of rescatterings, both p x and pu are almost insensitive to fi [205]. Therefore, these two 
quantities are considered reliable predictions of the model 7 . Note, however, that they both depend on the 
cm. energy yfs and on the pseudorapidity rj. On the other hand, Rm is quite sensitive to the IR regula- 
tor. This sensitivity may be traced back to the inverse-power dependence of the target opacity Ui on //, 
Ui oc 1//A where a > 2 is energy and rapidity dependent. The divergence of as \i — > indicates 
the need for unitarization of the parton-nucleon cross section and deserves further study. Therefore, \i 
can be here considered only as an effective scale which simulates nonperturbative physics [295,298], 
the nonlinear evolution of target PDFs [82] and physical effects so far neglected, e.g., collisional and 
radiative energy losses [299]. 

Another model which implements hard partonic rescatterings, labeled "hard GV", is found in 
Ref. [271]. This model is discussed in more detail in the Chapter on jet physics, see Ref. [301]. In the 
hard GV model, the transverse momentum broadening of a parton which experiences semihard rescat- 
terings is evaluated, using Eq. (80), to be 



The IR regulator is fixed to /i = 0.42 GeV and identified with the medium screening mass. It represents 
the typical momentum kick in each elastic rescattering of a hard parton. The factor c and the constant term 
1 are introduced in order to obtain no broadening for px ~ 0, as required by kinematic considerations. 
The average value of the opacity in the transverse plane, to m (0.4/fm)i? J 4 (where Ra is the nuclear 
radius), and the factor c//i 2 = 0.18 are fixed in order to reproduce the experimental data at ^fs = 27 A 
GeV and y/s = 38.8 GeV [291]. With these values the logarithmic enhancement in Eq. (81) is of order 
1 for px ~ 3 GeV. Note that oj and c are assumed to be independent of ^fs and 77. Finally, the transverse 
spectrum is computed by using Eq. (74) and adding to the semihard broadening of Eq. (81) the intrinsic 
momentum of the projectile partons, (k^) = 1.8 GeV 2 : 



with shadowing and antishadowing corrections to target partons from the EKS98 parametrization [27]. 
6.53 Predictions and conclusions 

In Table 12, I list the values of the Cronin parameters pm and Rm computed in various models of pA 
collisions at ^fs = 27 A GeV (for low-energy experiments), 200 GeV (RHIC) and 5.5 TeV (LHC). The 
targets considered in the Cronin ratio are also listed. I do not include p x in the tables since in almost 
all models p x ~ 1 GeV, independent of energy, which lies at the border of validity of the models. 
Uncertainties in the model calculations are included when discussed in the original references, see the 
table footnotes. In the case of the soft-saturated model, the uncertainty due to the choice of shadowing 

7 This result is very different from the conclusion that p x oc po, Ref. [287], based on a single-rescattering approximation. 
Hence p x cannot be used to "measure" the onset of hard dynamics as proposed in that paper. 
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Table 12: Comparison of models of the Cronin effect at r\ — 0. 





model 


charged pions 
pu (GeV) R M 


Ref. 


p M (GeV) 


partons 

Rm 


Ref. 




Soft 


4.0 


1.55* 


[269] 








27.4 GeV 


Soft-saturated 


4.5*; 4.4 


1.46*; 1.46 


[293] 


5.1*; 5.1° 


1.50*; 1.51° 


[293] 


pW j pBe 
data [291] 


Color dipole 


4.5 


1.43±0.08® 


[270] 








Hard AT 








6 ±0.8 + 


l.l b ; 1.4>> 




HardGV 


4 


1.4 


[271] 










Soft 


3.5 


1.35±0.2l 


[269] 








200 GeV 

pAu/pp 


Soft-saturated 


2.9*; 2.7 


1.15*; 1.47 


[293] 


4.4*; 4.2 


1.29*; 1.7O 


[293] 


Color dipole 
Hard AT 
HardGV 


2.7 
3.0 


1.1 
1.3 


[270] 
[271] 


7±lt 


1.25 b ; 1.2>> 






Soft 


3.5 


1.08±0.02* 


[294] 








5500 GeV 

pPb jpp 


Soft-saturated 


2.4*; 2.2 


0.78*; 1.36 


[293] 


4.2*; 4.2 


0.91*; 1.6O 


[293] 


Color dipole 
Hard AT 
HardGV 


2.5 

«40^ 


1.06 
1.05 


[270] 
[271] 


ll±1.3t 


2.1 b ; 1.2^ 





* With HIJING shadowing [282]. Without shadowing. ® Error estimated by varying dE/dz within error bars [292]. 

* Central value with multiple scattering effects only; the error estimated by using different shadowing parametrizations. 

* Numerical errors mainly. b Using /i = 1.5 GeV. ^ Using fi — 0.060 + 0.283 log(\/s), see text. 

* Very sensitive to the shadowing parametrization, see text. 



parametrization is illustrated by giving the results obtained with no shadowing, and with the HIJING 
parametrization [282]. Using the "new" HIJING parametrization [66] would change the mid-rapidity 
results only at LHC energies, where a 15% smaller Cronin peak would be predicted [293] 8 . In the case 
of "hard AT" rescattering, the major theoretical uncertainty is in the choice of fi, as discussed at the end 
of Sec. 6.52. In the table, two choices are presented: (a) an energy-independent value fi = 1.5 GeV, 
which leads to an increasing Cronin effect with energy; (b) fi is identified with the IR cutoff po discussed 
in Ref. [300], in the context of a leading order pQCD analysis of pp collisions. That analysis found po to 
be an increasing function of yfs. By performing a simple logarithmic fit to the values extracted from data 
in Ref. [300], we find \x = Po(y/s) = 0.060 + 0.283 log (-y/s), which decreases the Cronin effect with 
energy. Note that a scale increasing with ^fs also appears naturally in the so-called "saturation models" 
for hadron production in AA collisions [295-297]. In the "hard GV" model at the LHC, the remnants of 
the Cronin effect at pr ~ 3 GeV are overwhelmed by shadowing and the calculation is not considered 
reliable in this region. The value Rm = 1-05 at px — 40 in Ref. [271] is a result of antishadowing in the 
EKS98 parametrization and is not related to multiple initial state scatterings. 

As discussed in the introduction, experimental reconstruction of minijets with ALICE may be 
possible in pA collisions for minijet transverse momenta px > 5 GeV. The p^-spectrum of the partons 
which hadronize into the observed minijet may be obtained by setting D^=S(z— 1) in Eqs. (74) and (79). 
Minijet reconstruction may be very interesting because pQCD computations suffer from large uncertain- 
ties in the determination of the large-z fragmentation functions, where they are only loosely constrained 
by existing data [290]. For this reason I also listed the Cronin parameters for parton production in Ta- 
ble 12. However, jet reconstruction efficiency should be accurately evaluated to assess the usefulness of 
this observable. 

Table 12 shows that there are large theoretical uncertainties in the extrapolation of the Cronin ef- 
fect from lower energies to the LHC. A major source of uncertainty for most of the models is nuclear 
shadowing and antishadowing at small x, see Refs. [52, 53] and Sec. 4. for a detailed discussion and 
comparison of the existing parametrizations. For example, the HIJING parametrizations [66, 282] sug- 

8 Note however that the new parametrization [66], which predicts a much larger gluon shadowing at x < 10 -2 than the 
"old" one [282], seems ruled out by data on the Sn/Ca F2 ratio [52, 53]. 



gest more gluon shadowing than EKS98 [27] at small x < 10~ 2 . At the LHC the small- x region is 
dominant at midrapidity and small to medium px- On the other hand, the HIJING parametrizations give 
less antishadowing than EKS98 at x > 10 _1 , which is the dominant region at large pr for all energies. 
At the LHC, all these effects may lead to up to a factor of 2 uncertainty in the height of the Cronin peak 
Ru- 
in conclusion, a pA run at the LHC is necessary both to test theoretical models of particle pro- 
duction in a cleaner experimental environment than AA collisions, and to make reliable extrapolations 
to AA collisions, which are the key to disentangling known effects from new physics. Since, as we have 
seen, the nuclear effects are potentially large, it would even be preferable to have a pA run at the same 
energy as AA. In addition, the study of the A dependence, or collision centrality, would be interesting 
since it would change the opacity of the target - thus the size of the Cronin effect - in a controllable way. 
Finally, note that the r/-systematics of the Cronin effect has so far been considered only in Ref. [205]. 
However, as also discussed in Ref. [301], given the large pseudorapidity coverage of CMS, the 77 depen- 
dence might be a very powerful tool to understand the effect. It would systematically scan nuclear targets 
in the low-x region and would help test proposed models where the rapidity influences the Cronin effect 
in different ways. 



6.6 Quarkonium total cross sections 

Ramona Vogt 

To better understand quarkonium suppression, it is necessary to have a good estimate of the ex- 
pected yields. However, there are still a number of unknowns about quarkonium production in the pri- 
mary nucleon-nucleon interactions. In this section, we report the quarkonium total cross sections in the 
CEM [226,227] and NRQCD [229], both briefly described in Section 5.5. We also discuss shadowing 
effects on the p^-integrated rapidity distributions. 

Recall that at leading order, the production cross section of quarkonium state C in the CEM is 
a CEM = F C Y, H ds d Xl dx 2 ff( Xl ,Q 2 ) fj(x 2 , Q 2 ) a l3 (s) 5{s - Xl x 2 s) . (82) 

i,j Um Q J 

In apA collision, one of the proton parton densities is replaced by that of the nucleus. Since quarkonium 
production in the CEM is gluon dominated, isospin is a negligible effect. However, since shadowing 
effects on the gluon distribution are large, it could strongly influence the total cross sections. We will 
use the same parton densities and parameters that fit the QQ cross section data, given in Table 9 in 
Section 6.4, to determine Fc for J/ip and T production. 

In the left-hand side of Fig. 48, we show the resulting fits to the forward (xp > 0) J/ip total cross 
sections. These data also include feeddown from Xc and ip' decays. The right-hand plot is the sum of all 
three T S states, including the branching ratios to lepton pairs. 

All the fits are equivalent for y/s = 100 GeV but differ by up to a factor of two at 5.5 TeV. The 
high energy T results seem to agree best with the energy dependence of the MRST calculations with 
m b = 4.75 GeV and 4.5 GeV. 

The pp cross sections obtained for the individual states are shown in Tables 13 and 14 along with 
the values of Fc determined from the fits. Two values are given for the T family. 

The first is obtained from the combined S state fits shown in the right-hand side of Fig. 48, includ- 
ing the branching ratios to lepton pairs, and the second is that of the total T IS state after the branching 
ratios have been extracted. The cross sections given in the tables are, with the exception of the total J/ ijj 
cross section including feeddown, all for direct production. To obtain these direct cross sections, we use 




Fig. 48: Forward J/tp (left) and combined T + T' + T" inclusive (right) cross sections calculated to NLO in the CEM. On the 
left hand side, the solid curve employs the MRST HO distributions with m c = Q/2 = 1.2 GeV, the dashed, MRST HO with 
m c = Q = 1.4 GeV, the dot-dashed, CTEQ 5M with m c = Q/2 = 1.2 GeV, and the dotted, GRV 98 HO with m c = Q = 1.3 
GeV. On the right hand side, the solid curve employs the MRST HO distributions with m b — Q = 4.75 GeV, the dashed, 
m b = Q/2 = 4.5 GeV, the dot-dashed, m b = 2Q = 5 GeV, and the dotted, GRV 98 HO with m b = Q = 4.75 GeV. 



the production ratios given by data and branching ratios to determine the relative cross sections, as in 
Refs. [224,231]. 

The direct J/ijj and T rapidity distributions in pPb interactions at 5.5 TeV/nucleon are shown in 
Fig. 49 for all the fits along with the pVblpp ratios at the same energy. Note how the average rapidity 
in the pPb distributions is shifted to negative rapidities due to the low x shadowing at positive rapidity. 
The 'corners' in the J/ip rapidity distributions at \y\ ~ 4 are equivalent to x ~ 10~ 5 , the lowest x for 
which the MRST and CTEQ5 densities are valid. The pVblpp ratios show antishadowing at large negative 
rapidities, large x<i, and 30 — 40% shadowing at large positive rapidities for the T and J/ip respectively 
(small X2). At y ~ 0, the effect is ~ 30% for J/ij; and ~ 20% for T. We plot pPb/pp here but the Vbplpp 
ratio would simply be these results reflected around y = 0. 

The NRQCD total cross sections, calculated with the CTEQ3L densities [303], are also shown in 
Tables 13 and 14. The parameters used for the NRQCD calculations were determined for fixed-target 
hadroproduction in Ref. [302] with m c = 1.5 GeV, m& = 4.9 GeV and Q = Iuiq. Calculations of 
shadowing effects predicted for NRQCD can be found in Ref. [174]. The dependence on the shadow- 
ing parameterization is also discussed. It was found that the EKS98 parameterization [27, 28] predicted 
weaker shadowing than the HPC parameterization [173], mainly because of the larger gluon antishad- 
owing in the EKS98 parameterization [174]. The predicted shadowing effects were similar for the CEM 
and NRQCD approaches. 

Another approach to quarkonium production has been presented recently. Ref. [304] argues that 
the existing quarkonium data suggest the presence of a strong color field, V, surrounding the hard QQ 
production vertex in hadroproduction. This field is assumed to arise from the DGLAP evolution of the 
initial colliding partons. The effect of rescatterings between the quark pair and this comoving field was 
shown to be in qualitative agreement with the main features of the pp and pA data [304]. Quantitative 
predictions within this model are difficult because the form of the field F is unknown but some predictions 



Table 13: The production fractions obtained from fitting the CEM cross section to the J/ip forward cross sections (x F > 0) 
and y = cross sections as a function of energy. The parton distribution function (PDF), charm quark mass (in GeV), and 
scales used are the same as those obtained by comparison of the cc cross section to the pp data. The charmonium cross sections 
(in /ib) obtained from the CEM fits for NN collisions at 5.5 TeV are also given. The production fractions are then multiplied 
by the appropriate charmonium ratios determined from data. Note that the Xc cross section includes the branching ratios to 
J/xjj. The last row gives the NRQCD cross sections with the parameters from Ref. [302]. 



PDF 
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Q/m c 




tot 

a J/lp 


dir 




dir 


MRST HO 


1.2 
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0.0144 


30.6 


19.0 


9.2 


4.3 


MRST HO 


1.4 
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0.0248 


20.0 


12.4 


6.0 


2.8 


CTEQ 5M 


1.2 


2 


0.0155 


36.0 


22.2 


10.8 


5.0 


GRV 98 HO 


1.3 


1 


0.0229 


32.1 


19.8 


9.6 


4.5 


CTEQ 3L 


1.5 


2 




83.1 


48.1 


27.6 


13.6 



Table 14: The production fractions obtained from fitting the CEM cross section to the combined T cross sections to muon pairs 
at y — as a function of energy. The PDF, bottom quark mass (in GeV), and scales used are the same as those obtained by 
comparison of the bb cross section to the ir~p data. The direct bottomonium cross sections (in nb) obtained from the CEM fits 
for AW collisions at 5.5 TeV in each case above. The production fractions for the total T are multiplied by the appropriate 
ratios determined from data. In this case, the \b cross sections do not include any branching ratios to lower S states. The last 
row gives the NRQCD cross sections with the parameters from Ref. [302]. 
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assuming an isotropic and transversely polarized field were given in Ref. [305] . 

So far we have not discussed nuclear absorption in cold matter. The A dependence of J /if), ip' and 
the T S states have all been measured in pA interactions. (See Ref. [225] for a compilation of the data.) 
However, there is still no measurement of the P state A dependence. The CEM makes no distinction 
between singlet and octet QQ production so that e.g. the J/ip, ip' and Xc A dependencies should all be the 
same. The NRQCD approach, on the other hand, predicts that direct ip' and J/ ijj are 60 — 80% octet but 
the Xc states are 95% singlet [307]. If octets and singlets are absorbed differently, the \c A dependence 
should be quantitatively different than the quarkonium S states, see Ref. [307] for a detailed discussion 
at fixed-target energies. In addition, the hard rescattering approach leads to a larger Xci A dependence 
at high xp than for the Xc2 [308]. Thus, measuring the Xc A dependencies could help distinguish the 
underlying quarkonium production mechanism. 

NA60 [241] and HERA-B [309] will measure the Xc A dependence for the first time at y/s = 29.1 
and 41.6 GeV respectively. At the LHC, it may also be possible to extract the Xb A dependence in pA 
interactions as well. Comparing the A dependence of all quarkonium states at the LHC with earlier 
results will also provide a test of the energy dependence of the absorption mechanism, suggested to be 
independent of energy [310]. Since the absorption mechanism depends on impact parameter, centrality 
measurements of the A dependence could also be useful in the comparison to AA interactions. 

'Thanks to P. Hoyer and S. Peigne for comments on their work. See also [306]. 



Fig. 49: The J/tp and T rapidity distributions in pPb collisions (left) and the pPblpp ratios (right). The J/tp calculations 
employ the MRST HO distributions with m c = Q/2 = 1.2 GeV (solid), MRST HO with m c = Q = 1.4 GeV (dashed), 
CTEQ 5M with m c = Q/2 = 1.2 GeV (dot-dashed) and GRV 98 HO with m c = Q = 1.3 GeV (dotted). The T calculations 
employ the MRST HO distributions with m b = Q = 4.75 GeV (solid), m 6 = Q/2 = 4.5 GeV (dashed), m b = 2Q = 5 GeV 
(dot-dashed), and GRV 98 HO with m b = Q = 4.75 GeV (dotted). 



6.7 px broadening in quarkonium production 

Ramona Vogt 

The broadening of the quarkonium transverse momentum distribution was first observed by NA3 
[311]. Additional data at similar energies showed that the change in the J/tp px slope with A was larger 
than for the Drell-Yan slope [312]. The T px distribution was shown to broaden still more than the J/ip> 
[313]. 

The pt broadening effect was explained as multiple elastic scattering of the initial partem before 
the hard collision [314-316]. The difference between the J/ip and Drell-Yan broadening can be at- 
tributed to the larger color factor for gluons (J/ip) than for quarks (Drell-Yan). The difference between 
T and J/tp broadening was more difficult to explain. The hard rescattering model of quarkonium produc- 
tion [304, 305] is able to describe the effect [317]. By Coulomb rescattering in the nucleus, the QQ pair 
acquires additional transverse momentum kx with an upper limit given by tuq. Indeed, for kx > rriQ 
the color octet QQ pair can be 'probed' and cannot rescatter as a gluon any longer. Since mb > m c , the 
T transverse momentum broadening should be larger than that of the J/ip [317] 10 . 

In this section, we investigate the combined effect of shadowing and transverse momentum broad- 
ening on J/ip and T production. The double differential distribution for QQ pair production in pp 
collisions is 

s 2f^M = J2f ^^ci 2 A ;T1 d 2 fc r25p (4 1 )5p(42)/f(^,Q 2 )/J(^,Q 2 )^(^)^^P(83) 

cltau J x\ X2 ' dtdu 
hj 

10 Thanks to P. Hoyer and S. Peigne for comments on their work. See also [306]. 



where Jk is a kinematics-dependent Jacobian and s 2 d 2 Oij(s) / didu is the partonic cross section. These 
partonic cross sections are difficult to calculate analytically beyond LO and expressions only exist near 
threshold for s 2 d 2 cfij(s,i,u)/didu, see e.g. Ref. [213]. However, the exclusive QQ pair NLO [214] 
code calculates double differential distributions numerically. 

The bare quark px distributions have often been sufficient to describe inclusive D meson produc- 
tion. However, the QQ pair distributions, particularly the pair px and the azimuthal angle distributions 
at fixed target energies, are broader than can be accounted for by the bare quark distributions without 
fragmentation and intrinsic transverse momentum kx- Heavy quark fragmentation, parameterized from 
e + e~ — ► HX, which reduced the average momentum of the heavy quark in the final state hadron, is 
included. This effect, along with momentum broadening with (k T ) = 1 GeV 2 , can describe the QQ 
measurements [318]. We include the broadening here but not fragmentation since, in the CEM, both the 
Q and Q hadronize together to make the quarkonium state, inconsistent with independent fragmentation. 

The intrinsic kx of the initial partons has been used successfully to describe the px distributions of 
other processes such as Drell-Yan [280] and hard photon production [281]. Some of the low px Drell-Yan 
data have subsequently been described by resummation to all orders but the inclusion of higher orders 
has not eliminated the need for this intrinsic kx- 

The implementation of intrinsic kx in the QQ code in Ref. [214] is not the same as in other pro- 
cesses because divergences are canceled numerically. Since including additional numerical integrations 
would slow this process, the kx kick is added in the final, rather than the initial, state. In Eq. (83), the 
Gaussian g p (k T ), 

9p( k r) = /, 2 v exp(-4/(4) pp ) , (84) 

X' PP 

with (k T ) pp = 1 GeV 2 [318], assumes that the x and kx dependencies completely factorize. If this is 
true, it does not matter whether the k T dependence appears in the initial or final state, modulo some 
caveats. In the code, the QQ system is boosted to rest from its longitudinal center of mass frame. 
Intrinsic transverse momenta of the incoming partons, kxi and kx2, are chosen at random with k T1 and 
k T2 distributed according to Eq. (84). A second transverse boost out of the pair rest frame changes the 
initial transverse momentum of the QQ pair, px to px + kxi + kx2- The initial kx of the partons could 
have alternatively been given to the entire final-state system, as is essentially done if applied in the initial 
state, instead of to the QQ pair. There is no difference if the calculation is only to LO but at NLO a light 
parton may also appear in the final state. In Ref. [319], it is claimed that the difference between adding 
the kx to the initial or final state is rather small if k T < 2 GeV 2 . However, this may not turn out to be 
the case for nuclei [320]. 

The average intrinsic kx is expected to increase in pA interactions. This broadening is observed 
in Drell-Yan [312, 313], J/ip [311], and T [313] production and has been used to explain high px pion 
production in nuclear interactions [321]. We follow the formulation of Ref. [321], similar to Refs. [314- 
316] for J/ip and Drell-Yan production, where the kx broadening in pA interactions is 

(k 2 T ) pA = (k 2 T ) pp + - 1)A 2 (Q) (85) 

with (k T ) pp = 1 GeV 2 . The number of collisions in apA interaction, v, averaged over impact parameter, 
is [322, 323] 

fd 2 bT 2 (b) 3 

{U) = aNN J d 2 bT A (b) = 2 aNNpoRA (86) 



where T A (b) = f dzp A (b, z) is the nuclear profile function, ct^n is the inelastic nucleon-nucleon scat- 
tering cross section, po is the central nuclear density, and R A is the nuclear radius. Our calculations 



assume = 1.2A 1 / 3 and p = 0.16/fm 3 . The strength of the nuclear broadening, A 2 , depends on Q, 
the scale of the interaction [321] 

^ — iT^, ^ • 

Thus A 2 is larger for bb production than cc production. This empirically reflects the larger kr broadening 
of the T [313] relative to the J/ip [311]. We evaluate A 2 (Q) at Q = Iiuq for both charm and bottom 
production. We find ({v) - 1)A 2 (Q) = 0.35 GeV 2 for cc and 1.57 GeV 2 for bb production in pA 
collisions at b = and A = 200. We can change the centrality by changing (u) — 1. 
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Fig. 50: The energy dependence of {pt)pp for J/tp (left) and T (right) production assuming {k%) pv = 1 GeV 2 . The open 
circles are the calculated energies at ^/s = 0.2, 5.5, 8.8, and 14 TeV. 

In Fig. 50, we show the increase in {p 2 ^) for pp collisions at heavy ion collider energies using the 
MRST parton densities with m c = 1.2 GeV and nif, = 4.75 GeV from Table 9. The average pj, is 
obtained by integrating over all pt and rapidity. Increasing itiq increases {p\) pp - At these energies, the 
difference between (pfypp with and without {k?^) pp = 1 GeV 2 is small, a factor of ~ 1.3 for J/ip and 
~ 1.05 for T at 5.5 TeV. 
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Fig. 51: The pr distributions with xf > for J/ip (left) and T (right) production in pA interactions at 5.5 TeV/nucleon. 
The solid curve is the bare distribution, the dashed curve employs (k^) pp = 1 GeV 2 , and the dot-dashed curve uses a (k^) p A 
appropriate for Pb. 



Figure 51 compares the bare distributions, (k^,) pp = 1 GeV 2 , and (kfypPb calculated using 
Eq. (85), integrated over all rapidity. Without intrinsic transverse momentum, the value of the cross 
section in the lowest px bin is less than zero due to incomplete cancelation of divergences, as discussed 
in Ref. [324]. This point is left out of the plots. It is interesting to note that, while shadowing is a small 
effect at large pr, the average p\ increases with A, even without any broadening included. This arises 
because shadowing is a 25-50% effect on J/ip production at low pr, increasing with A, and a 10-20% 
effect on the T. The resulting increase in (pj) from pp to pPb is ~ 1.2 for J/ip and ~ 1.13 for T. 
Including (k^) pp = 1 GeV 2 increases (p\) over all px as well, but only by ~ 1.15 for both J/ip and T. 
Including the A dependence of the broadening as in Eq. (85) gives about the same increase in (p T ) as 
the bare distributions. 

Note that if a pt cut was imposed, e.g. p? > 5 GeV, (p 2 ^) increases significantly but the effect of 
broadening is greatly reduced. In fact, (pj) actually decreases with increasing (k^) p A when such a cut 
is imposed. This feature is clear from an inspection of Fig. 51. 

The energy and A dependence of A(pj,) p A = (Pt)pA — (Pt)pp * s gi yen hi Fig. 52 for y/s = 
200 GeV/nucleon at RHIC as well as 5.5 and 8.8 TeV/nucleon at the LHC. We find that A(p\) pA for 
J/ip and T increase by factors of ~ 3 and ~ 5.3 over yfij = 19.4 GeV [311] and 38.8 GeV [313] 
respectively. Thus, in this calculation, the increase in A(p^) p A with ^/s is stronger for T than J/ip. 
Such a phenomenological parameterization should soon be tested at RHIC for the J /if). 




Fig. 52: The y/s and A dependencies of A(py) p A for J/ip (left) and T (right) production in pA interactions. The dashed curve 
is at yfs = 200 GeV/nucleon, the solid at 5.5 TeV/nucleon and the dot-dashed at 8.8 TeV/nucleon. We use A =0, Si, I and Au 
at 200 GeV/nucleon, A =Q, Ar, Kr, and Pb at 5.5 TeV/nucleon and A =Pb at 8.8 TeV/nucleon. 



7. NOVEL QCD PHENOMENA IN pA COLLISIONS AT THE LHC 

Leonid Frankfurt and Mark Strikman 

7.1 Physics motivations 

An extensive discussion of the novel pA physics at the LHC was presented in Ref. [325]. Here we 
summarize few key points of this analysis and extend it to consider effects of the small x high-field 
regime. 

The proton-ion collisions at the LHC will be qualitatively different from those at fixed target 
energies. In soft interactions two prime effects are a strong increase of the mean number of "wounded" 
nucleons - from ~ 6 to ~ 20 in pA central collisions (due to increase of <7i ne i(i\W)) and a 50 % increase 
of the average impact parameters in pp collisions (as manifested in the shrinkage of the diffractive peak 
in the elastic pp collisions with increasing energy) making it possible for a proton to interact, often 
simultaneously, with nucleons at the same longitudinal coordinate but different impact parameters. At 
the same time, the geometry of pA collisions at the LHC should be closer to that in classical mechanics 
due to a strong increase of the pp total cross section and smaller fluctuations of the interaction strength 
in the value of the effective NN cross section. (The dispersion of the fluctuations of the cross section 
is given by the ratio of the inelastic and elastic diffractive cross sections at t = 0, which is expected 
to be a factor of 3 ... 4 smaller at the LHC relative to RHIC.) Among other things this will lead to a 
much weaker increase of the inelastic coherent diffractive cross section with A relative to fixed target 
energies [327]. In hard interactions, the prime effect is a strong increase of the gluon densities at small 
x in both projectile and ion. Gluon densities can be measured at the LHC down to x oc 10~ 6 . . . 10" 7 
where pQCD will be probed in a new domain [325]. As a consequence of blackening of the interaction 
of partons from the proton with nuclear partons at small xa, one should expect a disappearance of low 
Pt hadrons in the proton fragmentation region and, hence, a strong suppression of nonperturbative QCD 
effects 11 . Thus the high pt proton fragmentation region is a natural place to search for unusual QCD 
states of quark-gluon matter which should be long lived due to the Lorentz time dilation. In spite of 
a small coupling constant, gluon densities of heavy ions become larger over a wide range of x and 
Q 2 and impact parameters than that permitted by the conservation of probability within the leading- 
twist approximation [157,298,325,329,333]. Thus a decrease of a s with Q 2 is insufficient for the 
applicability of leading-twist pQCD to hard processes in the small x regime where gluon fields may 
be strong enough for the nonperturbative QCD vacuum to become unstable and to require modification 
of QCD dynamics including modification of the QCD evolution equations. The possibility for a novel 
QCD regime and therefore new phenomena are maximized in this QCD regime 12 . At the same time 
interpretation of these novel effects would be much more definitive in pA collisions than in heavy ion 
collisions. 

New classes of strong-interaction phenomena in the short-distance regime, near the light cone, 
which could be associated with perturbative QCD or with the interface of the perturbative and nonpertur- 
bative QCD regions are much more probable in pA than pp collisions. The parton and energy densities 
may be high enough to "burn away" the nonperturbative QCD vacuum structure in a cylinder of ra- 
dius ~ 1.2 fm and length ~ 2Ra, leaving behind only a dense partonic fluid governed by new QCD 
dynamics with small a s but large effective short-distance couplings. If there is a large enhancement 
of hard-collision, gluon-induced processes, then there should be an enhancement of hard multi-parton 
collisions, like production of several pairs heavy-flavor particles [325]. The study of the x, p\ and A- 
dependence of multi-parton interactions, and of forward charm and bottom production (in the proton 

n A similar phenomenon of the low pr suppression of another quantity, the nuclear light cone wave function at small x, is 
discussed in the saturation models [157, 298]. 

12 It is important to distinguish between parton densities which are defined within the leading-twist approximation only [328] 
and may increase with energy forever and physical quantities - cross sections/structure functions - whose increase with energy 
at a given impact parameter is restricted by probability conservation. 



fragmentation region) should be especially useful for addressing these questions. 

We estimate the boundary of the kinematical regime for the new QCD dynamics to set in in pA 
collisions based on the LO formulae for the dipole cross section and the probability conservation. For 
central impact parameters with A ~ 200 and Q 2 ~ 25 GeV 2 , we find xa < 5 • 10~ 5 for quarks and 
x A < 10~ 3 forgluons [329]. 

7.2 Measurements of nuclear parton densities at ultra small x 

One of the fundamental issues in high energy QCD is the dynamics of hard interactions at small x. 
Depending on the Q 2 , one investigates either leading-twist effects or strong color fields. Since the cross 
sections for hard processes increase roughly linearly with A (except at very small x and relatively small 
Q 2 where the counting rates are high anyway), even short runs with nuclear beams will produce data 
samples sufficient to measure parton distributions with statistical accuracy better than 1% practically 
down to the smallest kinematically allowed x: x ~ 10~ 7 for quarks and x ~ 10~ 6 for gluons [330]. The 
systematic errors on the ratios of nuclear and nucleon structure functions are also expected to be small 
since most of these errors cancel in the ratios of the cross sections, cf. Ref. [90]. The main limitations on 
such studies are the forward acceptance and the transverse momentum resolution of the LHC detectors. 

There are several processes which could be used to probe the small x dynamics in pA collisions. 
They include the Drell-Yan pairs, dijets, jet + photon, charm production, diffractive exclusive 3-jet pro- 
duction, and multijet events. It is worth emphasizing that, for these measurements, it is necessary to 
detect jets, leptons and photons at rapidities which are close to the nucleon rapidity. In this y-range, 
the accompanying soft hadron multiplicities are relatively small, leading to a soft particle background 
comparable to that in pp scattering. In fact, for \y max — y\ < 2 — 4 the background level is likely to 
be significantly smaller than in pp collisions due to suppression of the leading particle production in pA 
scattering, see Ref. [325]. 

Measurements of dimuon production appear to be more feasible at forward angles than hadron 
production due to small muon energy losses. In the kinematics where the leading-twist dominates, one 
can study both the small- x quark distributions at relatively small px, pt(^ + aO <C M(fi + /j 1 ^), and the 
gluon distribution at pr{^ + ~ M(/j, + fj,~) [194]. Note also that in the forward region, the dimuon 
background from charm production is likely to be rather small since, in average, muons carry only 1/3 
of the charm decay momentum. 

7.3 Signals for onset of the black body limit 

A fast increase of nPDFs for heavy nuclei with energy contradicts probability conservation for a wide 
range of parton momenta and virtualities for the LHC energy range. On the contrary, the rapid, power- 
like growth of the nucleon structure function observed at HERA may continue forever because of the 
diffuse edge of a nucleon [331]. At present, new QCD dynamics in the regime of high parton densities 
where multi-parton interactions are inhibited remains a subject of lively debate. Possibilities include a 
pre-QCD black-body regime for the structure functions of heavy nuclei [332], a QCD black-body regime 
for hard processes with heavy nuclei [331], saturation of parton densities at a given impact parameter 
[157], a color glass condensate [298], turbulence and related scaling laws. To estimate the pattern of 
the novel phenomena we shall use the formulae of the black-body regime because they are generic 
enough to illustrate the expectations of the onset of new QCD dynamics within the existing theoretical 
approaches 13 . 

High parton density effects can be enhanced in nuclei because of their uniform density and because 

13 In the black-body regime, a highly virtual probe interacts at the same time with several high kr partons from the wave 
function of a hadron [331]. This is qualitatively different from the asymptotic freedom limit where interactions with only one 
highly virtual parton are important. Hence the interacting parton is far from being free when probed by an object with a large 
coherence length, l c = 1/2itinx. 



it is possibile to select scattering at central impact parameters. Indeed, in the small x regime, a hard 
collision of a parton of the projectile nucleon with a parton of the nucleus occurs coherently with all 
the nucleons at a given impact parameter. The coherence length l c rs Xjlm^x by far exceeds the 
nuclear size. In the kinematic regime accessible at the LHC, lc can be up to 10 5 fm in the nuclear rest 
frame. Here the large coherence length regime can be visualized in terms of the propagation of a parton 
through high density gluon fields over much larger distances than is possible with free nucleons. In the 
Breit frame, this corresponds to the fact that small x partons cannot be localized longitudinally within 
the Lorentz-contracted thickness of the nucleus. Thus low x partons from different nucleons overlap 
spatially, creating much larger parton densities than in the free nucleon case. This leads to a large 
amplification of multi-parton hard collisions, expected at small x in QCD, see e.g. Refs [90, 325]. The 
eikonal approximation [26, 150] and constraints from the probability conservation [333, 334] indicate 
that new QCD dynamics should be revealed at significantly larger x in nuclei than in the proton and 
with more striking effects. Recently, extensive studies of hard diffraction were performed at HERA. For 
review see Ref. [335]. Information obtained in these studies allows to evaluate the leading-twist nuclear 
shadowing which turns out to be large both in the gluon and quark channels and significant up to very 
large Q 2 . For a recent discussion see Ref. [329]. This leading-twist shadowing tames the growth of 
the nPDFs at small x. However, even such large shadowing is still insufficient to restore probability 
conservation, leaving room for violation of the leading-twist approximation. This violation may be 
more serious than merely an enhancement of higher twist effects since the entire interaction picture can 
qualitatively change. The common wisdom that a small a s guarantees asymptotic freedom, confirmed 
experimentally in many hard processes, is violated at small x. This phenomenon can be considered in 
terms of effective violation of asymptotic freedom in this limit. This regime is expected to be reached 
in heavy nuclei at x values at least an order of magnitude larger than in the proton. Hence in the LHC 
kinematics, the region where novel QCD phenomena are important will extend at least two orders of 
magnitude in x\ over a large range of virtualities. 

The space-time picture of the black-body limit in pA collisions can be seen best in the rest frame of 
the nucleus. A parton belonging to the proton emits a hard gluon (virtual photon) long before the target 
and interacts with the target in a black regime releasing the fluctuation, a jet/qq (a Drell-Yan pair). This 
leads to a qualitative change in the picture of pA interactions for the partons with x p and pr satisfying 
the condition that 

x A = (88) 

XpSNN 

is in the black-body kinematics for the resolution scale px < p t u (xa)- Here p t u (xa) is the maximum 
Pt at which the black-body approximation is applicable. Hence by varying xa for fixed x v and looking 
for the the maximum value of p T hl we can determine p T hl (x p ). In the kinematics of the LHC, Q 2 
4(p T h1 ) 2 can be estimated by using formulae derived in Ref. [333]. At x p ~ 0.3 Q 2 may be as large 
as (p T hl (x p )) 2 = 15 GeV 2 . All the partons with such x p will obtain pr(jet) ~ Pj? l (x p ), leading to 
multi-jet production. The black body regime will extend down in x p with increasing y / sj\F/v. For the 
LHC with pr < 3 GeV/c, this regime may cover the whole region x p > 0.01 where of the order ten 
partons reside. Hence, in this limit, most of the final states will correspond to multi-parton collisions. 
For pt < 2 GeV/c the region extends to x p > 0.001. At the LHC, collisions of such partons correspond 
to central rapidities. Dynamics of the conversion of high-py partons with similar rapidities to hadrons is 
a collective effect requiring special consideration. 

Inclusive observables. For the total cross sections, logic similar to that employed in 7* - nucleus 
scattering [331] should be applicable. In particular, for the dimuon total cross section, 

da{pA^ n + n~ + X) _ 4ira 2 cm K(x A ,Xp,M 2 ) 2 M 2 2 f x 



[ F 2p ( Xp , Q 2 )^2ttR 2 a In , (89) 



dx A dx p 9 M 2 

for large but not too large dimuon masses, M. In particular we estimate M 2 (bbl)( X A = 10~ 7 ) « 
60 GeV 2 . Here the iv'-factor has the same origin as in the leading-twist case but it should be smaller 



since it originates from gluon emissions only from the partem belonging to the proton. We define x as 
the maximum x for which the black body limit is valid at the resolution scale M. Hence the expected 
M 2 dependence of dimuon production is qualitatively different than in pp scattering where scattering 
at large impact parameters may mask the contribution from the black-body limit. In this limit, the xa 
dependence of the cross section becomes weak. 

The study of the dimuon pt distribution may also signal the onset of the black-body regime. As in 
the case of leading partons in deep inelastic scattering [331], a broadening of the dimuon px distribution 
is expected in the black-body limit compared to the DGLAP expectations. See Ref. [336] for a calculation 
of this effect in the color glass condensate model. 

At still lower xa, the black-body formalism will probably overestimate the cross section because 
the interaction the sea quarks and gluons in the proton will also be in the black-body limit. 

The onset of the black body regime will also lead to important changes in hadron production such 
as a much stronger drop of the xp spectrum in the proton fragmentation region accompanied by pr 
broadening and enhancement of hadron production at smaller rapidities, see the discussion in section 
7.43. 

It is worth emphasizing that onset of the black-body limit for ultra small xa will places limits on 
using forward kinematics for PDF measurements at larger xa- At the very least it would be necessary 
to satisfy the condition pr(jet) S> p T u (x p ). Indeed, as soon as large enough x p are used for the nPDF 
measurements one would have to take into account that such collisions with nuclear partons are always 
accompanied by significant px broadening, due to the interaction with "the black component" of the 
nuclear wave function. At a minimum, this would lead to significant broadening of the pt distribution 
of the produced hard system. One can also question the validity of the leading-twist expansion for the 
PT-integrated cross sections. To investigate these phenomena, the ability to make measurements at fixed 
xa and different x p would be very important. 

7.4 Mapping the three-dimensional nucleon parton structure 

Systematic studies of hard inclusive processes during the last two decades have led to a reasonably good 
understanding of the single parton densities in nucleons. However, very little is known about multiparton 
correlations in nucleons. These correlations can provide critical new insights into strong interaction 
dynamics, as well as discriminate between models of the nucleon. Such correlations may be generated, 
for example, by fluctuations in the transverse size of the color field in the nucleon leading, via color 
screening, to correlated fluctuations of the gluon and quark densities. 

QCD evolution is a related source of correlations since selection of a parton with a given x and 
Q 2 may lead to a local enhancement of the parton density in the transverse plane at different x values. 
Also, practically nothing is known about possible correlations between the transverse size of a particular 
configuration in the nucleon and the longitudinal distribution of partons in this configuration. 

7.41 Multi-jet production and double parton distributions 

It was recognized already more than two decades ago [337] that the increase of parton densities at small 
x leads to a strong increase in the probability of NN collisions where two or more partons from each 
nucleon undergo pairwise independent hard interactions. Although multijet production through double- 
parton scattering was investigated in several experiments [207, 208, 338] at pp and pp colliders, the inter- 
pretation of the data was hampered by the need to model both the longitudinal and the transverse partonic 
correlations simultaneously. Studies of pA collisions at the LHC will provide a feasible opportunity to 
study the longitudinal and transverse partonic correlations in the nucleon separately as well as to check 
the validity of the underlying picture of multiple collisions. 

The simplest case of a multi-parton process is double parton collisions. Since the momentum 
scale pt of a hard interaction corresponds to much smaller transverse distances, ~ l/pr, in coordinate 



space than the hadronic radius, in a double-parton collision the two interaction regions are well separated 
in transverse space. Also, in the center of mass frame, pairs of partons from the colliding hadrons are 
located in pancakes of thickness less than (1/xi + l/x2)/p cm - Thus two hard collisions occur almost 
simultaneously as soon as x\ and X2 are not too small. Hence there is no cross talk between two hard 
collisions. A consequence is that the different partem processes add incoherently in the cross section. 
The double partem scattering cross section, proportional to the square of the elementary parton-parton 
cross section, is therefore characterized by a scale factor with dimension of inverse length squared. The 
dimensional quantity is provided by nonperturbative input, namely by the multi-parton distributions. 
In fact, because of the localization of the interactions in transverse space, the two pairs of colliding 
partons are aligned in such a way that the transverse distance between the interacting target partons 
is practically the same as the transverse distance between the projectile partons. The double parton 
distribution, F(x, x', b), is therefore a function of two momentum fractions and of the transverse distance, 
b. Although T also depends on the virtualities of the partons, Q 2 and Q' 2 , to make the expressions 
more compact, this Q 2 dependence is not explicitly expressed. Hence the double-parton scattering cross 
section for the two "2 — > 2" parton processes a and (3 in an inelastic interaction between hadrons (or 
nuclei) a and b can be written as: 

ao(a,P) = y J r a (xi,X2]b)a a (xi,x[)a j 3(x2,X2)T b (x[,X2]b)dxidx' 1 dx2dx2d 2 b (90) 

where m = 1 for indistinguishable and m = 2 for distinguishable parton processes. Note that though 
the factorization approximation of Eq. (90) is generally accepted in the analyses of the multijet processes 
and appears natural based on the geometry of the process, no formal proof exists. As we will show below, 
the study of the ^4-dependence of this process will be a stringent test of this approximation. 

To simplify the discussion, we neglect small non-additive effects in the parton densities, a rea- 
sonable approximation for 0.02 < x < 0.5. In this case, we have to take into account only b- space 
correlations of partons in individual nucleons. 

There are thus two different contributions to the double-parton scattering cross section (A = 
a,p = b): ctd = a x D + a 2 D . The first one, <rf\ the interaction of two partons from the same nucleon, is 
the same for nucleons and nuclei except for the enhancement of the parton flux. The corresponding cross 
section is 



is the nuclear thickness as a function of the impact parameter B of the hadron-nucleus collision. 

The contribution to T^x^, x' 2 , b) from partons in different target nucleons, a 2 D , can be calcu- 
lated solely from the geometry of the process by observing that the nuclear density does not change 
within the transverse scale (b) <C Ra- It rapidly increases with A as / T\{B)d 2 B. Using infor- 
mation from the CDF double scattering analysis [207, 208] on the mean transverse separation of par- 
tons in a nucleon, one finds that the contribution of the second term should dominate in pA collisions: 
a \>l a \) ~ 0.68(yl/12) 39 for A > 40 [339]. Hence one expects a stronger than oc A increase of multijet 
production in pA collisions at the LHC. Measurements with a range of nuclei would probe the double- 
parton distributions in nucleons and also check the validity of QCD factorization for such processes. 
Factorization appears natural but is so far not derived in pQCD. An important application of this pro- 
cess would be the investigation of transverse correlations between the nuclear partons in the shadowing 
region. A study of these correlations would require selection of both nuclear partons in the shadowing 
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region, xa < ~ 10 . 

As discussed in section 7.3, partons with sufficiently large x p , satisfying Eq. (88), are expected to 
interact with small x partons in the nucleus with a probability of order one. As a result, a hard collision 
of partons with sufficiently large x p and xa > 0.01 will be accompanied by production of one or more 
minijets in the black-body regime with px ~ p^(xa)- Most sufficiently fast partons from the proton 
will generate minijets, leading to a strong suppression of the cross section of events with only two jets. 

Detectors with sufficiently forward rapidity acceptance could detect events originating from triple 
parton collisions at xa > £ s h an d large enough pt where pr broadening related to the blackbody limit 
effects are sufficiently small. These triple parton collisions provide stringent tests of the hard interaction 
dynamics as well as provide additional information on two-parton correlations and unique information 
on triple-parton correlations. 

Other opportunities with multi-jets include: 

• Probing correlations between partons in the nucleus at high densities, i.e. for x±a,X2A *C 10~ 3 , 
which would provide qualitatively new information about the dynamics of nuclear shadowing and 
the presence of possible new parton condensates. 

• Studying the accompanying soft hadron production (cf. the discussion in section 7.42), which 
would measure the transverse size of a proton configuration containing partons x\ and xi- In 
particular, the production of 4 jets with two of the jets at large x p , x\ + X2 > 0.5, may provide 
another way to look for point-like configurations in nucleons, see section 7.5 and Ref. [325] for 
discussion of other possibilities. 

7.42 Proton-ion collisions probe transverse nucleon structure 

As we discuss in section 7.5, in some subclass of events the distribution of constituents in the initial 
proton may be unusually local in the transverse (impact parameter) plane when the proton collides with 
the ion. If this is so, its effective cross section per nucleon will be greatly reduced, perhaps all the way to 
the perturbative-QCD level. If the effective cross section of such a point-like configuration goes below 
20 mb, there will be an appreciable probability that it can penetrate through the center of a lead ion and 
survive, leading to a strongly enhanced diffractive yield of the products of the point-like configuration in 
collisions with heavy ions. 

Not only might the properties of the final-state collision products depend upon the nature of the 
transverse structure of the proton primary on arrival at the collision point, but even the conventional 
parton distributions may also be affected. For example, let us consider a nucleon as a quark and small 
diquark connected by a narrow QCD flux-tube. It should be clear that if the flux-tube is at right angles 
to the collision axis at arrival, then the valence partons will have comparable longitudinal momentum 
or x. On the other hand, if the flux tube is parallel to the direction of motion, then one of the valence 
systems will have very large x and the other very small. This happens, because in this case, the internal 
longitudinal momenta of the valence systems, in the rest frame of the projectile proton, are in opposite 
directions. Therefore the smallness of the configuration is correlated with the joint x-distribution of 
its constituents. This kind of nonfactorization may be determined by the study of perturbative QCD 
processes such as dilepton, direct-photon, or dijet production as a function of the centrality 15 , soft hadron 
multiplicity, and A. Naturally such studies would also require investigation of soft hadron production as 
a function of impact parameter. At LHC energies it may differ quite strongly from fixed target energies. 

One possible kinematics where a strong correlation is expected is when a parton with large x > 
0.6, is selected in the proton. The presence of such a parton requires three quarks to exchange rather 

14 The A-dependence of the ratio of a% ja\, in the kinematics where only one of the nuclear partons has xa < x B h is 
practically the same as for the case when both nuclear partons have x > x^. 

15 The impact parameter dependence is accessible in the pA collisions via a study of nuclear fragmentation into various 
channels. For an extensive discussion see Ref. [325]. 



large momenta. Hence one may expect that these configurations have a smaller transverse size and 
therefore interact with the target with a smaller effective cross section, a e s(x). Suggestions for such x 
dependence on the size are widely discussed in the literature. Using a geometric (eikonal type) picture of 
pA interactions as a guide and neglecting (for simplicity) shadowing effects on nPDFs, one can estimate 
the number of wounded nucleons, u(x, A), in events with a hard trigger (Drell-Yan pair, 7-jet, dijet...) as 
a function of cr e fr [340] : 



where the nuclear density per unit area T^{B) is defined in Eq. (92). At the LHC, for average inelastic 
pPb collisions and cj c fj ~ cr mc \(pp), Eq. (93) leads to u m 10, somewhat larger than the average number 
of wounded nucleons in pA collisions due to the selection of more central impact parameters in events 
with a hard trigger. 

A decrease of the effective cross section for large x, say, by a factor of 2, would result in a 
comparable drop of the number of particles produced at central rapidities as well as in a smaller number 
of nucleons produced in the nuclear fragmentation region. 

7.43 A-dependence of the particle production 

Proton fragmentation region. The A-dependence of hadron production in the proton fragmentation 
region remains one of the least understood aspects of hadron-nucleus interactions. Practically all avail- 
able data are inclusive and correspond to energies where the inelastic NN cross section is about three 
times smaller than at the LHC. They indicate that the cross section is dominated by the production of 
leading particles at large impact parameters where the projectile interacts with only one or two nucleons 
of the nucleus. As a result, very little information is available about hadron production at the central 
impact parameters most crucial for the study of AA collisions. Theoretical predictions for this region are 
also rather uncertain. 

In eikonal-type models, where the energy is split between several soft interactions, one may expect 
a very strong decrease of the leading particle yield. The dependence is expected to be exponential with 
path length through the nucleus: the mean energy is attenuated exponentially. On the other hand, if the 
valence partons of the projectile do not lose a significant amount of their initial momentum, as suggested 
by pQCD motivated models, see the review in Ref. [187], the spectrum of leading particles may approach 
a finite limit for large A and central impact parameters [341]. Indeed, in this case, the leading partons will 
acquire significant transverse momenta and will not be able to coalescence back into leading baryons and 
mesons, as seemingly happens in nucleon-nucleon collisions. As a result, they will fragment practically 
independently, leading to much softer longitudinal momentum distributions for mesons and especially 
for baryons [341]. 

We discussed previously that the interaction strength of fast partons in the nucleon (x p > 10~ 2 for 
LHC) with heavy nuclei at central impact parameters may approach the black-body limit. In this scenario, 
all these partons will acquire large transverse momenta ~ p^ 1 (xa), 16 leading to a very significant pt 
broadening of the leading hadron spectrum [120]. Probably the most feasible way to study this effect is 
measurement of the pr distribution of leading hadrons (neutrons, Kl,tt°) in central collisions. These 
hadrons will closely follow the pt of the leading quarks [343]. The results of the calculation [343] 
using the color glass condensate model for modeling p? broadening are shown in Fig. 53 for different 
values of the saturation parameter Q s . It is worth emphasizing that strong pr broadening, increasing 
with xp, clearly distinguishes this mechanism of leading hadron suppression from the soft physics effect 
of increasing the total pp cross section by nearly a factor of three over fixed target energies. 

16 Hence in the black-body limit, pr broadening of the partons should be much larger than at low energies where it is 
consistent with QCD multiple rescattering model [187]. See Ref. [342] for a discussion of matching these two regimes. 
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Fig. 53: Left: Transverse momentum distribution of neutrons (in fact, n/2 + n/2) from the breakup of an incident proton at 
various longitudinal momentum fractions z and target saturation momenta Q s (bottom six curves). The top three curves depict 
the underlying quark distributions. Right: Longitudinal distributions of n/2 + n/2 and 7r°, integrated over kr- From Ref. 
[343]. 

Central region. The Gribov-Glauber model, including inelastic screening effects but neglecting the 
final state interactions between the hadrons produced in different Pomeron exchanges, as well as ne- 
glecting the interactions of these hadrons with the nuclear target, predicts that, at high energies and 
rapidities y such that \y p — y\ 3> 1, the inclusive pA spectrum should be proportional to that of pN 
scattering: 

da^ A ^ x (y,p T ) = A daP+ N - h + x (y,p T ) 
dy d 2 px dy d 2 px 

Since (j- ine \(pA) / <Ji ne \{pN) oc A 2 / 3 , Eq. (94) implies that the multiplicity of particles produced in the 
inelastic pA collisions should increase as oc A 1 / 3 . Data at fixed target energies do not contradict this 
relation but the energy is too low for an unambiguous interpretation. Effects of splitting the energy be- 
tween multiple soft NN interactions are much larger in AA collisions than in pA collisions so checking 
whether a similar formula is valid for AA collisions would require energies much higher than those 
available at RHIC. 

The opposite extreme is to assume that parton interactions in the black body regime become im- 
portant in a large range of x p > 10~ 2 . In this case, multiple scatterings between the partons, their in- 
dependent fragmentation together with associated QCD radiation (gluon bremsstrahlung), may produce 
most of the total entropy and transverse energy and will enhance hadron production at central rapidities 
relative to the Gribov-Glauber model. 



Nuclear fragmentation region. For the rapidities close to that of the nucleus, there are indications of 
slow hadron rescatterings, leading to increased nucleon multiplicity forl.0>pAr>0.3 GeV compared 
to no final state reinteractions, a factor of ~ 2 for heavy nuclei [344]. As discussed above, it is doubtful 
that the Gribov-Glauber picture is correct at LHC energies, so that significant deviations from Eq. (94) 
maybe expected. It is likely that a colored quark-gluon system with rapidities close to yA will be excited 
along a cylinder of radius ~ 1 fm and length 2Ra, resulting in unusual properties of hadron production 
in the nuclear fragmentation region, including production of exotic multi-quark, multi-gluon states. 

We emphasize that a quantitative understanding of particle production in pA collisions in this 
phase-space region is a prerequisite for understanding the corresponding phenomenology in ^4^4 colli- 
sions. 



7.5 A new hard QCD phenomenon: proton diffraction into three jets 

7.51 Introduction 

During the last ten years, a number of new hard small-x phenomena have been observed and calculated 
in QCD based on the QCD factorization theorems for inclusive DIS and hard exclusive processes. These 
phenomena include: (i) observation of a fast increase of the PDFs with energy at HERA, reasonably 
described within the QCD evolution equation approximation 17 ; (ii) discovery of color transparency in 
pion coherent dissociation into two high-p^ jets [345], consistent with the predictions of Refs. [346, 347]; 
(iii) observation of various regularities in exclusive vector-meson electroproduction at HERA induced 
by longitudinally polarized photons and photoproduction of mesons with hidden heavy flavor consistent 
with the predictions of Ref. [348]. These processes, in a wide kinematic range of small x, have provided 
effective ways to study the interaction of small colorless dipoles with hadrons at high energies and to 
study the hadron wave functions in their minimal Fock space configurations. Here we outline new QCD 
phenomena [325] which may become important at the LHC because their cross sections are rapidly 
increasing with energy. 

A fast increase of hard diffractive cross sections with energy, predicted within the leading-twist 
(LT) approximation, would be at variance with unitarity of the S'-matrix for the interaction of spatially 
small wave packet of quarks and gluons with a hadron target at a given impact parameter. This new 
phenomenon may reveal itself at sufficiently small x and/or sufficiently heavy nuclei. In the proton 
structure functions, this physics is masked by a significant contribution to the hard collisions of the 
nucleon periphery. Thus, unitarity of the S'-matrix does not preclude an increase of nucleon structure 
functions as fast as oc ln(l/x) 3 [331]. Hence no dramatic signals for the breakdown of the DGLAP 
regime are expected in this case. On the contrary, a breakdown of the LT approximation in hard diffractive 
processes off nuclei leads to drastic changes in the cross section dependence on A, incident energy and 
jet p T . 

7.52 Three-jet exclusive diffraction 

A nucleon (meson) has a significant amplitude to be in a configuration where the valence partons are 
localized in a small transverse area together with the rest of the partons. These configurations are usually 
referred to as minimal Fock space configurations - 1 3q) . Experimental evidence for a significant prob- 
ability for these configurations in mesons include the significant decay amplitude for it — ► leptons as 
well as the significant observed diffractive vector meson electroproduction cross sections and the process 
7r + A — > 2 jets + A [345]. In the proton case, the amplitude of the three quark configuration can be 
estimated within QCD-inspired models and, eventually calculated in lattice QCD. Note that knowledge 
of this amplitude is important for the unambiguous calculation of proton decay within Grand Unification 
models. 

Hadrons, in small size configurations with transverse size d, interact with a small cross section 
oc c1 2 xGn(x,Q 2 ), where Q 2 ~ W/d 2 . The factor xGn(x,Q 2 ) leads to a fast increase of the small 
cross section with energy. In the case of "elastic scattering" of such a proton configuration off another 
proton, this three-quark system with large relative momenta should preferentially dissociate diffractively 
into a system of three jets with large transverse momenta pxi ~ n/d, where d is the transverse size of 
the minimal configuration. The kinematics of this process is presented in the LEGO plot of Fig. 54 for 
the case of coherent scattering off a nucleus. 

The 3-jet production cross section can be evaluated based on the kind of QCD factorization the- 
orem deduced in Ref. [347]. The cross section is proportional to the square of the gluon density in the 
nucleon at x ps M 2 - ]Cts /s and virtuality ~ (1 — 2)p\ [349]. The distribution over the fractions of the 
longitudinal momentum carried by the jets is proportional to the square of the light cone wave function 

I7 For a comprehensive review of the HERA data on inclusive and exclusive processes and the relevant theory, see Ref. [335]. 




Fig. 54: LEGO plot for diffraction of proton into three jets. 



of the \3q) configuration, ip(zi, z 2 , z 3 ,pti,PT2,PT3) and large transverse momenta where 
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(95) 



The proper renormalization procedure accounting for cancellation of infrared divergences is im- 
plied. Otherwise the light cone gauge p^^A^ = 0, should be used to describe the diffractive fragmenta- 
tion region where there are no infrared divergences. Here p^ is the momentum of the diffracted proton. 
The summation over collinear radiation is included in the definition of the jets. 

Naively, one can hardly expect that exclusive jet production would be a leading-twist effect at 
high energies since high-energy hadrons contain a very large number of soft partons in addition to three 
valence quarks. However application of the QCD factorization theorem demonstrates that in high-energy 
exclusive processes these partons are hidden in the structure function of the target. Hence proton diffrac- 
tion into three jets provides important information about the short-distance quark structure of the proton 
and also provides unique information about the longitudinal momentum distribution in the 1 3q) configu- 
ration at high pt- The amplitude of the process can be written as 



A 



-Alf)(zi,Z 2 ,Z 3 ,pTl,PT2,PT3) 



2tt< 



a s (xG(x,Q 2 ,t)) 



(96) 



Here A = £? =1 f^. In Eq. (96) we simplify the expression by substituting a skewed gluon density 
reflecting the difference in the light-cone fractions of the gluons attached to the nucleon by the gluon 
density which is numerically a reasonable approximation for rough estimates. An additional factor of 2 
in the numerator relative to the pion allows two gluons to be attached to different pairs of quarks. The 
Pt and longitudinal momentum fraction, Zi, dependence of the cross section is given by 



da 
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(97) 



The coefficient cn is calculable in QCD and 4>n(z\, z 2 , z 3 ) = z\z 2 z 3 . The two-gluon form factor of 
a nucleon is known to some extent from hard diffractive processes [351]. A numerical estimate of the 



three-jet production cross section with p? of one of the jet greater than a given value and integrated over 
all other variables gives, for px > 10 GeV/c, 



. a s xG(x,Q 2 2 fi 7 /10GeV\ 8 

a(pp -► 3 jets + p) oc — V « !0 • • • 10~ 7 mb (98) 

Pt V Pt J 

for LHC energies. An urgent question is whether this cross section will continue to grow up to LHC 
energies, as assumed here based on extrapolations of Gn(x,Q 2 ) to x ~ 1CT 5 . The probability of 
the 1 3q) configuration is estimated using a phenomenological fit to the probability of configurations of 
different interaction strengths in a nucleon, cf. Refs. [349, 350]. 

A study of the same process in pA collisions would provide an unambiguous test of the dominance 
of hard physics in this process. The cross section should grow with A as -^k-G^x, Q 2 ) / A 2 G 2 N {x , Q 2 ). 

A 

Accounting for the skewedness of the gluon distribution will lead to some enhancement in the absolute 
value of the cross section but will not strongly influence the A-dependence. 

Our estimates indicate that for the relevant LHC kinematics, x ~ 10~ 5 and Q 2 S ~ p\, nuclear 
shadowing in the gluon channel is a rather small correction, reducing the A-dependence of the cross 
section by A 0A , leading to a(pA — > 3jets + A) oc A 1 - 25 . The A-dependence of total inelastic cross 
section is oc A - 7 . Hence we expect that the counting rate per inelastic interaction will be enhanced 
in pA collisions for heavy nuclear target by a factor oc A ' 55 . At the same time, the background from 
soft diffraction should be strongly reduced. Indeed, if we estimate the A-dependence of soft diffraction 
based on the picture of color fluctuations, which provides a good description of total cross section of 
the inelastic diffraction off nuclei at fixed target energies [327], we find cr^f 1 ' dlff - oc A - 25 . Scattering 
off nuclei also has an obvious advantage in terms of selecting coherent processes since practically all 
inelastic interactions with nuclei would lead to the emission of neutrons at 0°. 

A competing process is proton diffraction into 3 jets off the nuclear Coulomb field. It follows from 
the theorem proved in Refs. [347, 352] that to leading order in a s the dominant contribution is given in 
the Weizsacker- Williams representation by transverse photon interactions with external quark lines. The 
amplitude of the process is [352] : 

where is the electric charge of the quark in the units of the electron charge, q 2 is the virtuality of the 
exchanged photon and D(p T ) is the renormalization factor for the quark Green's function. 

In the kinematics where LT dominates the QCD mechanism for 3-jet production, the electromag- 
netic process is a small correction. However if the screening effects (approach to the black body limit) 
occurs in the relevant pr range, the electromagnetic process may compete with the QCD contribution. 
In the case of nuclear breakup, the electromagnetic contribution remains a relatively modest correction, 
even in the black body limit. This is primarily because the cross section for electromagnetic nuclear 
breakup does not have a singularity at q 2 = 0. Note also that, in the case of QCD interactions, one can 
consider processes with a large rapidity gap which have the same Z{ and pn dependence as the coherent 
process but with a significantly larger cross section. 

Another interesting group of hard processes is proton diffraction into two high pt jets and one 
collinear jet. These processes are, in general, dominated at high energies by collisions of two nucleons 
at large impact parameters (otherwise multiple soft interactions would destroy the coherence). Hence the 
A-dependence of the dominant term should be similar to that of cr^f 1 ' . However fluctuations with 
cross sections of ~ 40 mb would have a much stronger A-dependence, ~ A - 7 . Thus the study of the 
A-dependence of coherent diffraction can filter out the smaller than average components in the nucleons 
at LHC energies. 



7.53 Violation of the LT approximation and new QCD strong interactions 

The energy dependence of cross section in the LT approximation follows from the properties of the tar- 
get gluon distribution. Hence for the virtualities probed in 3-jet production, the energy dependence is 
oc x 2 G 2 {x, Q 2 cS ) oc s 8 for Q 2 cS ~ 100 - 1000 GeV 2 and 1(T 5 < x < 10~ 2 . In the black-body limit, 
the inelastic diffractive processes are strongly suppressed and hence the energy dependence of the 3-jet 
production cross section integrated over pr > PTmm ~a few GeV should slow down since only scatter- 
ing off the grey edge of the nucleus (nucleon) will effectively survive. For Q 2 relevant to this process, 
the black-body limit may only be an interesting hypothesis for scattering off heavy nuclei. It would man- 
ifest itself in a strong decrease of the ^-dependence of 3-jet production at the pr corresponding to the 
black-body limit. At sufficiently small x for a heavy nuclear target, da(pA — > 3jets + A)jdp\ oc A 1 / 3 , 
in striking contrast with the LT expectations, oc A 1 - 2 . 

7.6 Conclusions 

In conclusion, studies of pA collisions at the LHC will reveal strong interactions in the high field domain 
over an extended rapidity region. New phenomena will be especially prominent in the proton fragmen- 
tation region but will extend to the central region and the nuclear fragmentation region as well. 
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